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Whenever  moist  air  exists  in  the  vicinity  of  a surface  whose  temperature  is  both 
below  the  dew-point  temperature  of  water  vapor  in  air  and  below  the  freezing  point,  frost 
will  form.  This  frost  contributes  to  coil  heat-transfer  performance  degradation  due  to 
both  the  insulating  effect  of  the  frost  layer  and  the  resulting  coil  blockage.  While  the 
frosting  problem  has  been  extensively  researched,  no  study  has  investigated  the  problem 
when  the  freezer  is  operating  near  saturated  conditions.  Field  observations  of  industrial 
freezer  operation  indicate  that  this  mode  of  operation  is  much  more  common  than 
originally  thought;  and  that  the  transition  to  supersaturated  operation  can  easily  occur  if 
proper  care  is  not  exercised  by  either  the  refrigeration  system  designer  or  the  freezer 
operator.  Field  observations  also  reveal  that  operating  freezers  in  supersaturated 
conditions  can  significantly  accelerate  the  occurrence  of  the  negative  aspects  of  frost 
formation  vis-a-vis  coil  heat  transfer  performance.  For  this  and  other  reasons,  the  study 
presented  in  this  dissertation  focused  on  that  aspect  of  freezer  operation  with  the  hope  of 
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developing  tools  and  protocols  that  may  help  guide  refrigeration-system  designers  and 
freezer  operators  avoid  the  conditions  that  contribute  to  this  mode  of  operation;  and  thus 
help  improve  the  performance  of  freezer  coils  that  are  particularly  prone  to  the  frosting 
problem. 
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CHAPTER  1 
INTRODUCTION 

Literature  Review 

Commercial  freezers  of  the  type  found  commonly  in  food  distribution  warehouses 
and  supermarkets  are  prone  to  frost  formation  inside  and  around  the  freezer  door, 
especially  if  the  freezer  has  a high  volume  of  traffic.  The  problem  is  the  exchange  of 
warm  and  cold  air  that  takes  place  when  the  freezer  door  is  opened  to  deliver  or  retrieve 
goods.  Basically,  cold  air  from  the  freezer  rushes  out  along  the  floor  of  the  vestibule  or 
anteroom,  while  warm  and  humid  outside  air  rushes  inside  the  freezer  near  the  ceiling. 
The  moisture  introduced  into  the  freezer  space  is  likely  to  cause  the  air  to  become 
supersaturated1  at  the  prevailing  air  temperature.  If  that  temperature  is  above  the  freezing 
point  of  water,  moisture  in  the  air  would  exist  in  the  form  of  tiny  liquid  droplets 
suspended  in  the  air  stream — a familiar  condition  known  as  fog.  If  the  same  scenario 
exists,  but  with  the  air  temperature  below  the  freezing  point,  ice  fog  forms.  Air-borne  ice 
crystals  that  normally  constitute  ice  fog  have  great  affinity  for  cold  surfaces;  and  thus 
tend  to  deposit  on  any  surface  they  can  find  inside  the  freezer. 

While  frost  formation  on  cooling  coils  in  domestic  refrigeration,  air  conditioning, 

and  heat  pump  applications  has  been  recognized  as  a problem  that  contributes  to  heat 

transfer  coil  performance  degradation  (Sherif  1991),  the  problem  in  industrial  walk-in 

freezers  is  much  more  significant  and  potentially  more  expensive  to  deal  with. 

1 In  this  investigation  the  term  supersaturated  air  is  defined  as  a mixture  of  saturated  air  and  suspended 
liquid  water  droplets  or  ice  crystals  depending  on  whether  the  prevailing  air  temperature  is  larger  or  smaller 
than  the  freezing  point  temperature  of  water,  respectively. 


1 


2 


Reasons  that  contribute  to  this  fact  are  many  and  include  the  following: 

• The  temperature  difference  between  the  freezer  and  the  anteroom  is  larger  than  in  air 
conditioning  applications. 

• The  temperature  of  the  freezer  is  significantly  lower  than  those  temperatures 
encountered  in  air  conditioning  and  heat  pump  applications. 

• The  magnitude  of  potential  energy  waste  in  industrial  walk-in  freezers  is  significantly 
larger  than  the  corresponding  problem  in  either  domestic  refrigeration  or  heat  pump 
and  air  conditioning  applications. 

The  combination  of  the  first  two  facts,  coupled  with  the  manner  in  which  the 
saturation  curve  on  the  psychrometric  chart  is  shaped  in  the  low  temperature  region, 
contribute  to  increasing  the  likelihood  that  the  freezer  would  operate  under  supersaturated 
conditions  and  that  airborne  ice  crystals  would  form. 

While  literature  dealing  with  the  frosting  problem  on  different  geometries  is 
abundant,  our  review  is  confined  to  finned  coils,  cylindrical  geometries,  and  studies 
dealing  with  frost  property  formulations.  Among  the  investigators  who  reported  frost 
formation  on  cylinders  are  Andrichak  (1962),  Chung  and  Algren  (1958a,  1958b),  Padki  et 
al.  (1989),  Parish  and  Sepsy  (1972),  Raju  and  Sherif  (1993),  and  Sengupta  et  al.  (1998). 
Chung  and  Algren  (1958a,  1958b)  reported  experimental  mass  transfer  results  for 
cylinders  in  cross  flow  of  humid  air  for  a host  of  environmental  and  surface  conditions, 
represented  by  the  free  stream  Reynolds  number,  the  ambient  dry-bulb  temperature,  the 
humidity  ratio,  and  the  cylinder  surface  temperature.  The  results  presented  were  in  terms 
of  the  frost  total  mass  accumulated  on  the  cylinder  as  a function  of  time  for  a range  of 
Reynolds  numbers  from  2592  to  9965,  ambient  dry-bulb  temperatures  of  15°C  to  21.4°C, 
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and  cylinder  surface  temperatures  of  -1 1.9°C  to  -23°C.  In  a similar  manner,  Andrichak 
(1962)  experimentally  studied  the  frost  formation  on  a cylinder  in  cross  flow  for  different 
angular  positions  from  the  stagnation  point  to  the  separation  point,  which  was  considered 
to  be  80°  for  laminar  flow  over  a cylinder.  He  found  that  the  frost  thickness  decreases 
with  the  angular  position  on  the  cylinder.  Sengupta  et  al.  (1998)  performed  an 
experimental  investigation  in  which  the  emphasis  was  placed  on  obtaining  correlations 
for  the  thickness-time  history  and  the  heat  transfer-time  history  for  a cylinder  in  cross 
flow  of  humid  air.  They  used  a low-velocity  wind  tunnel  comprising  a rectangular  test 
section,  a transition  section,  and  a honeycomb  placed  at  the  tunnel  entrance.  Humid  air 
was  forced  over  the  coil  using  a suction  fan  located  at  the  downstream  end  of  the  tunnel. 
They  derived  correlations  in  terms  of  nondimensional  parameters  representing  the  free- 
stream  velocity,  humidity  ratio,  and  temperature,  as  well  as  time.  They  found  these 
correlations  to  agree  well  with  the  experimental  data.  On  the  other  hand,  Parish  and 
Sepsy  (1972)  and  Raju  and  Sherif  (1993)  derived  a transient  model  to  predict  the  frosting 
problem  over  the  same  geometry  as  Andrichak  (1962).  They  used  a numerical  solution 
for  the  momentum,  energy,  and  diffusion  boundary-layer  equations  along  with  the 
continuity  equation,  using  a finite  difference  scheme.  Both  groups  of  researchers 
reported  the  variation  of  frost  thickness  with  time  at  two  specific  locations  on  the 
cylinder:  the  stagnation  point  and  the  separation  point  (that  was  considered  to  be  80°C). 
Raju  and  Sherif  (1993)  also  reported  the  variation  of  the  frost  surface  temperature  with 
the  angular  position  along  the  circumference  of  the  cylinder  for  the  same  range 


mentioned  above  (0°C  to  80°C). 
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Among  the  investigators  who  studied  finned  coils  are  Chen  et  al.  (2000a,  b) 
Kondepudi  and  O’Neal  (1988,  1989,  1991),  and  Tao  et  al.  (1994).  Kondepudi  and 
O’Neal  (1988,  1989,  1991)  studied  the  effects  of  frost  growth  on  the  thermal  performance 
of  finned  tube  heat  exchangers  with  triangular  spine  fins,  with  flat  fins,  and  with  louvered 
fins,  respectively.  In  their  results  the  frost  accumulation,  pressure  drop  across  the  coils, 
and  energy  transfer  coefficient  were  quantified  with  respect  to  air  humidity,  refrigerant 
temperature,  and  fin  geometry.  They  found  that  higher  air  humidities,  lower  refrigerant 
temperatures,  and  higher  fin  densities  lead  to  increased  frost  accumulation,  higher 
pressure  drops,  and  higher  energy  transfer.  They  also  found  that  the  louvered  fins  had 
more  frost  accumulation  and  air-side  pressure  drop  across  the  coils  compared  to  the  flat 
fins. 

Tao  et  al.  (1994)  summarized  some  of  the  studies  performed  in  this  field  and,  in 
particular,  reported  on  measurements  for  frost  growth  during  the  initial  stage  on  the 
surfaces  of  different  materials.  They  also  reported  on  new  numerical  results  of  frost 
formation  in  the  freezer  temperature  range. 

Chen  et  al.  (2000a,  b)  presented  a numerical  model  for  predicting  frost  growth  on 
plate-fin  heat  exchangers.  The  model  they  presented  treats  the  frost  growth  as  a porous 
media  problem  that  is  transient  and  one-dimensional  and  includes  transient,  two- 
dimensional  heat  conduction  in  the  fins.  They  found  that,  for  a constant  air  flow  rate 
through  the  heat  exchanger,  the  blockage  of  air  flow  resulting  from  the  frost  growth 
increases  the  pressure  drop  across  the  heat  exchanger  by  up  to  a factor  of  8,  while  the 
heat  transfer  rate  decreases  by  20%.  In  another  investigation,  the  same  authors  compared 
the  numerical  model  results  with  experimental  data  for  frost  height,  frost  density,  airflow 
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pressure  drop  through  the  finned  section,  and  heat  transfer  rate  through  the  cold  base, 
finding  that  the  simulation  results  and  experimental  data  agree  within  the  experimental 
uncertainty  with  few  exceptions. 

Studies  emphasizing  empirical  or  theoretical  frost  property  formulations  were 
reported  by  Biguria  and  Wenzel  (1970),  Brian  et  al.  (1969),  Hayashi  et  al.  (1977a, 
1977b),  Marinyuk  (1980),  Yonko  and  Sepsy  (1967),  among  others.  Yonko  and  Sepsy 
(1967)  performed  an  experimental  investigation  and  proposed  an  empirical  correlation  for 
the  frost  thermal  conductivity  based  on  the  frost  density.  Their  correlation  is  valid  for 
frost  densities  less  than  573  kg/m3.  In  a similar  manner,  Brian  et  al.  (1969)  proposed  an 
empirical  correlation  for  the  frost  thermal  conductivity  as  a function  of  the  frost  density 
and  average  frost  temperature.  His  correlation  was  developed  for  a small  range  of  frost 
densities  (25  to  130  kg/m3),  which  limited  its  applicability.  On  the  other  hand,  Biguria 
and  Wenzel  (1970)  presented  equations  to  determine  the  frost  density  and  frost  thermal 
conductivity  based  on  the  humidity  of  the  free  stream,  the  velocity  of  the  air,  and  the 
temperatures  of  the  frost  and  the  plate  surface.  Their  correlation  provided  results  for  frost 
density  and  thermal  conductivity  with  an  error  of  ±7%  and  ±25%,  respectively.  Hayashi 
et  al.  (1977a,  1977b)  classified  the  frost-formation  process  into  three  discrete  periods: 
crystal  growth,  frost  layer  development,  and  frost  layer  full  growth.  They  studied  the 
frost  crystal  growth,  density,  and  thermal  conductivity  in  each  of  the  periods.  From  their 
studies,  they  proposed  an  empirical  correlation  for  the  frost  density  as  a function  of  the 
frost  surface  temperature,  for  airstream  velocities  between  2 m/s  and  6 m/s,  frost  surface 
temperature  between  -25°C  and  0°C,  and  an  airstream  humidity  ratio  equal  to  0.0075 


kgw/kga. 
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Marinyuk  (1980)  performed  an  experimental  investigation  from  which  a 
correlation  for  the  frost  thermal  conductivity  was  proposed,  based  on  frost  density,  frost 
surface  temperature,  and  the  wall  temperature.  His  correlation  was  limited  to  frost 
densities  up  to  300  kg/m  . 

Despite  the  large  effort  expended  on  investigating  the  problem  of  frost  formation, 
very  few  studies  looked  into  the  formation  mechanism  under  supersaturated  conditions 
while  accounting  for  psychrometric  effects.  Thus,  it  appears  that  any  new  effort  to  study 
the  frosting  problem  should  try  to  address  the  aforementioned  aspect  of  the  issue. 

While  there  is  a significant  amount  of  information  pertaining  to  coil  frosting, 
studies  dealing  with  coil  defrosting  are  relatively  few.  Among  the  investigators  who 
examined  defrosting  of  freezer  coils  are  Coley  (1983),  Kerschbaumer  (1971),  Niederer 
(1976),  and  Stoecker  et  al.  (1983).  Kerschbaumer  (1971)  defined  the  defrost  efficiency 
as  the  ratio  of  the  total  heat  required  to  melt  the  frost  alone  to  the  total  amount  of  defrost 
heat  input.  Niederer  (1976)  attempted  to  determine  the  amount  of  heat  required  for 
defrost  by  measuring  the  amount  of  hot-gas  condensing  and  the  power  input  during  the 
defrost  cycle.  He  found  that  only  15  to  25%  of  the  total  heat  went  to  melt  the  frost.  He 
observed  that  the  remaining  75  to  85%  of  the  heat  went  to  heat  the  surrounding 
environment  and  the  coil  surfaces.  Stoecker  et  al.  (1983)  performed  a study  to  analyze 
hot-gas  defrosting  of  coils  in  industrial  refrigeration  systems.  Coley  (1983)  suggested 
that  during  a defrost  cycle  in  freezers,  at  least  15%  of  the  frost  sublimates  into  the 
surrounding  refrigerated  air  space.  He  noted  that  this  represented  an  additional  amount 
of  load  to  be  subsequently  removed  during  the  refrigeration  portion  of  the 
refrigeration/defrost  (R/D)  cycle. 
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While  there  is  a reasonable  amount  of  interest  in  coil  defrosting  studies,  it  is 
obvious  that  there  are  areas  that  can  benefit  from  additional  research.  For  example,  there 
seems  to  be  no  quantitative  information  on  the  performance  of  freezer  coils  at  low 
temperatures  during  the  defrosting  process  near  saturated  freezer  air  conditions. 
Developing  a reliable  and  extensive  body  of  data  for  those  types  of  scenarios  should 
therefore  prove  a worthy  cause. 

Objectives  and  Scope 

The  main  objective  of  this  investigation  is  to  develop  tools  and  protocols  that  may 
help  refrigeration  system  designers  and  freezer  operators  minimize  the  negative  impact  of 
the  coil  frosting  problem  in  industrial  refrigeration  applications.  Since  the  frosting 
problem  has  been  extensively  researched,  our  emphasis  was  on  freezer  operation  near 
saturated  conditions  (a  mode  of  operation  that  has  never  been  studied  before  in  the 
freezer  temperature  range).  To  achieve  this  objective,  the  following  milestones  had  to  be 
accomplished: 

• Establish  if  there  is  a demarcation  line  between  the  type  of  frost  that  forms  under 
subsaturated  conditions,  and  the  type  that  forms  under  supersaturated  conditions. 
Once  the  existence  of  the  demarcation  line  was  established,  its  location  was 
determined.  This  is  important  since  we  need  to  know  where  the  transition  from  one 
condition  to  the  other  occurs  since  the  nature  of  frost  forming  under  both  modes  is 
different  as  suggested  by  field  observations. 

• Investigate  the  issues  pertaining  to  the  coil  defrosting  process  using  the  hot-gas 
refrigerant  method.  This  method  of  defrosting  was  selected  because  it  is  the  most 
widely  used  defrosting  technique  in  industrial  freezers.  This  part  of  our  study 
included  both  physical  and  infrared  analysis  of  a coil  being  defrosted.  The  infrared 
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portion  of  the  analysis  provided  the  spatial  and  temporal  distributions  of  temperature 
on  the  coil  surface,  a key  variable  in  understanding  the  defrosting  dynamics. 
Throughout  this  study,  we  also  attempted  to  compare  the  performance  of  the  defrosted 
coil  with  the  performance  of  similar  coils  operating  under  subsaturated  conditions. 

• Study  alternative  energy  conservation  technologies  in  industrial  freezers  that  are  prone 
to  coil  frosting  problems.  The  feasibility  of  three  distinct  technologies  was  explored 
both  experimentally  and  analytically.  These  techniques  are:  the  use  of  heat-pipe 
assisted  heat  exchangers,  use  of  dampered  coils,  and  operating  freezer  coils  at  high 
coil-face  velocities.  These  technologies  were  selected  while  pursuing  the 
experimental  investigation  described  in  this  dissertation.  Their  selection  was  primarily 
based  on  observations  made  while  exploring  how  freezer  coils  perform  near  saturated 
conditions.  While  these  techniques  were  not  initially  identified  as  potentially 
beneficial  technologies,  experimental  observations  suggested  a potential  benefit,  and 
hence  it  was  decided  to  explore  their  merit. 

• Demonstrate  the  development  and  use  of  a calculation  procedure  for  coil  surfaces 
exposed  to  supersaturated  air.  While  every  effort  should  be  made  to  avoid  this  mode 
of  operation,  it  may  be  advantageous  to  develop  a computational  capability  to  predict 
frost  formation  and  heat  transfer  rates  on  such  a coil.  For  demonstration  purposes,  a 
cylinder  in  humid  air  cross  flow  was  used.  The  cylinder  surface  was  selected  because 
it  could  represent  the  tube  in  a freezer  coil.  While  this  is  not  identical  to  a finned-tube 
freezer  coil  per  se,  it  shows  how  the  techniques  developed  in  this  study  can  be  used  to 
predict  the  performance  of  a real  freezer  coil  in  the  supersaturated  zone  of  the 
psychrometric  chart.  This  geometry  has  also  been  selected  so  that  results  of  the 
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supersaturated  model  could  be  compared  to  other  published  works  in  the  subsaturated 
zone  of  the  psychrometric  chart.  To  do  this,  psychrometric  formulations  and 
thermodynamic  properties  for  supersaturated  air  were  developed.  Achieving  this 
objective  could  potentially  break  new  ground  in  the  area  of  coil  frosting  as  the  lack  of 
psychrometric  data  in  the  supersaturated  zone  of  the  psychrometric  chart  has 
historically  prevented  such  a capability  from  being  developed.  These  properties  will 
expand  the  available  psychrometric  data  as  well  as  allow  us  to  develop  new 
psychrometric  charts  incorporating  moist  air  properties  in  the  supersaturated  zone. 

• The  final  objective  is  geared  toward  practicing  refrigeration  engineers  and  freezer 
operators.  The  intent  here  is  to  develop  two  pieces  of  software  that  can  be  used  in  the 
design  and  operational  phases  of  refrigerated  facilities.  The  first  piece  of  software  was 
based  on  the  analytical  model  described  in  the  previous  objective.  This  software  was 
used  to  determine  the  frost-formation  and  heat-transfer  rates  on  a cylinder  surface 
under  forced  convection.  It  is  important  to  point  out  that  although  the  results  obtained 
using  the  software  are  for  a cylinder,  they  are  a good  first  approximation  for  tubes  in  a 
freezer  coil.  This  software  will  perform  analysis  under  both  subsaturated  and 
supersaturated  freezer  conditions.  The  other  piece  of  software  can  be  used  as  a tool  to 
prevent  freezers  from  operating  under  supersaturated  conditions.  This  would  allow 
operators  to  modify  any  parameter  in  the  system  to  avoid  operating  the  freezer  in  the 
supersaturated  zone. 


CHAPTER  2 

EXPERIMENTAL  FACILITIES,  INSTRUMENTATION  AND  DATA  ACQUISITION 

Experimental  Facilities 

The  experimental  facility  utilized  in  this  investigation  is  shown  in  Figure  2-1.  A 
schematic  of  the  psychometric  freezer  room  is  shown  in  Figure  2-2.  The  figure  shows  a 
laboratory-size  freezer  with  an  industrial-size  freezer  coil  (having  four  fins  per  inch) 
located  at  the  center;  and  a water-vapor  generator  (WVG)  facing  the  coil.  Each  of  the 
freezer  doors  has  good  rubber  seals  to  protect  the  freezer  and  the  test  results  from  the 
effect  of  air  infiltration.  The  heat  transfer  rate  of  the  test  enclosure  was  determined 
experimentally  while  it  was  clean,  dry,  and  empty. 

The  finned-tube  freezer  coil  is  a liquid  overfeed  recirculating  evaporator  with  an 
overfeed  ratio  of  three.  This  coil  has  a refrigerating  capacity  of  about  2 tons  at  a coil 
suction  temperature  of  -40°C  and  is  part  of  a complete  refrigeration  system  discussed  in 
detail  by  Al-Mutawa  et  al.  (1998a,  b,  c).  A photograph  of  the  coil  front  surface  is  shown 
in  Figure  2-3,  while  coil  details  are  provided  in  Figure  2-4.  The  coil  finned  tubes  are 
arranged  in  eight  rows  in  the  direction  of  airflow  and  in  a staggered  pattern  of  38  mm  x 
33  mm,  where  the  tube  material  is  copper  having  15.9  mm  outside  diameter  and  0.46  mm 
thickness.  The  fins  are  made  of  aluminum  and  have  a flat  pattern  with  flat  edges.  The 
fins  have  a thickness  of  0.25  mm;  and  their  spacing  is  four  fins  per  inch.  The  coil  has  a 
finned  height  of  533  mm  and  a finned  length  of  737  mm;  its  outside  dimensions  are  1016 
mm  L x 610  mm  W x 627  mm  H. 
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Figure  2-1.  Experimental  facility  used  in  this  investigation 

The  hot-gas  defrosting  method  is  used  with  this  coil  when  the  coil  is  operating  in 
the  defrosting  mode.  The  finned-tube  freezer  coil  is  classified  as  a draw-through  unit 
since  the  fan  draws  the  air  against  the  refrigerant  in  a crossflow  direction  where  each 
fluid  flows  at  right  angle  to  the  other.  However,  the  tubes  are  circuiting  such  that  the  two 
fluids  will  approach  in  a counterflow-type  heat-exchanger  arrangement.  The  freezer  coil 
is  employed  with  one  fan  which  has  a diameter  of  508  mm  and  a speed  of  18  rev/s 
(1075  rpm)  and  is  operated  by  a 186-W  motor.  The  unit  has  a face  area  of  0.39  m2  and  a 
capability  to  change  the  coil  face  velocity  from  1.02  m/s  to  7.62  m/s.  The  fan  material  is 
aluminum,  while  its  hub  is  steel.  The  full-load  amperage  (FLA)  of  the  fan  motor  is 
1.7  Amps,  while  the  maximum  circuit  amperage  capacity  of  the  motor  is  15  Amps.  A 
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demand  defrost  system  was  used  to  measure  the  frost  thickness.  This  system  consists  of  a 
commercial  frost  detector  head  to  accommodate  15.9  mm  OD  tubes;  and  a defrost  control 
unit  switch.  The  detector  head  is  a U-shaped  infrared  sensor  with  optics  on  the  inside  of 
the  open  end.  An  infrared  light  beam  crosses  the  open  end  of  the  sensor  to  detect  the 
frost  thickness.  The  detector  head  is  mounted  vertically  on  a front  freezer-coil  tube  after 
cutting  a portion  of  the  coil  fins  (Figure  2-3).  The  detector  head  is  employed  with  an 
adjusting  screw  to  set  the  head  to  the  desired  frost  thickness.  The  detector  senses  the 
thickness  of  0.78  mm  when  the  adjusting  screw  is  turned  all  the  way  in,  and  senses  the 
thickness  of  9.5  mm  when  the  adjusting  screw  is  turned  all  the  way  out.  For  the  other 
frost  thickness  ranging  from  0.78  mm  to  9.5  mm,  the  detector’s  operating  manual 
provides  a table  showing  the  number  of  adjusting  screw  turns  needed  for  each  thickness. 
The  defrost  control  unit  switch  is  enclosed  in  a junction  box  and  mounted  on  the  outside 
of  the  freezer,  where  it  is  connected  to  the  detector  head  through  a 4.5  m polarized  plug 
cable.  The  unit  switch  is  a normally-opened  contact  that  closes  when  the  frost  build-up 
on  the  detector  head  reaches  the  desired  adjusted  frost  thickness,  at  which  point  the 
infrared  light  beam  is  interrupted  by  the  frost  build-up.  The  closing  of  the  unit  switch 
causes  a small  red  light  located  on  the  junction  box  to  signal  the  beam  crossing  (i.e.,  the 
desired  frost  thickness);  and  also  sends  a 5-Volt  DC  analog  signal  to  the  computer  used  in 
this  investigation. 

The  artificial  load  generation  system  is  probably  one  of  the  most  critical  systems 
in  the  experimental  program  (Figure  2-5).  It  was  designed  to  provide  the  freezer  with  the 
required  latent  heat  load  in  order  to  be  able  to  manipulate  the  moisture  content  inside  the 
freezer  during  the  testing  period. 
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Figure  2-2.  Schematic  of  the  Psychometric  Freezer  Room  showing  the  test  and  auxiliary 
coils,  the  artificial  load  generator,  and  the  location  of  thermocouples  for 
temperature  measurements 


Figure  2-3.  Photograph  showing  the  front  face  of  the  Test  Fan-Coil  Unit  (FCU-T) 
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Figure  2-4.  Schematic  of  the  Test  Fan-Coil  Unit  (FCU-T). 

This  latent  heat  load  is  generated  by  the  water- vapor  generator  (WVG)  located 
outside  the  freezer.  City  water  is  allowed  to  flow  into  the  WVG  through  an  electronic 
diaphragm  metering  pump  in  order  to  control  the  mass  flow  rate  of  the  steam  to  be 
injected  inside  the  freezer.  The  metering  pump  capacity  ranges  from  0.011  to  1.1  mm/s 
(maximum  capacity  per  day  is  0.091  m3).  The  pump’s  maximum  injection  pressure  is 
758  kPa.  This  metering  pump  has  an  adjustable  speed,  which  ranges  from  5 to  100 
strokes  per  minute,  while  it  also  has  an  adjustable  stroke  length  that  ranges  from  0 to 
100%.  The  metering  pump  can  be  operated  manually  and  by  a computer. 

The  WVG  is  a liquid-injection  type  water-vapor  generator  equipped  with  three 
heating  elements  strapped  to  its  side  and  bottom.  The  vaporizer  is  insulated  with 
25.4  mm  thick  high-temperature  insulation  enclosed  in  a galvanized  steel  housing.  The 
WVG  is  also  equipped  with  a thermometer,  a thermostat,  and  a pressure-relief  valve.  The 
bimetal  dial  thermometer  is  accurate  to  within  0.5°C  and  has  a reading  range  of  10°C  to 
288°C.  The  thermostat  adjusting  screw  was  used  to  obtain  the  desired  set  point.  The 
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pressure  relief  valve  has  a cracking  pressure  range  of  0 to  138  kPa.  The  cracking 
pressure  can  be  adjusted  to  the  desired  set  point  using  the  valve’s  adjustment  screw.  The 
water  injected  into  the  WVG  is  heated  to  the  desired  temperature  inside  the  vaporizer  and 
is  allowed  to  leave  the  WVG  as  steam  through  a copper  tube  that  passes  through  the 
freezer  wall  to  the  steam  outlet  inside  the  artificial  load  generator.  To  prevent  steam  from 
freezing  inside  the  copper  tube,  an  electric  heating  cable  is  wrapped  around  the  copper 
tube  inside  the  freezer.  The  copper  tube  and  the  heating  cable  are  both  covered  by  a 19 
mm  thick  Armaflex  insulation.  A clean-up  strip  heater  of  150-W  is  also  attached  to  the 
bottom  of  the  steam  outlet  in  order  to  prevent  the  steam  from  freezing  at  the  outlet  that 
will  then  stop  the  flow  of  steam  to  the  freezer.  The  artificial  load  generation  system  uses 
direct-drive  blower  which  distributes  the  steam  inside  the  freezer.  The  blower  has  a 
capacity  of  0.73  m3/s  at  0.025  m static  pressure  when  its  motor  is  operating  at  a speed  of 
17.5  rev/s  (1050  rpm). 

Instrumentation 

Temperatures  have  been  measured  at  twenty  locations  using  copper-constantan 
(type-T)  thermocouples  having  an  uncertainty  of  ±0.2°C.  This  type  is  suitable  for  low 
temperature  applications  as  well  as  up  to  370°C.  All  thermocouples  have  been  calibrated 
using  a constant-temperature  water  bath  (model  Polyscience-80).  Locations  of  some  of 
the  thermocouples  used  are  shown  in  Figure  2-2.  Temperature  measurement  using 
thermocouple  TC9  is  identified  as  coil  entering  air  temperature,  while  temperature 
measurement  using  thermocouple  TC10  is  identified  as  coil  leaving  air  temperature. 
Infrared  images  were  recorded  employing  a Thermacam  PM  695,  manufactured  by  FLIR 
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systems.  This  camera  can  make  accurate  non-contact  measurements  from  -40°C  to 
2000°C  with  an  uncertainty  of  ±2°C. 

The  relative  humidity  ( RH)  of  the  air  inside  the  freezer  was  measured  both 
upstream  and  downstream  of  the  test  coil.  Relative  humidity  was  measured  using  two 
Mamac  humidity  transducers  (HU-224)  connected  to  a remote  probe  hung  inside  the 
freezer.  The  two  transducers  have  a humidity  range  of  0 to  100%  and  are  accurate  to 
±2%  of  the  full-scale.  They  use  a DC  power  supply  of  12  to  28  Volts  to  operate,  sending 
output  signals  of  4-20  mA. 

The  differential  pressure  of  the  air  across  the  coil  was  measured  using  a low- 
pressure  transducer  (Mamac-PR-274)  with  two  static  pressure  probes  (Mamac-A-520). 
Both  of  the  203-mm-long  aluminum  probes  have  two  orifices  opposing  each  other  on  a 
vertical  position.  Each  probe  is  attached  to  a 50.8  mm  OD  flange  (used  to  mount  the 
probe  on  the  coil  sheet  metal  frame)  by  drilling  a 6.4  mm  hole  on  the  latter  and  inserting 
the  probe  into  that  hole.  The  probes  are  inserted  at  right  angles  to  the  airflow,  but  with 
the  orifices  parallel  to  the  airflow  so  that  static  (not  stagnation)  pressure  is  recorded. 
Each  probe  is  connected  to  the  low-pressure  transducer  by  a 6.4  mm  OD  copper  pilot  line 
having  PVC  end  connections.  The  probes  are  capable  of  operating  under  a maximum 
pressure  of  69  kPa  and  a maximum  temperature  of  121°C. 

Data  Acquisition  System 

All  measurements  were  recorded  using  a computer  equipped  with  a state-of-the-art 
data  acquisition  system.  The  data-acquisition  software  called  Lab  VIEW  for  Windows 
was  used  in  this  project.  This  software  is  a graphical  programming  language  that  has 
been  widely  adopted  throughout  industry,  academia  and  research  labs  as  the  standard  for 
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data  acquisition  and  instrumentation  control.  Figure  2-6  shows  the  screen  used  to 
monitor  the  data  during  the  experiments. 


Figure  2-5.  Schematic  of  the  artificial  sensible  and  latent  load  generator 


Temperatures  (°C),  Pressure  (kPa),  Sensible  heat  (kW) 


TCI:  FCU-T  melt  line  upstream  of  the  valve 

TC2:  FCU-T  melt  line  downstream  of  the  valve 

TC3:  FCU-T  wet  retum/hot-gas  line 

TC4:  FCU-T  pan  hot-gas  inlet  line 

TC5:  FCU-T  pan  hot-gas  outlet  line 

TC6:  FCU-T  refrigerant  liquid/defrost  refrigerant  line 

TC7:  FCU-T  wet  retum/hot-gas  line 

TC8:  FCU-T  refrigerant  liquid/defrost  refrigerant  line 

TC9:  FCU-T  coil  entrance 

TC10:  FCU-T  coil  exit 


PI:  Defrost  refrigerant  line  downstream  of  the  O/M 
P2:  FCU-T  wet  return/hot-gas  line 
P3:  FCU-T  refrigerant  liquid/defrost  refrigerant  line 
P4;  FCU-T  pan  hot-gas  inlet  line  upstream  of  the  O/M 


Sensible  heat 


TC11:  FCU-T  blower  outlet 

TCI 2:  Upstream  of  the  FCU-T  coil 

TC13:  Downstream  of  the  FCU-T  blower 

TC14:  Inside  the  freezer  next  to  the  right  wall 

TC15:  Inside  the  freezer  next  to  the  front  wall 

TC16:  FCU-A  melt  line  upstream  of  the  valve 

TCI 7:  FCU-A  melt  line  downstream  of  the  valve 

TC18:  FCU-A  refrigerant  liquid  line 

TC19:  FCU-A  wet  return  line 

TC20:  Outside  the  psychrometric  freezer  room 


P5:  Main  liquid  line  downstream  of  the  liquid  pump 
P6:  Main  suction  line  upstream  of  the  LCU  PRV 
P7:  Main  suction  line  downstream  of  the  LCU  PRV 
P8:  Main  liquid  line  downstream  of  the  main  LSV 


Figure  2-6.  Screen  for  monitoring  the  data  during  the  experiment 


CHAPTER  3 

ISSUES  PERTAINING  TO  DEMARCATION  LINE  IN  THE  FROST  ZONE 

Objectives 

One  of  the  objectives  of  our  study  was  to  experimentally  investigate  frost  formation 
on  freezer  coils  under  supersaturated  conditions;  and  to  relate  the  findings  to 
psychrometric  theory.  The  object  was  to  determine  if  there  is  a demarcation  line  between 
the  type  of  frost  that  forms  under  subsaturated  conditions  and  the  frost  that  forms  under 
supersaturated  conditions.  Other  objectives  included  providing  a critical  evaluation  of 
the  current  state  of  knowledge  of  the  frosting  problem  in  general;  and  discussing  the 
differences  and  similarities  as  they  pertain  to  freezer  entryways  in  particular.  Frost 
formation  in  freezer  entryways  usually  occurs  under  supersaturated  conditions  because  of 
peculiar  infiltration  and  air  cooling  processes  (discussed  later).  Thus  we  analyzed  the 
frost  formation  under  both  subsaturated  and  supersaturated  conditions.  In  cases  where 
the  current  state  of  knowledge  is  sketchy,  an  attempt  was  made  to  hypothesize  the 
phenomenon  in  question.  The  hypotheses  presented  in  this  dissertation  were  partially 
based  on  preliminary  results  of  this  experimental  investigation. 

Establishing  the  Existence  of  the  Demarcation  Line 
As  mentioned  above,  one  of  the  objectives  of  this  experimental  program  was  to 
demonstrate  the  existence  of  a demarcation  line  between  the  frost  that  forms  under 
subsaturated  conditions  to  that  that  forms  under  supersaturated  conditions.  This  was 
achieved  by  starting  with  a predetermined  coil  entering  air  temperature,  a predetermined 
room  (freezer)  sensible  heat  ratio  ( SHR ),  and  a predetermined  coil  surface  temperature, 
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and  then  steadily  decreasing  the  coil  sensible  heat  ratio  until  the  cooling  process  line  has 
crossed  the  saturation  curve  on  the  psychrometric  chart  into  the  supersaturated  zone. 
Decreasing  the  coil  sensible  heat  ratio  was  achieved  by  injecting  known  amounts  of 
steam  into  the  freezer  while  keeping  the  coil  entering  air  temperature,  and  the  coil  surface 
temperature  constant.  In  this  experiment,  steam  injection  was  conducted  at  a rate  of  1.5 
kg/h.  This  resulted  in  increasing  the  relative  humidity  of  the  air  inside  the  freezer 
(decreasing  the  coil  sensible  heat  ratio),  while  keeping  constant  the  coil  entering  air 
temperature  at  -8.3°C.  Freezer  temperature  and  sensible  heat  ratio  were  controlled  by 
either  adding  an  artificial  sensible  load  to  the  freezer  or  by  increasing  the  cooling  rate 
inside  the  freezer  employing  the  refrigeration  system  (without  changing  the  coil  surface 
temperature).  The  starting  relative  humidity  of  the  coil  entering  air  was  64%,  while  the 
coil  surface  temperature  was  kept  constant  at  -18.3°C.  As  steam  injection  continued,  the 
relative  humidity  of  the  coil  entering  air  was  increased  until  it  reached  99%.  Four 
different  scenarios  were  tried  employing  the  same  entering  air  temperature,  and  coil 
surface  temperature  (see  Figure  3-1). 

In  Scenario  1,  the  inlet  dry-bulb  temperature  ( DBT)  was  kept  at  -8.3°C  and  the 
inlet  relative  humidity  ( RH)  was  64%,  while  the  leaving  air  DBT  was  -16.2°C  and  the 
leaving  air  RH  was  92%.  The  degree  of  saturation  (jud)  for  this  scenario  was  91.9%.  In 
Scenario  2,  the  inlet  DBT  was  -8.3°C  and  the  RH  was  84%,  while  the  leaving  DBT  was 
-16.2°C  and  the  degree  of  saturation  (/^>sup(?r)was  102%.  In  Scenario  3,  the  inlet  DBT 

was  -8.3°C  and  RH  92%,  while  the  leaving  DBT  was  -16.1°C  and  judsuper  was  105%.  In 


Scenario  4,  the  inlet  DBT  was  -8.3°C  and  RH  was  99%,  while  the  outlet  DBT  was 
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-16. 1°C  and  fid  super  was  109%.  It  is  important  to  observe  that  the  relative  humidity  is  not 

defined  in  the  supersaturated  zone  of  the  psychrometric  chart.  An  alternative  term  that  is 
equally  representative  of  the  amount  of  moisture  present  in  the  air  in  that  zone  is  the 
degree  of  saturation  y/djSuper).  This  is  the  ratio  between  the  supersaturated  humidity  ratio 

(O  and  the  saturated  humidity  ratio  (Ws)  at  the  prevailing  dry-bulb  temperature. 


Figure  3-1.  Psychrometric  depiction  of  air-cooling  processes  as  measured  inside  the 
freezer 

While  psychrometric  depictions  of  all  four  processes  are  shown  in  Figure  3-1, 
actual  freezer  conditions  corresponding  to  the  four  scenarios  in  question  are  shown  in 
Figures  3-2  through  3-5,  respectively.  For  Scenario  1,  where  the  cooling  process  line 
(Line  1-ADP)  is  completely  located  in  the  subsaturated  zone,  clear  conditions  (no  fog) 
can  be  observed  inside  the  freezer.  For  Scenario  2,  where  a small  part  of  the  cooling 
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process  line  (Line  la-ADP)  is  located  in  the  supersaturated  zone  on  the  psychrometric 
chart,  a slight  amount  of  fog  can  be  detected  (Figure  3-3).  For  Scenario  3,  where  a bigger 
portion  of  the  cooling  process  line  (Line  lb-ADP)  is  located  in  the  supersaturated  zone,  a 
significant  amount  of  fog  can  be  detected  in  the  freezer  (Figure  3-4).  Finally,  for 
Scenario  4,  where  the  cooling  process  line  (Line  lc-ADP)  is  fully  located  in  the 
supersaturated  zone,  fog  inside  the  freezer  is  so  heavy  that  visibility  is  impaired  (Figure 
3-5). 

A close  examination  of  the  nature  of  frost  forming  under  the  aforementioned 
scenarios  reveals  an  interesting  phenomenon.  Frost  forming  under  subsaturated 
conditions  was  found  to  be  more  fluffy  and  dryer  than  that  which  deposited  under 
supersaturated  conditions.  Frost  forming  under  subsaturated  conditions  originated  as  a 
result  of  a transformation  of  the  water  vapor  in  air  from  a vapor  phase  to  a solid  phase 
bypassing  the  liquid  phase.  It  is  a diffusion-dominated  phenomenon.  Frost  forming 
under  severe  supersaturated  conditions  continued  until  the  coil  was  completely  clogged 
(Figure  3-6).  Frost  formation  under  those  conditions  is  thought  to  be  a combination  of 
two  processes:  a deposition  process  (dominated  by  convection)  in  which  air-borne  ice 
crystals  deposit  on  the  freezer  coil  surface,  and  a diffusion-dominated  process  similar  to 
the  one  observed  under  subsaturated  conditions.  The  wetter  and  denser  nature  of  this 
type  of  frost  is  partly  due  to  the  fact  that  frozen  water  droplets  (or  ice  crystals)  exhibit 
wetter  and  denser  characteristics  than  the  more  traditional  frost.  Frost  forming  under 
these  conditions  is  thought  to  be  more  snow-like  and  is  not  as  favorable  from  a heat 
transfer  performance  standpoint. 
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Figure  3-2.  Visualization  of  freezer  conditions  during  Scenario  1 as  measured  inside  the 
freezer 


Scenario  2 
Inlet:  DBT=- 8.3°C 

RH= 84% 

Outlet:  DSr=-16.2°C 

Md, super  = 102% 

Mid-Coil:  DB7=-1 2.2  °C 

Md, super  ~ 105% 
Face  Velocity:  3.8  m/s 


Figure  3-3.  Visualization  of  freezer  conditions  during  Scenario  2 as  measured  inside  the 
freezer 
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Scenario  3 


Inlet: 
Outlet: 


DBT=- 8.3°C 
RH=  92% 

DB  T=-1 6. 1°C 

A#, super  “ 105% 

Mid-Coil:  DBT=A2.2°C 

Md.super~  115% 

Face  Velocity:  3.8  m/s 


SSI 


Figure  3-4.  Visualization  of  freezer  conditions  during  Scenario  3 as  measured  inside  the 
freezer 


Figure  3-5.  Visualization  of  freezer  conditions  during  Scenario  4 as  measured  inside  the 
freezer 
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Figure  3-6.  Frosted  coil  after  severe  supersaturated  conditions  have  existed  in  the  freezer 
for  a while 

In  order  to  demonstrate  the  significance  of  the  existence  of  supersaturated 
conditions  during  the  deposition  process,  both  steam  injection  and  the  refrigeration 
process  were  abruptly  stopped  after  the  coil  has  been  completely  blocked.  Freezer  air 
was  shortly  restored  to  subsaturated  conditions  and,  after  a relatively  short  time  period, 
the  coil  started  to  become  unclogged  again  (Figure  3-7).  What  was  actually  happening 
was  that  the  frost  on  the  coil  was  beginning  a sublimation  process  (transformation  from 
solid  to  vapor  bypassing  the  liquid  phase). 

This  was  primarily  a reflection  of  the  fact  that  the  vapor  pressure  of  the  air  in  the 
vicinity  of  the  frost  surface  was  getting  higher  than  that  in  the  bulk  freezer  air. 
Sublimation  is  a diffusion-dominated  phenomenon  not  unlike  the  transformation  of  water 
vapor  to  frost.  Keeping  freezer  air  subsaturated  is  essential  if  frequent  defrosting  of  coils 
is  to  be  avoided.  Based  on  the  above,  a frost  formation  theory  will  be  articulated  for  the 
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formation  mechanisms  in  both  the  subsaturated  and  the  supersaturated  zones.  This  will 
be  introduced  after  a model  that  describes  the  cooling  process  through  a frosted  coil  has 
been  developed. 


Figure  3-7.  Frosted  coil  due  to  sublimation  of  the  frost  after  subsaturated  conditions  have 
been  restored  in  the  freezer  after  severe  supersaturation  of  freezer  air 

Modeling  the  Cooling  Process  Path  through  a Frosted  Coil 
In  our  continuing  effort  to  provide  more  insight  into  the  frost-formation 
mechanism  and  how  it  relates  to  the  prevailing  psychrometric  conditions,  we  will  attempt 
to  calculate  the  path  on  an  actual  industrial-size  finned-tube,  multi-row  freezer  coil.  This 
will  be  done  employing  basic  heat  and  mass  transfer  principles.  The  model  will  then  be 
validated  employing  measured  data  on  the  same  coil.  Development  of  a systematic 
mechanism  for  calculating  the  cooling  process  path  on  a finned  freezer  coil  will  prove 
very  significant  towards  enhancing  our  understanding  of  the  frost-formation  process  and 
its  relationship  to  the  psychrometric  theory  of  moist  air.  Determining  the  actual 
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conditions  leaving  a given  row  in  a multi-row  freezer  coil  is  a crucial  step  in  identifying 
the  coil  location  in  the  vicinity  of  which  the  transformation  from  subsaturated  conditions 
to  supersaturated  conditions  occurs.  This  is  a key  step  in  identifying  a demarcation  line 
between  the  unfavorable  snow-like  frost  and  the  more  traditional  (and  more  favorable) 
frost-formation  patterns.  Thus,  the  objective  of  this  effort  is  to  calculate  the  air  path  on 
an  actual  industrial-size  finned-tube,  multi-row  coil  utilizing  experimentally  derived  data 
and  correlate  the  shape  of  the  path  with  the  prevailing  psychrometric  conditions  in  the 
freezer  in  the  hope  of  identifying  the  demarcation  line  in  question. 

Analysis 

The  primary  purpose  of  the  analysis  described  in  this  section  is  calculating  the 
shape  of  the  cooling  process  line  as  the  air  passes  through  a finned-tube,  multi-row 
dehumidifying  coil  under  frosting  conditions.  In  order  to  achieve  this  goal,  the  slope  of 
the  enthalpy-humidity  ratio  curve  on  the  psychrometric  chart  must  be  determined.  In 
order  to  calculate  this  slope,  several  other  quantities  such  as  the  overall  heat  transfer 
coefficient  and  the  frost-air  interface  temperature  have  to  first  be  computed.  These 
quantities  are  determined  according  to  what  follows. 

External  heat  transfer  coefficient.  To  determine  the  external  heat  transfer 
coefficient  a staggered  array  is  considered  as  shown  in  Figure  3-8.  For  this  configuration, 
the  maximum  velocity  may  occur  at  either  the  transverse  plane  ( A, ) or  the  diagonal  plane 
(A2)  of  Figure  3-8.  Here  we  use  the  diameter  of  the  frosted  tube  [pf]  instead  of  the 
outside  tube  diameter.  The  maximum  velocity  occurs  at  the  diagonal  plane  A2  if 
SD  <0.5 (Sj-  + Z)y). 
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In  this  case,  the  maximum  velocity  is  given  by 

ST 


2(S„  -D,?' 


(3-1) 


max 


If  the  maximum  velocity  occurs  at  the  transverse  plane  A, , it  can  be  expressed  as 


Figure  3-8.  Staggered  array  of  a coil 

The  Reynolds  number  based  on  the  fluid  velocity  occurring  within  the  tube  bank 
is  given  by 


The  Nusselt  number  for  a tube  bank  composed  of  10  or  more  rows  (NL  > 10 ) can 
be  determined  using  the  correlation  proposed  by  Grimison  (1937) 


(3-2) 


u, 


max 


(3-4) 


^>10 

2000  < Re0  < 40000 
Pra  > 0.7 
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If  the  number  of  rows  ( NL ) is  less  than  1 0 a correction  factor  may  be  applied  such  that 


Nu 


= C,  Nu 


(3-5) 


where  C,  and  j are  listed  in  Table  3-1  while  C2  is  listed  in  Table  3-2. 


Table  3-1.  Constants  Ct  and  j for  a staggered  arrangement  (Equation  3-4)  for  air  flow 


over  a tube  bank  of  10  or  more  rows. 


SJD, 

11 

cT 

1.25 

Sr/Df 

1.50 

11 

cT 

2.00 

11 

cT 

to" 

3.00 

C\ 

j 

C\ 

j 

c, 

j 

C, 

j 

0.60 

- 

- 

- 

- 

- 

- 

0.213 

0.636 

0.90 

- 

- 

- 

- 

0.446 

0.571 

0.401 

0.581 

1.00 

- 

- 

0.497 

0.558 

- 

- 

- 

- 

1.25 

0.518 

0.556 

0.505 

0.554 

0.519 

0.556 

0.552 

0.562 

1.50 

0.451 

0.568 

0.460 

0.562 

0.452 

0.568 

0.488 

0.568 

2.00 

0.404 

0.572 

0416 

0.568 

0.482 

0.556 

0.449 

0.570 

3.00 

0.310 

0.592 

0356 

0.580 

0.440 

0.562 

0.428 

0.574 

Table  3-2.  Correction  factor  C2  for  a staggered  arrangement  as  provided  for  in  Equation 


Number  of  rows  (NL) 

Correction  factor  ( C2 ) 

1 

0.68 

2 

0.75 

3 

0.83 

4 

0.89 

5 

0.92 

6 

0.95 

7 

0.97 

8 

0.98 

9 

0.99 

The  external  heat  transfer  coefficient  is  then  expressed  by 


h 


c,o 


(3-6) 


Internal  heat  transfer  coefficient.  The  properties  are  evaluated  at 


the  mean 


refrigerant  temperature  between  the  tube  inlet  and  outlet  TRm 


)/2 


The  Reynolds  number  for  a flow  in  a circular  tube  is  given  by 
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Re,  = 


_ pRumA 


Mr 


(3-7) 


where  um , is  the  mean  velocity  inside  the  tube  and  Z>  is  the  inside  tube  diameter.  The 
mean  velocity  is  defined  by 


u 


m,i 


PrAc 


(3-8) 


where  Ac  is  the  cross-sectional  area  of  the  tube  and  mR  is  the  mass  flow  rate.  The 
Reynolds  number  can  be  reduced  to 


r.  4 mR 
Re , = R 


7lD^R 


(3-9) 


For  turbulent  flow,  the  Nusselt  number  can  be  calculated  using  the  Dittus-Boelter 
correlation  (1930) 

Nut  = 0.023 Re,°Vr°4  (3-10) 


The  internal  heat  transfer  coefficient  is  then  given  by 


h . = Nu(  — 
c *'  ' D 


(3-11) 


Frost  surface  temperature.  Referring  to  Figure  3-9,  the  heat  transfer  rate  can  be 
expressed  by 

Q = KM>r*fc'-Tf)  (3'12) 

as  well  as  by 


2nkfl(Tf-Tp)  (3-13) 

L°g\D,lD.) 


which  if  combined  can  give  the  following  expression  for  the  frost  surface  temperature: 
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T kriJ>fLog'(p{/D')r+tt,Tr  O'14) 

’ *„fl/io&(0//D.)+2*/ 

where  kf  is  the  thermal  conductivity  of  the  frost  in  W/(m  K). 


a 

\ 


V/ 

: Ir 

Figure  3-9.  Physical  model  of  pipe 

Several  empirical  and  theoretical  correlations,  relating  the  frost  thermal 
conductivity  and  density  to  the  frost  surface  temperature,  exist.  In  this  analysis,  the 
correlation  by  Yonko  and  Sepsy  (1967)  will  be  used  to  calculate  the  frost  thermal 
conductivity,  while  that  by  Hayashi  et  al.  (1977a,b)  will  be  used  to  compute  the  frost 
density.  This  gives  the  following  equations,  respectively: 

kf  = 0.024248  + 0.00072311/^+0.000001183^/  (3-15) 

pf  =650exp(o.227t/)  (3-16) 

Equation  3-16  is  valid  for  a frost  surface  temperature  in  the  range  -25°C  < tf< 

0°C  and  for  an  air  stream  velocity  between  2 and  6 m/s.  In  this  analysis,  the  temperature 
of  the  frost  is  calculated  twice  for  the  purpose  of  calculating  the  overall  heat  transfer 
coefficient,  one  time  at  the  entrance  and  another  time  at  the  exit  of  the  coil. 

Overall  heat  transfer  coefficient  for  a frosted  finned- tube  coil.  In  this  section  we 
will  determine  the  overall  heat  transfer  coefficient  for  the  case  when  the  fins  and  tubes 


are  frosted.  To  achieve  this  we  assume  that  the  thermal  resistance  of  the  tube  is 
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negligible,  that  the  tube  mean-surface  temperature  is  tp , and  that  the  refrigerant  mean 
temperature  is  tR . It  will  also  be  assumed  that  a thin  layer  of  frost  having  an  average 
thickness  equal  to  X f covers  the  fins  and  tubes.  The  local  rate  of  heat  transfer  inside  the 
tube  can  be  written  as 

Q = KMb~h)  <3'17> 


The  enthalpy  of  saturated  moist  air  can,  in  general,  be  expressed  as  is  = a + bts  for 
a small  temperature  range,  such  as  5°C.  The  term  b is  the  slope  of  the  saturated  enthalpy- 
temperature  curve  (i.e.,  b = Ais/Ats  )■  Thus,  a quantity  b'R  can  be  defined  according  to 
the  equation 


b\  = 


ls,P  ls,R 
lp  ~lR 


(3-18) 


In  Equation  3-18,  the  subscripts  P and  R indicate  conditions  pertaining  to  the  pipe 
surface  and  refrigerant  mean  temperatures,  respectively. 

From  Equations  3-17  and  3-18  we  get 


Q = M/v 


{h,P  h,R ) 


b\ 


(3-19) 


The  rate  of  heat  transfer  from  the  tube  and  fins  to  the  air  can  be  expressed  by 


q = Ap  o ( i - is  p )+  Af  (i  _ iF  m ) 

bf,p  bf  m 


(3-20) 


where  bf  p and  bf  m are  the  slopes  of  the  saturated  enthalpy-temperature  curve  evaluated 


at  the  frost  temperature  on  the  tube  surface  and  on  the  fin,  respectively.  Both  of  these 
quantities  can  be  evaluated  employing  moist  air  tables.  The  quantity  iF  m is  the  enthalpy 
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of  saturated  air  at  the  mean  fin  temperature.  The  term  ha  f is  the  external  heat  transfer 
coefficient  for  a frosted  surface  and  is  given  by 

, 1 (3-21) 

Koj  ~ 


( c ) 

P,“ 

M 

\bfho) 

kr) 

where  kf  is  calculated  by  Equation  3-15. 

The  overall  heat  transfer  coefficient  can  be  expressed  in  terms  of  the  fin  efficiency.  The 
efficiency  of  a frosted  fin  is  defined  as 

(3-22) 


i -i 


Vf 


F,m 


l -l 


F,B 


where  iFm  and  iFB  are  the  enthalpies  of  saturated  air  evaluated  at  the  mean  fin 

temperature  and  fin  base  temperature,  respectively. 

Assuming  that  the  temperature  of  the  frost  surface  for  the  part  of  the  frost  on  the 
tube  to  be  approximately  the  same  as  the  surface  temperature  of  the  frost  for  the  part  of 
the  frost  on  the  fin  base,  the  term  bfp  can  be  made  equal  to  the  term  bf  m . Furthermore, 
the  temperature  of  the  pipe  surface  can  be  assumed  equal  to  the  temperature  of  the  fin 
base,  thus  resulting  in  the  following  two  enthalpy  terms  being  equal:  is  P = iFB . This 
simplifies  Equation  3-20  to  the  following  form  after  incorporating  Equation  3-22: 


Q = ~7~~  Wo  + Vf^f  X*  - h,p ) 


(3-23) 


f,m 


By  definition,  the  overall  heat  transfer  coefficient  based  on  the  enthalpy 
difference,  Ua  f , may  be  expressed  as 

Q = U0jA0{i-W 


(3-24) 
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Employing  Equations  3-19,  3-23,  and  3-24  we  get 


(3-25) 


The  overall  heat  transfer  coefficient  for  a frosted  finned  tube  coil  based  on  the 
enthalpy  difference  is  calculated  using  the  above  procedure  at  two  locations,  one  at  the 
entrance  and  the  other  at  the  exit  of  the  coil.  An  average  value  is  then  computed. 

Calculated  path  of  the  cooling  and  dehumidifving  process.  As  has  been  explained 
earlier,  when  simultaneous  cooling  and  dehumidification  occurs,  the  overall  heat  transfer 
coefficient,  UoJ,  is  defined  in  terms  of  an  enthalpy  (and  not  a temperature)  driving 

force.  Furthermore,  the  enthalpy  “z”  is  the  true  air  enthalpy.  This  is  to  be  contrasted  with 
the  quantity  isR  which  can  be  thought  of  as  a fictitious  saturated  air  enthalpy  evaluated  at 

the  refrigerant  temperature.  The  mean  enthalpy  difference  A im  is  defined  by  the  equation 


The  mean  enthalpy  difference  Azmis  defined  by  [modified  from  Incropera  and  Dewitt 


were  z,  and  z'2  are  the  true  enthalpy  values  of  the  entering  and  leaving  air  conditions, 
respectively,  and  is  R , and  is  R 2 are  the  fictitious  enthalpies  of  saturated  air  calculated  at 


(3-26) 


(1996)] 


(3-27) 


the  refrigerant  entering  and  leaving  conditions,  respectively. 
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After  calculating  the  inlet  air  conditions,  the  inlet  refrigerant  state,  the  various 
heat  transfer  coefficients,  and  the  cooling  coil  surface  data,  the  path  of  the  cooling 
process  can  be  determined.  To  get  the  enthalpy  of  saturated  air  at  the  mean  frost  surface 
temperature,  we  use  the  following  equation: 


* “ h,m  = VF  (*  - h,P  ) = Tf~pr-  (*  ~ is,f,nt  ) 
'lc,oj'~'p,a 


(3-28) 


And  from  Equations  3-19  and  3-24, 


i-i 


■s,P 


\ 


bc,Ap,i  ) 


(3-29) 


Using  Equations  3-28  and  3-29,  the  enthalpy  of  saturated  air  evaluated  at  the 
mean  frost  surface  temperature  can  be  expressed  by 


it  —i  — 

s,f,m 


C h ftlrf 

• p,a  c,o,f  IF 


bf,A,o 


{ b\UoJA0 
KAp>  J 


(3-30) 


For  a given  coil  entering  air  condition,  the  quantity  can  be  calculated.  This 
determines  State  “s,fm,  1”  on  the  psychrometric  chart.  This  allows  us  to  compute  the 
quantity  i,  -z'j/m  JWy  -WsfmX.  The  slope  at  the  coil  entering  air  state,  ( di/dW\ , can 


then  be  calculated  according  to  the  equation  (Kuehn  et  al.  1998) 

f di  ^ 


= Le  h ^L+| 

KdWjx 


fc.i  ~i0gLe) 


(3-31) 


where  Le  is  the  Lewis  number,  i , is  the  enthalpy  of  saturated  water  vapor  at  the  air  dry- 


bulb  temperature,  and  i°  is  the  enthalpy  of  saturated  water  vapor  at  a reference 
temperature  (0°C).  With  the  aid  of  the  chart,  we  draw  a short  segment  of  direction 
(di/dW\  through  State  1.  With  an  arbitrary  value  of  enthalpy  we  can  locate  a new  point 
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on  this  segment.  Then  the  procedure  is  repeated  until  the  condition  line  meets  with  the 
state  of  saturated  air  at  the  coil  apparatus  dew-point  temperature  (approximately  equal  to 
the  coil  refrigerant  temperature).  Having  determined  the  path  of  the  cooling  process,  we 
can  predict  the  condition  downstream  of  any  row  using  the  following  equation  (Kuehn  et 
al.  1998): 

. (l-e-(l-Cl^+dl-C,>-(,-C,)C'  0-32) 

where  C3  =mabR/mRCpR  and  C4  =UoJA0/ma . In  computing  the  quantity C4,  it 
should  be  kept  in  mind  that  the  numerical  value  of  Aa  encompasses  the  areas  of  the  row 

being  calculated  as  well  as  all  rows  upstream  of  that  row. 

Results  and  Discussion 

As  was  reported  earlier,  the  objective  of  this  section  was  to  calculate  the  path  the 
air  would  take  as  it  passes  through  an  industrial-size,  finned-tube,  multi-row  freezer  coil. 
The  calculated  path  is  to  be  then  compared  to  the  path  obtained  from  the  straight-line 
principle  of  applied  psychrometrics  in  order  to  assess  the  goodness  of  the  theoretical 
model.  Table  3-3  provides  a summary  of  the  geometrical  data  of  the  coil  under  study, 
whereas  the  values  of  the  parameters  employed  in  the  model  are  given  in  Table  3-4. 
Table  3-5  gives  a summary  of  the  numerical  results  obtained  by  applying  the  equations 
used  to  compute  the  overall  heat  transfer  coefficient. 
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Table  3-3.  Geometrical  data  of  the  finned  coil 


Tube  outside  diameter,  mm 

15.88 

Tube  inside  diameter,  mm 

14.97 

Tube  spacing  between  rows,  mm 

38.10 

Tube  spacing  across  coil-face,  mm 

32.80 

Thickness  of  fins,  mm 

0.25 

Net  outside  surface  area,  m2 

6.23 

Inside  surface  area,  m2 

6.09 

Fin  surface  area,  m2 

43.37 

Total  outside  surface  area,  m2 

49.60 

Table  3-4.  Parameter  values  employed  in  the  model 

Air  inlet  temperature,  °C 
Air  inlet  relative  humidity,  % 

Air  face  velocity,  m/s 

Pipe  surface  temperature  at  the  inlet,  °C 

Pipe  surface  temperature  at  the  outlet,  °C 

Pipe  mean  surface  temperature,  °C 

Frost  thickness,  mm 

Frost  temperature  at  the  coil  inlet,  °C 

Frost  temperature  at  the  coil  outlet,  °C 

Frost  mean  temperature,  °C 

Fins  efficiency,  % 


0 

65 

3.8 

-10 

-2.2 

-6.1 

1 

-6.1 

-1.8 

-3.9 

85 


Table  3-5.  Heat  transfer  coefficients  using  Equations  3-6,  3-11,  3-21,  and  3-27 




Internal  heat  transfer  coefficient,  W/(m  K)  13471 

External  heat  transfer  coefficient,  W/(m2  K)  107.8 

Heat  transfer  coefficient  for  a frosted  surface,  W/(m  K)  71.5 

Overall  heat  transfer  coefficient  for  a frosted  finned-tube  coil  based  on  an  0.075 

enthalpy  driving  force,  (kW  kg)/(m  . kJ) 


In  order  to  determine  the  calculated  path  for  the  conditions  in  Table  3-4,  the 
surface  area  of  the  coil  analyzed  is  arbitrarily  divided  into  nine  segments,  thus  producing 
ten  state  points  through  the  coil.  A summary  of  the  conditions  of  these  ten  points  is 
provided  in  Table  3-6.  Figure  3-10  compares  both  the  calculated  and  the  straight-line 
dehumidifying  paths.  It  is  obvious  that  a path  based  on  a straight-line  process  is  in  error 
to  varying  degrees,  depending  on  the  location  of  the  point  of  interest  relative  to  the  coil. 
This  can  be  easily  demonstrated  by  comparing  the  relative  humidity  values  of  both  the 
straight-line  and  calculated  paths  at  a specific  dry-bulb  temperature. 


37 


Figure  3-11  displays  both  the  calculated  path  showing  the  outlet  conditions  from 
successive  rows  in  the  coil  as  well  as  the  entering  and  leaving  air  conditions  based  on 
experimental  measurements.  For  example,  experimental  results  for  the  leaving  air 
conditions  are:  DBT=  -8.2°C  and  RH=91%,  while  calculated  leaving  conditions  are: 
DBT=  -8.4°C  and  RH=91 .4%,  as  is  evident  from  examining  Figure  3-11.  Thus,  there  is  a 


slight  difference  of  0.2°C  in  the  dry-bulb  temperature  and  0.4%  in  the  relative  humidity. 


The  differences  in  both  temperature  and  humidity  are  insignificant  and  lie  within  the 
uncertainty  bounds  of  both  variables.  As  can  be  observed,  the  outlet  conditions  become 
more  closely  spaced  towards  the  last  rows  of  the  coil,  implying  that  more  cooling  is 


performed  in  the  first  half  of  the  coil.  This  is  because  the  enthalpy  driving  force  assumes 
smaller  values  as  the  air  passes  through  the  coil  rows. 


Table  3-6.  Calculated  cooling  process  path  through  the  coil 


Point 

*(°  c) 

RH  (%) 

W (kgw/kga) 

i (kJ/kg) 

h,/.m  (kJ/kg) 

(di/dW) 

1 

0 

65.0 

0.00240 

6.04 

-1.35 

21310.36 

2 

1 

70.0 

0.00235 

4.92 

-1.79 

19340.68 

3 

2 

74.2 

0.00229 

3.76 

-2.24 

16783.44 

4 

3 

78.5 

0.00223 

2.60 

-2.69 

14823.96 

5 

4 

82.6 

0.00216 

1.40 

-3.16 

14094.55 

6 

5 

86.2 

0.00206 

0.18 

-3.63 

13675.00 

7 

6 

89.8 

0.00198 

-1.06 

-4.11 

11942.40 

8 

7 

92.8 

0.00188 

-2.32 

-4.61 

11427.26 

9 

8 

95.8 

0.00178 

-3.58 

-5.10 

9790.76 

10 

9 

97.5 

0.00163 

-4.88 

-5.60 

- 

Now  that  the  cooling  process  path  has  been  found,  the  coil  location  in  the  vicinity 


of  which  the  transformation  from  a subsaturated  state  to  a supersaturated  state  occurs,  can 
be  determined.  In  order  to  achieve  this  goal,  we  will  assume  that  the  percentage  of 
energy  removed  by  a given  row  of  the  coil  relative  to  the  total  amount  of  energy  that  the 
coil  can  remove  remains  the  same  whether  the  air  exists  in  a subsaturated  state  or  a 


supersaturated  state.  Employing  Figure  3-1 1,  it  is  relatively  easy  to  observe  that  the  total 
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amount  of  cooling  per  unit  mass  performed  by  the  coil  is  10.1  kJ/kg.  In  Table  3-7,  the 
amount  of  cooling  per  unit  mass  performed  by  each  row  and  the  percentage  that  that 
amount  represents  relative  to  the  total  amount  of  cooling  are  provided.  As  reported 
earlier,  more  cooling  is  performed  in  the  first  row  relative  to  subsequent  rows  as  shown  in 
Figure  3-12.  Carrying  these  cooling  percentages  to  the  supersaturated  zone,  the  coil 
location  in  the  vicinity  of  which  transformation  from  a subsaturated  condition  to  a 
supersaturated  condition  occurs  can  be  calculated. 


Figure  3-10.  Representation  of  the  actual  and  theoretical  air  paths  through  the 
dehumidifying  coil  (points  on  the  chart  are  arbitrarily  chosen  for  purposes 
of  calculating  the  path) 
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Table  3-7.  Amount  of  cooling  and  percentage  of  the  total  cooling  for  each  row 


Row 

A i (kJ/kg) 

Percentage  of  total  cooling  (%) 

R1 

2.64 

26.14 

R2 

2.00 

19.80 

R3 

1.61 

15.94 

R4 

1.23 

12.18 

R5 

1.02 

10.10 

R6 

0.74 

7.33 

R7 

0.49 

4.85 

R8 

0.37 

3.66 

Figure  3-11.  Representation  of  the  actual  air-cooling  path  through  the  dehumidifying  coil 
showing  the  calculated  leaving  air  conditions  from  each  row  as  well  as  the 
measured  coil  inlet  and  outlet  states 


An  experiment  has  been  performed  to  verify  the  validity  of  the  aforementioned 
assumption.  In  the  experiment,  the  entering  air  temperature  was  kept  constant  at  -8.3°C, 
while  the  relative  humidity  was  steadily  increased  from  64%  to  99%.  The  apparatus 
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dew-point  was  kept  constant  at  -18.3°C  throughout  the  experiment.  The  apparatus  dew- 
point is  the  lowest  temperature  that  the  air  can  achieve  during  the  cooling  and 
dehumidifying  process  and  is,  thus,  approximately  equal  to  the  coil  refrigerant 
temperature.  Results  of  this  experiment  are  shown  on  the  psychrometric  chart  of  Figure 
3-13  for  both  the  straight-line  theory  and  the  calculated  path.  As  can  be  seen,  four 
scenarios  were  examined.  For  Scenario  1,  the  cooling  process  line  falls  completely  in  the 
subsaturated  zone.  For  Scenario  2,  transformation  to  the  supersaturated  zone  using  the 
straight-line  theory  occurs  at  Point  3 a (DBT=-  15.1°C),  while  it  occurs  at  Point  3a’ 
(DBT=-  12.3°C)  using  the  calculated  path.  For  this  scenario,  the  amount  of  cooling 
achieved  using  the  straight-line  theory  represents  84%  of  the  total  cooling  performed  by 
the  coil.  This  means  that  at  Point  3a  the  air  was  leaving  the  fifth  row  supersaturated. 
However,  the  amount  of  cooling  achieved  using  the  calculated  path  represents  only  50% 
of  the  total  cooling  performed  by  the  coil.  This  means  that  at  Point  3a’  the  air  was 
leaving  the  third  row  supersaturated  (as  opposed  to  the  fifth  row  for  the  straight-line 
theory).  For  Scenario  3,  supersaturated  conditions  using  the  straight-line  theory  began  at 
Point  3b  (DBT=- 1 1.3°C)  and  using  the  calculated  path  at  Point  3b  ’ (DBT=-  10.3°C).  For 
this  scenario,  the  amount  of  cooling  achieved  using  the  straight-line  theory  represents 
37%  of  the  total  cooling  performed.  This  means  that  transformation  to  the  supersaturated 
zone  occurred  between  the  second  and  third  rows  (Point  3b).  However,  the  amount  of 
cooling  achieved  using  the  calculated  path  represents  only  24%  of  the  total  cooling 
performed.  This  means  that  at  Point  3b’  the  air  was  leaving  the  first  row  supersaturated. 
For  Scenario  4,  supersaturated  conditions  started  at  Point  3c  (DBT=- 8.7°C)  for  both  the 
straight-line  theory  and  the  calculated  path.  For  this  scenario,  the  amount  of  cooling 
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achieved  from  Points  lc  to  3c  or  from  Points  lc  to  3c  ’ represents  5%  of  the  total  cooling 
performed.  This  means  that  transition  to  the  supersaturated  state  occurred  in  the  vicinity 
of  the  first  row.  Obviously,  as  the  relative  humidity  of  the  entering  air  increases, 
transition  to  a supersaturated  state  occurs  quicker  and  at  a more  upstream  location  of  the 
coil.  This  phenomenon  has  also  been  confirmed  by  other  investigations  performed  by  the 
authors  (see  Sherif  et  al.  2001). 
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Figure  3-12.  Percentage  of  the  total  cooling  for  each  coil  row 

As  was  pointed  out  before,  the  procedure  employed  in  this  section  is  based  on  the 
premise  that  the  percentage  of  cooling  achieved  between  two  consecutive  rows  of  the  coil 
will  remain  the  same  whether  the  air  exists  in  the  subsaturated  or  the  supersaturated  zone. 
In  order  to  determine  the  percentage  of  cooling  referred  to  above,  the  state  of  the  air  in 
the  vicinity  of  two  consecutive  rows  has  to  be  known.  The  procedure  described  for 
determining  the  enthalpy  at  the  exit  of  a given  row  is  thus  necessary  for  computing  the 
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percentage  of  cooling.  Employing  this  procedure  necessarily  produces  a curved  path  for 
the  air  passing  through  the  coil.  While  the  bulk  of  the  air  paths  lie  in  the  subsaturated 
zone  for  all  four  processes  shown,  the  procedure  outlined  is  inherently  capable  of 
handling  both  supersaturated  and  subsaturated  conditions.  Using  the  straight-line  theory 
to  determine  the  air  enthalpy  in  the  vicinity  of  a given  row  produces  results  that  are 


significantly  less  accurate  than  the  ones  produced  using  the  curved-path.  Figure  3-13 
illustrates  this  point. 


Figure  3-13.  Representation  of  processes  on  the  psychrometric  chart  showing  where  the 
transition  to  a supersaturated  state  occurs  by  keeping  the  same  entering  air 
temperature  (-8.3°C)  and  gradually  increasing  the  entering  air  relative 
humidity 
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While  the  straight-line  theory  and  the  calculated  path  approach  agree  reasonably 
well  at  the  entrance  and  exit  of  the  coil,  significant  differences  between  the  two  methods 
exist  in  the  middle  section  of  the  coil.  For  example,  Figure  3-13  shows  a 34%  deviation 
in  the  percentage  of  cooling  calculated  for  Scenario  2 employing  both  methods.  This  is 
manifested  in  the  differences  between  paths  la-3a  and  la-3a  ’,  where,  in  the  former,  the 
straight-line  theory  predicts  84%  cooling  achieved  between  Rows  5 and  6,  whereas  the 
corresponding  value  for  Path  la-3a  ’ is  only  50%. 

Conclusions 

This  chapter  described  a calculation  procedure  as  well  as  an  experiment  for 
determining  the  air  path  on  an  actual  industrial-size  finned-tube,  multi-row  coil  operating 
under  frosting  conditions.  A calculated  path  based  on  the  analysis  presented  and  a path 
based  on  the  straight-line  theory  of  applied  psychrometrics  were  correlated  with  the 
prevailing  psychrometric  conditions  in  the  freezer.  Determining  the  calculated  path 
enabled  accurate  calculation  of  the  conditions  leaving  a given  row  in  a multi-row  coil. 
This,  in  turn,  helped  identify  the  coil  location  in  the  vicinity  of  which  the  transformation 
from  the  subsaturated  zone  to  the  supersaturated  zone  occurred.  This  is  crucial  in 
identifying  a demarcation  line  between  the  unfavorable  snow-like  frost  and  the  more 
traditional  frost-formation  patterns.  Calculated  results  were  verified  against  experimental 
results  at  the  entrance  and  exit  of  the  coil  with  typical  differences  falling  well  within  the 
uncertainty  ranges  of  both  temperature  and  humidity.  In  the  middle  section  of  the  coil, 
on  the  other  hand,  significant  deviations  (up  to  34%)  between  the  calculated  and  the 
straight-line  path  were  observed.  This  underscores  the  significance  of  using  the  method 
outlined  in  this  chapter  if  reasonably  accurate  prediction  of  the  location  of  the 


demarcation  line  between  snow-like  and  conventional  frost  is  to  be  determined. 
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According  to  results  pertaining  to  the  location  of  the  demarcation  line,  transition  to  a 
supersaturated  state  occurs  quicker  and  at  a more  upstream  location  of  the  coil  as  the 
relative  humidity  of  the  entering  air  increases. 


CHAPTER  4 

ISSUES  PERTAINING  TO  HOT-GAS  COIL  DEFROSTING 

Objectives 

Previous  coil  defrosting  studies  have  largely  neglected  the  distinctions  between 
two  types  of  frost-formation  patterns.  As  was  mentioned  before,  the  frost  forming  near 
saturated  entering-air  conditions  is  more  snow-like  in  appearance  and  contributes  more 
adversely  to  coil  heat  transfer  performance  than  frost  forming  under  predominantly 
subsaturated  conditions.  While  the  mechanisms  of  formation  of  both  types  of  frost  are 
similar  in  some  aspects  and  different  in  others,  the  implications  on  the  defrosting  process 
are  markedly  different.  It  is  thus  the  objective  of  this  chapter  to  investigate  the  issues 
pertaining  to  the  coil  defrosting  process  using  the  hot-gas  refrigerant  method  in  order  to 
provide  a deeper  insight  and  a better  understanding  of  the  problem  in  question.  This 
method  of  defrosting  was  selected  because  it  is  the  most  widely  used  defrosting  technique 
in  industrial  freezers. 

Since  available  coil  defrosting  studies  lack  data  on  coils  operating  near  saturated 
entering-air  conditions,  and  since  this  seems  to  be  a prevalent  (but  undesirable)  mode  of 
operation  in  industrial  refrigeration  applications,  it  was  natural  to  investigate  this 
mechanism  further.  This  part  of  the  study  includes  both  physical  and  infrared  analysis  of 
a coil  being  defrosted.  The  infrared  portion  of  the  analysis  will  provide  the  spatial  and 
temporal  distributions  of  temperature  on  the  coil  surface,  a key  variable  in  understanding 
the  defrosting  dynamics.  Throughout  the  study,  an  attempt  will  also  be  made  to  compare 
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the  performance  of  the  defrosted  coil  with  those  of  similar  coils  operating  under 
subsaturated  conditions  (e.g.  Stoecker  et  al.,  1983). 


A procedure  used  to  calculate  the  defrost  efficiency  is  described  in  this  section. 
Defrosting  heat  input  rate. 


where  Qdef  in  is  the  defrosting  heat  input  rate  (kW),  mR  vap  is  the  refrigerant  mass  flow  rate 
(kg/s),  iin  is  the  refrigerant  inlet  enthalpy  (kJ/kg),  and  iout  is  the  refrigerant  outlet 
enthalpy  (kJ/kg).  Defrosting  heat  input  rates  computed  employing  Equation  4-1  have  a 
maximum  uncertainty  of  ±4%  . 

Defrosting  heat  input. 


where  Qdefin  is  the  defrosting  heat  input  (kJ)  and  rdef  is  the  defrosting  time  (s).  Defrosting 
heat  inputs  computed  employing  Equation  4-2  have  a maximum  uncertainty  of  ±4%. 

Refrigerant  vapor  mass  flow  rate.  To  calculate  the  defrosting  heat  input  it  is 
necessary  to  determine  the  mass  flow  rate  of  the  refrigerant.  The  following  procedure  is 
used  to  calculate  the  vapor  mass  flow  rate  of  the  refrigerant: 

a.  Calculate  the  volumetric  flow  rate  ( VR ) from  the  measured  value  of  the 
differential  pressure  across  an  orifice  meter  located  on  the  hot-gas  defrosting  line. 
Detailed  information  regarding  the  shape,  dimensions,  and  performance  of  the 


Analysis 


(4-1) 


(4-2) 
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orifice  meter  is  provided  by  Sherif  et  al.  (1997).  To  get  the  volumetric  flow  rate  a 
fifth  order  polynomial  equation  is  used 

VR  = 4. 187657x1  O'5  +(6.3786  lxl  O'4  Ap)- (2.6509  lxl  0‘4AP2)  (4‘3) 

+ (7.27778x  1 O'5  AP3 )-  (l  .0033x1  O'5  AP4 )+  (5 .3 1 96x1  O'7  A P5 ) 

where  VR  is  the  volumetric  flow  rate  (m3/s),  and  AP  is  the  differential  pressure 

across  the  orifice  meter  in  inches  of  water  column  differential  (inches  WCD) 

b.  Calculate  the  local  average  velocity  (umR ) of  the  refrigerant  vapor  upstream  of 

the  orifice  using  the  following  equation: 


Um,R  = 


VD 


1 orifice 


(4-4) 


2 

where  Aorifice  is  the  inside  cross  sectional  area  upstream  of  the  orifice  (m  ),  and 

where  the  inside  diameter  (f>,)of  the  pipe  is  13.39  mm. 
c.  From  the  properties  of  R-22  refrigerant  vapor,  calculate  the  Reynolds  number 
(Re)  as  follows: 

Re  = (4'5) 
Mr 


where  pR  is  the  density  of  the  refrigerant  vapor  (kg/m3)  and  pR  is  the  dynamic 
viscosity  of  the  refrigerant  vapor  [kg/(m  s)]. 
d.  Calculate  the  dimensionless  discharge  coefficient  ( Cd ) of  the  orifice  meter  using 
the  following  equation  (ASME  1991,  Miller  1981): 

2,  a 1 o/t, .,8  , a mn  V*  a„,«0,51V 


Cd  = 0.5959  + 0.0312^/  —0.1 84^  +0.039 


1-V/4 


0.0158^/  +- 


Re 


0.75 


(4-6) 
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where  (//  is  the  ratio  of  the  orifice  throat  diameter  to  the  pipe  inside  diameter  and 
is  equal  to  0.5. 

e.  Calculate  the  dimensionless  flow  coefficient  [Korifice)  of  the  orifice  meter  using  the 


following  equation  (Fox  and  McDonald  1992,  Miller  1981): 

A . - ' 


(4-7) 


orifice 


Vi7 


v 


f.  Calculate  the  mass  flow  rate  ( mR  vap ) of  the  refrigerant  vapor  using  the  following 


equation  (ASHRAE  2001  and  Fox  and  McDonald  1992): 


™R,vap  K orifice  ^orifice  P ^^SPr  P water  ^ 


(4-8) 


where  mR  vap  is  the  mass  flow  rate  (kg/s),  pwater  and  pR  are  the  density  of  water 

and  refrigerant  in  (kg/m3),  respectively,  A P is  the  differential  pressure  across  the 
orifice  meter  in  meters  of  water  column  differential,  and  g is  the  acceleration  of 
gravity  (9.81  m/sl 2). 

The  uncertainties  of  the  volumetric  flow  rate  (k^),  the  local  average  velocity  (umR ), 
the  Reynolds  number  ( Re ),  the  discharge  coefficient {Cd),  the  flow  coefficient  (k  onfic Jof 
the  orifice  meter,  and  the  refrigerant  vapor  mass  flow  rate  (mR  vap ) are,  respectively,  ±2%, 

±2%,  ±3.5%,  ±1%,  ±1%,  and  ±2.5%. 

Rate  of  Heat  Removed  with  the  FCU-T  Melt. 


Qmelt 


m 


ijliq  h ce  ) 


melt  \ liq 


(4-9) 


ldef 


where  Qmell  is  the  rate  of  heat  removed  with  the  FCU-T  melt  (kW),  mmelt  is  the  mass  of  the 
melt  during  the  defrosting  process  (kg),  iUq  is  the  enthalpy  of  liquid  water  at  the  leaving 
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melt  temperature  (kJ/kg),  and  iice  is  the  enthalpy  of  ice  at  the  coil  surface  temperature 

(which  is  equal  to  the  temperature  of  the  frost  at  the  beginning  of  the  defrosting  process) 
(kJ/kg).  The  rates  of  heat  removed  with  the  FCU-T  melt  computed  employing  Equation 
4-9  have  a maximum  uncertainty  of  ±4%. 


Defrost  efficiency. 

The  defrost  efficiency  is  defined  as  the  ratio  between  the  rate  of  heat  removed 
with  the  FCU-T  melt  and  the  total  heat  input  rate  for  the  defrosting  process 


(4-10) 


Defrost  efficiencies  computed  using  Equation  4-10  have  a maximum  uncertainty 


of  ±4% 


There  are  two  parameters  of  relevance  to  this  defrosting  problem;  the  average  coil 
performance  load  rate,  and  the  fan-defrost  contribution.  A procedure  used  to  calculate, 
the  average  coil  performance  load  rate,  and  the  fan-defrost  contribution  is  described  in 
this  section. 

Average  coil  performance  load  rate.  The  average  coil  performance  load  rate 
( ACPLR ) is  defined  as  the  ratio  of  the  sum  of  the  total  refrigeration  load  and  the  defrost 
heat  load  divided  by  the  total  combined  time  of  the  FCU-T  refrigeration  and  defrosting 
modes.  The  following  procedure  is  used  to  calculate  the  average  coil  performance  load 
rate,  which  has  a maximum  uncertainty  of  ±4%: 

a.  Calculate  the  total  refrigeration  load  (Qlotal) in  kJ  using  the  following  equation: 

Qlotal  ( Qlotal tS  Qlotal  ref 


(4-11) 
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where  Q,otal<s  is  the  total  sensible  heat  rate  gain  (kW),  QlolalJ  is  the  total  latent  heat 
rate  gain  (kW),  and  zref  is  the  duration  of  the  FCU-T  refrigeration  mode  (s). 

b.  Calculate  the  defrost  heat  load  ( QdefJoad ) in  kJ  using  the  following  equation: 

Qdef  .load  ^Qdeffn  Qmelt ) ^ def  ^ ^ 

where  zdef  is  the  duration  of  the  FCU-T  hot-gas  defrosting  mode  (s). 

c.  Calculate  the  total  time  duration  (ztotal  )for  a complete  FCU-T  refrigeration  and 


defrosting  cycle  using  the  following  equation: 


T total  ~ ^def  + Tref 


(4-13) 


d.  Calculate  the  average  coil  performance  load  rate  (ACPLR)  in  kW  using  the 
following  equation: 


^ Qtotal  QdefJoad 


(4-14) 


v total 


Fan  defrost  contribution  IFDQ.  In  many  refrigeration  systems,  the  theoretical 
contribution  of  the  fan  and  defrost  heat  loads  to  the  net  refrigeration  load  exceeds  1 5%  of 
the  net  refrigeration  load  (ASHRAE  1998),  where  the  net  refrigeration  load  equals  the 
total  refrigeration  load  less  the  fan  heat.  Considering  this  idea,  the  ratio  of  the  fan  and 
defrost  heat  loads  to  the  net  refrigeration  load  has  been  utilized  here  in  this  investigation 
for  comparative  analysis  among  the  various  FCU-T  refrigeration/defrost  cycles.  The 
following  procedure  is  used  to  calculate  the  Fan-Defrost  Contribution,  which  has  a 


maximum  uncertainty  of  ±4%: 
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a.  Calculate  the  fan  heat  (Qfan)  using  the  following  equation: 


fan  Q fan  ^ ref 


(4-15) 


where  Qfcn  is  the  fan  power  (kW)  and  xref  is  the  duration  of  the  FCU-T 
refrigeration  mode  (s). 


where  Qlolal  is  the  total  refrigeration  load  calculated  (kJ). 

c.  Calculate  the  fan-defrost  contribution  ( FDC)  in  (%)  using  the  following 
equation: 


As  was  stated  earlier,  one  of  the  objectives  of  this  section  was  to  examine  the 
dynamics  of  defrosting  freezer  coils  using  the  hot-gas  refrigerant  method  for  cases  in 
which  the  coil  entering-air  conditions  are  in  the  vicinity  of  the  saturation  curve  of  the 
psychrometric  chart.  In  order  to  achieve  this  objective,  both  infrared  and  physical  images 
of  the  coil  being  defrosted  are  recorded  at  successive  time  increments.  In  order  to 
provide  some  confidence  in  the  authenticity  of  the  data  collected,  a series  of 
qualifications  tests  have  been  performed.  In  these  tests,  key  defrost  variables  have  been 
measured  and  compared  with  corresponding  variables  found  in  the  open  literature  of 
similar  coils.  These  tests  are  described  and  discussed  below.  It  should  be  noted  here, 


b.  Calculate  the  net  refrigeration  load  (Q„el)i n kJ  using  the  following  equation: 


(4-16) 


(4-17) 


Results  and  Discussion 


however,  that  existing  data  on  coil  defrosting  in  the  open  literature  are  available  only  for 
conventionally-frosted  coils  operating  under  subsaturated  conditions.  So,  while  the 
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comparison  is  not  an  exact  one,  analysis  of  the  general  trends  could  still  prove  useful  and 
insightful. 

Comparisons  with  and  verification  of  existing  defrosting  data.  Figure  4-1  shows 
the  temperature-time  history  of  the  coil  refrigerant  temperatures  at  the  inlet  and  outlet  for 
the  present  experiment  along  with  the  experiments  reported  by  Stoecker  et  al.  (1983).  As 
can  be  seen  in  the  figure,  the  hot  gas  leaving  the  FCU-T  pan  then  enters  the  FCU-T  coil 
at  a lower  temperature  due  to  the  loss  of  heat  when  the  hot  gas  mixes  with  the  liquid 
refrigerant  at  and  inside  the  FCU-T  coil.  The  defrost  refrigerant  leaves  the  FCU-T  coil  at 
a lower  temperature  due  to  the  mixing  with  the  liquid  refrigerant  inside  the  FCU-T  and 
also  due  to  the  refrigerant  heat  lost  to  melt  the  frost  accumulated  in  the  FCU-T  coil.  As 
can  be  seen,  the  trends  are  the  same  for  both  Stoecker’ s data  and  the  present  data, 
however  the  numerical  values  are  different  since  the  conditions  used  for  the  experiment 
are  not  the  same  as  those  of  Stoecker  et  al.  (1983). 

A similar  comparison  has  also  been  conducted  with  Stoecker’ s mass  flow  rate 
data  as  shown  in  Figure  4-2  and  with  Stoecker’s  pressure  data  as  shown  in  Figure  4-3. 
The  mass  flow  rate  is  high  at  the  beginning  of  the  defrosting  process,  but  it  falls  off 
sharply  during  the  first  two  minutes.  Shortly  afterward,  the  mass  flow  rate  remains 
almost  constant  until  the  conclusion  of  the  defrosting  process.  The  trends  observed  in  the 
mass  flow  rate  data  agree  with  those  trends  presented  by  Stoecker  et  al.  (1983). 

In  Figure  4-3,  the  refrigerant  pressures  at  the  inlet  and  outlet  of  the  coil  during 
defrosting  have  been  measured.  From  the  trends  observed,  the  inlet  pressure  increases  at 
the  beginning  of  the  defrosting  process  until  it  reaches  the  value  of  the  inlet  pressure 
regulator.  The  difference  between  the  inlet  and  outlet  pressures  is  due  to  the  pressure 
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Figure  4-1 . Refrigerant  temperatures  at  the  inlet  and  outlet  of  the  coil  during  defrosting 
for  both  Stoecker  et  al.  (1983)  and  present  work 


Figure  4-2.  Mass  flow  rate  of  refrigerant  during  defrosting  for  both  Stoecker  et  al.  (1983) 
and  present  work 
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drop  across  the  coil.  Again,  the  trends  are  the  same  as  those  of  Stoecker  et  al.  (1983) 
even  though  the  numerical  values  of  pressures  are  different. 
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Figure  4-3.  Refrigerant  pressures  at  the  inlet  and  outlet  of  the  coil  during  defrosting  for 
both  Stoecker  et  al.  (1983)  and  present  work 

Defrosting  snow-like  frost  versus  conventional  frost.  In  order  to  experimentally 
investigate  the  complete  refrigeration-defrost  (R/D)  cycle  when  the  freezer  operates 
under  subsaturated  and  near  saturated  entering-air  conditions,  two  tests  were  performed 
under  the  same  initial  freezer  conditions,  for  5 hours  during  refrigeration  and  1 5 minutes 
during  defrosting.  Initial  freezer  air  conditions  were  £>57’=-15°C  and  RH=  75%,  with  an 
apparatus  dew-point  of  -22°C.  The  coil  face  velocity  for  both  tests  was  kept  constant  at 


—♦—Present work  (Inlet) 

— • — Present  work  (Outlet) 

Stoecker  etal.  (1983)  (Inlet) 

Stoecker  etal.  (1983)  (Outlet) 

0 2.5  5 7.5  10  12.5  1 

Defrosting  time  (min) 


3.7  m/s. 
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Case  a.  In  this  case  refrigeration  was  performed  near  saturated  entering-air 
conditions.  During  the  five  hours  of  refrigeration  0.4  kg/h  of  steam  were  injected.  Figure 
4-4  shows  the  room  and  coil  surface  temperatures  for  a complete  R/D  cycle.  The  room 
temperature  during  the  refrigeration  process  was  between  -12°C  and  -14°C  for 
approximately  four  hours.  As  frost  started  to  block  the  coil,  the  room  temperature  started 
to  increase  as  the  capacity  of  the  coil  to  refrigerate  was  beginning  to  decline.  At  that 
point,  the  freezer  air  temperature  has  risen  to  5°C.  A visual  image  of  the  frosted  coil  at 
this  stage  is  shown  in  Figure  4-5.  The  frost  thickness  before  starting  the  defrosting 
process  was  5 mm.  Figure  4-6  shows  an  image  of  the  coil  five  minutes  into  the 
defrosting  process.  In  this  figure,  chunks  of  frost  can  be  seen  falling  from  the  top  instead 
of  melting  and  collecting  into  the  pan  located  at  the  bottom  of  the  coil.  Figure  4-7  shows 
an  image  of  the  coil  ten  minutes  into  the  defrosting  process.  In  this  figure,  chunks  of 
frost  continue  to  fall  down  from  the  top  of  the  coil.  It  is  noteworthy  to  observe  how 
difficult  the  frost  is  to  remove.  Fifteen  minutes  into  defrosting,  the  coil  can  still  be  seen 
with  about  half  the  amount  of  frost  present  (Figure  4-8).  The  amount  of  melt  collected 
after  performing  the  defrosting  process  was  0.36  kg.  The  total  heat  input  into  the 
defrosting  process  was  649.5  kJ  and  the  total  heat  used  in  the  melting  process  was  141.7 
kJ,  yielding  a defrost  efficiency  of  21.8%.  This  is  a rather  poor  defrost  efficiency,  which 
can  be  blamed  on  the  small  amount  of  melt  collected  and  the  corresponding  large  amount 
of  evaporated  melt  and  sublimated  frost. 
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Figure  4-4.  Freezer  air  and  coil  surface  temperatures  during  refrigeration  and  defrosting 


Figure  4-5.  Coil  condition  after  5 hours  of  refrigeration  under  supersaturated  conditions 
in  the  refrigeration  mode  (same  as  beginning  of  defrosting) 
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Figure  4-6.  Coil  after  5 minutes  into  the  defrosting  process 


Figure  4-7.  Coil  after  10  minutes  into  the  defrosting  process 
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Figure  4-8.  Coil  after  15  minutes  into  the  defrosting  process 

Case  b.  In  this  case  the  refrigeration  process  was  conducted  under  subsaturated 
conditions  for  the  same  length  of  time  as  the  previous  experiment.  Figure  4-9  displays 
both  the  freezer  air  and  coil  surface  temperatures  during  the  refrigeration  portion  of  the 
R/D  cycle.  Freezer  air  temperature  was  approximately  constant  at  -15°C,  while  the  coil 
surface  temperature  was  approximately  constant  at  -22°C.  Defrosting  was  not  conducted 
in  this  case  as  very  little  frost  could  be  detected  after  the  lapse  of  five  hours  in  the 
refrigeration  mode  (see  Figure  4-10).  This  is  in  stark  contrast  to  Case  a and  highlights 
the  need  to  keep  freezer  air  in  a subsaturated  condition. 

As  can  be  seen  from  Figures  4-5  and  4-10,  operating  freezers  near  saturated  coil 
entering-air  conditions  produces  a more  dense,  snow-like  frost  capable  of  blocking  the 
coil  in  relatively  short  time  periods  and  contributing  adversely  to  the  coil  heat  transfer 


performance. 
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Figure  4-9.  Freezer  and  coil  surface  temperatures  during  refrigeration 


Figure  4-10.  Coil  after  5 hours  of  operation  under  subsaturated  conditions  during  the 
refrigeration  process 


60 


Infrared  analysis  of  the  defrosting  process.  In  order  to  provide  another  insight  of 
the  changes  taking  place  during  defrosting,  infrared  images  are  shown  as  well  (see 
Figures  4-11  through  4-15).  These  images  are  intended  to  show  how  the  frost  surface 
temperature  changes  with  time  as  the  defrost  heat  penetrates  the  frost  layer  on  the  coil.  In 
order  to  get  the  numerical  values  of  the  frost  surface  temperature  at  key  locations,  five 
spots  (Tl,  T2,  T3,  T4  and  T5)  were  identified. 

The  average  frost  surface  temperature  before  the  beginning  of  the  defrosting 
process  was  -13.3°C  (Figure  4-11).  One  minute  into  defrosting  (Figure  4-12),  the 
temperature  of  the  frost  surface  at  the  top  left  comer  increases  to  -2.9°C.  Three  minutes 
into  defrosting  (Figure  4-13),  the  top  of  the  coil  is  seen  to  have  an  average  temperature  of 
-2.5°C,  while  the  rest  of  the  frost  has  an  average  temperature  of  8.6°C.  After  eight 
minutes  (Figure  4-14),  frost  at  the  top  part  of  the  coil  had  already  melted.  In  this  figure 
one  can  observe  the  temperature  of  the  metal  surface  of  the  coil.  This  temperature 
changes  from  26.4°C  at  the  top  left  comer  to  1 1 .7°C  at  the  top  right  comer,  while  the 
average  frost  surface  temperature  for  the  remainder  of  the  surface  is  only  -2.1°C. 

In  order  to  create  a more  comprehensive  picture  of  the  quantitative  variation  in 
the  frost  surface  temperature  during  defrosting,  the  temperature-time  history  of  all  five 
spots  selected  is  plotted  in  Figure  4-15.  In  this  figure  one  can  observe  that  the  spots 
located  at  the  top  of  the  coil  (Tl  and  T2)  reach  higher  temperatures  faster  than  the  others. 
Obviously,  this  thermal  picture  is  consistent  with  the  images  of  the  physical  coil  while 
being  defrosted  shown  earlier  in  this  investigaion.  Five  minutes  into  defrosting  (Figure 
4-16),  one  can  observe  that  Spots  Tl  and  T2  now  represent  the  temperatures  of  the  bare 


coil  surface. 
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Figure  4-11.  Frost  surface  temperature  before  the  beginning  of  the  defrosting  process 


6.7°C 


-10 


-1  9.6°C 


Figure  4-12.  Frost  surface  temperature  after  1 minute  into  defrosting 
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Figure  4-13.  Frost  surface  temperature  after  3 minutes  into  defrosting 
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Figure  4-14.  Frost  surface  temperature  after  8 minutes  into  defrosting 
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Figure  4-15.  Frost  surface  temperature  vs.  time  during  defrosting 
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Figure  4-16.  Frost  and  pipe  surface  temperatures  vs.  time  during  defrosting 
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Conclusions 

An  experimental  study  of  a finned-tube  freezer  coil  during  hot-gas  defrosting  was 
conducted.  The  study  included  quantitative  comparisons  with  similar  coils  as  well  as 
both  quantitative  and  qualitative  evaluations  of  the  spatial  and  temporal  distributions  of 
both  the  frost  surface  temperature  and  the  frost  mass.  Throughout  the  study,  comparisons 
between  two  distinct  types  of  frost  were  made  both  in  the  formation  pattern  and  in  the 
manner  with  which  they  melt  during  the  defrosting  process.  Thermal  images  taken  of  the 
coil  while  being  defrosted  are  consistent  with  the  physical  images  of  the  defrosted  coil. 
Comparisons  with  existing  defrosting  data  were  presented  in  order  to  provide  some 
confidence  in  the  experimental  results  presented. 


CHAPTER  5 

ALTERNATIVE  ENERGY  CONSERVATION  TECHNOLOGIES  IN  INDUSTRIAL 
FREEZERS  OPERATING  UNDER  FROSTING  CONDITIONS 

Objective 

Operating  freezers  near  saturated  coil  entering-air  conditions  produce  a more 
dense,  snow-like  frost  capable  of  blocking  the  coil  in  relatively  short  time  periods  and 
contributing  adversely  to  the  coil  heat  transfer  performance.  While  every  effort  should  be 
made  to  avoid  operating  the  freezers  and  refrigerated  facilities  under  theses  adverse 
conditions,  some  technologies  may  be  available  to  lessen  the  negative  impact  of  such 
conditions  once  operating  in  this  mode  becomes  inevitable.  This  chapter  will  particularly 
target  three  distinct  techniques;  namely,  the  use  of  heat-pipe  assisted  heat  exchangers, 
employing  dampered  coils,  and  operating  freezer  coils  at  high  face  velocities.  The 
feasibility  of  these  technologies  will  be  explored  both  experimentally  and  analytically. 
The  results  will  be  supported  by  visual  observations  and  images  taken  inside  the 
laboratory  freezer  under  study. 

Feasibility  of  Heat  Pipes 

One  of  the  technologies  is  the  use  of  heat-pipe  assisted  freezer  coils  as  a way  to 
reduce  the  amount  of  energy  required  during  the  refrigeration  portion  of  the 
refrigeration/defrost  (R/D)  cycle  of  operation.  The  main  idea  here  is  to  let  the  sensible 
cooling  be  performed  by  one  leg  of  a heat-pipe  heat  exchanger  while  letting  the 
dehumidification  part  of  the  refrigeration  process  be  performed  by  the  actual  freezer  coil. 
This  may  provide  some  energy  savings  potential.  In  this  section  a discussion  of  the  use 
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of  heat  pipes  for  dehumidification  will  be  presented,  followed  by  performing  a 
comparative  analysis  with  and  without  heat  pipes. 

Heat  Pipes  for  Dehumidification 

A heat  pipe  consists  of  a sealed  container  (made  of  aluminum  or  copper)  whose 
internal  surfaces  have  a capillary  wicking  material.  This  container  has  liquid  inside 
(under  its  own  pressure)  that  enters  the  pores  of  the  capillary  material,  wetting  all  internal 
surfaces.  When  heat  is  applied  at  any  point  of  the  heat  pipe  the  liquid  at  that  point  boils 
and  changes  to  vapor.  Then,  the  liquid  picks  up  the  latent  heat  of  vaporization.  The  gas 
having  a higher  pressure,  moves  inside  the  sealed  container  to  a colder  location  where  it 
condenses.  Thus,  the  gas  gives  up  the  latent  heat  of  vaporization  and  moves  heat  from 
the  input  to  the  output  end  of  the  heat  pipe.  The  difference  between  heat  pipes  and 
conventional  heat  exchangers  is  the  capacity  to  transfer  energy.  For  example,  a heat  pipe 
is  capable  of  transferring  1000  times  more  energy  than  a solid  copper  bar,  the  best-known 
commercial  conductor,  without  moving  parts,  external  fans,  connections  to  a fluid  supply, 
noise,  and  external  power. 

In  its  in-line  (wrap  around)  configuration,  the  heat  pipe  wraps  around  the 
dehumidification  coil,  with  one  section  of  the  heat  pipe  coil  upstream  and  one  section 
downstream  (Figure  5-1).  Hot  and  humid  air  enters  through  the  first  heat  pipe  coil 
section  where  it  is  precooled  prior  to  entering  the  dehumidification  coil.  This  shifts  the 
burden  of  sensible  cooling  to  the  upstream  leg  of  the  heat  pipe  coil  instead  of  the 
dehumidification  coil.  Dehumidified  air  emerges  from  the  cooling  coil,  passes  through 
the  second  section  of  the  heat  pipe  coil,  where  it  is  reheated  to  a temperature  with  a more 
comfortable  relative  humidity.  Precooling  and  reheating  by  the  heat  pipe  are  performed 
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without  any  expenditure  of  external  energy.  The  two  sections  of  the  heat  pipe  coil  simply 
exchange  precooling  energy  in  the  first  section  for  reheating  energy  in  the  second  section. 
A psychrometric  depiction  of  dehumidification  processes  using  heat  pipes  is  shown  in 
Figure  5-2. 


Heat  Pipe  Cooling  Coil  Heat  Pipe 


Figure  5-1 . Heat  pipe  technology  used  in  dehumidification  processes 
Calculation  Software 

In  order  to  simulate  the  performance  of  the  heat  pipe-assisted  coil,  a calculation 
software,  called  HeatpipePRO,  was  used  (Nelson  1999  private  communications).  This 
software  was  developed  by  Colmac  Coil  Manufacturing,  Inc.  Figure  5-3  shows  the  main 
screen  of  the  software. 

The  data  needed  to  operate  this  software  are  the  inlet  air  dry-bulb  temperature  and 
either  the  wet-bulb  temperature  or  the  relative  humidity,  the  altitude,  the  volumetric 
supply  air  flow  rate  in  cfm  (also,  either  scfm  for  standard  conditions  or  acfm  for  altitude 
effects  can  be  used),  exhaust  air  dry-bulb  temperature  and  either  the  wet-bulb 
temperature  or  the  relative  humidity,  and  the  volumetric  exhaust  airflow  rate  in  cfm. 
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Figure  5-2.  Representation  of  the  dehumidification  process  using  heat  pipes 


Figure  5-3.  Main  screen  of  the  FieatpipePRO  vl.O  software 
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Analysis  of  the  Dehumidification  Process 

In  order  to  demonstrate  the  advantages  of  using  heat  pipes,  data  from  an 
experimental  test  are  utilized  in  a computer  code  that  simulates  the  performance  of  wrap- 
around heat  pipe  heat  exchangers.  The  entering  air  conditions  for  the  test  were 
DBT-  3°C  and  RH=  50%,  while  the  apparatus  dew-point  temperature  was  ADP=-  11°C. 
The  volumetric  airflow  rate  was  90  m /min  and  the  room  sensible  heat  ratio  0.81.  The 
leaving  air  conditions,  on  the  other  hand,  were  DBT-  -8.2°C  and  RH=  88%.  The  cooling 
and  dehumidification  processes  with  and  without  heat  pipes  are  shown  in  Figure  5-4. 


Figure  5-4.  Coil  and  dehumidification  processes  with  and  without  heat  pipes 
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Line  1-3  represents  the  cooling  and  dehumidification  process  without  heat  pipes 
and  Lines  1-2-3-4  represent  the  processes  using  heat  pipes.  For  the  cooling  and 
dehumidification  process  using  heat  pipes,  Line  1-2  represents  the  precooling  heat  pipe, 
Line  2-3  represents  the  cooling  process  through  the  coil,  Line  3-4  represents  the  reheating 
heat  pipe,  and  Line  4-1  represents  the  supply  air  warm-up  path. 

The  total  energy  required  for  the  cooling  and  dehumidification  process  without  a 
heat  pipe  (Line  1-3 ) is  24.94  kW.  However,  the  total  energy  required  for  the  cooling  and 
dehumidification  process  with  a wrap-around  heat  pipe  (Line  2-3 ) is  18.71  kW.  The 
energy  savings  using  a wrap-around  heat  pipe  are  represented  by  the  energy  associated 
with  Line  1-2  (precooling  process).  This  energy  has  a value  of  6.23  kW.  The  total 
energy  required  for  the  cooling  and  dehumidification  process  decreases  by  25%  when  a 
wrap-around  heat  pipe  is  used. 

The  same  analysis  can  be  conducted  if  two  or  more  heat-pipes  of  the  same  model 
are  used  in  parallel.  This  arrangement  corresponds  to  stacked  heat  pipes,  where  the 
supply/exhaust  air  flow  rates  are  distributed  over  the  total  amount  of  heat  pipe  surface 
area.  In  this  case,  the  energy  associated  with  Line  2-3  (cooling  process  through  the  coil) 
will  be  16.8  kW,  while  that  associated  with  Line  1-2  (precooling  process)  will  be  8.14 
kW  and  that  associated  with  Line  3-4  (reheating  process)  will  be  8.2  kW.  Furthermore, 
the  room  sensible  heat  ratio  will  be  0.79  and  the  heat  pipe  effectiveness  will  be  39.1%.  If 
three  wrap-around  heat  pipe  units  are  used  in  parallel,  on  the  other  hand,  the  energy 
associated  with  Line  2-3  (cooling  process  through  the  coil)  becomes  15.92  kW,  while  that 
associated  with  Line  1-2  (precooling  process)  will  be  9.02  kW.  Also,  the  energy 
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associated  with  Line  3-4  (reheating  process)  will  be  9.1  kW,  with  a room  sensible  heat 
ratio  of  0.78  and  a heat  pipe  effectiveness  of  43.2%. 

According  to  the  above  results,  the  heat  pipe  effectiveness  increases  and  the 
energy  required  for  cooling  and  dehumidification  decreases  when  two  or  more  wrap- 
around heat  pipes  are  used  in  parallel.  However,  an  important  parameter  to  consider  here 
is  cost.  The  cost  of  two  units  in  parallel  is  approximately  1.8  times  the  cost  of  a single 
heat  pipe  unit;  while,  the  cost  of  three  units  in  parallel  is  approximately  2.5  times  the  cost 
of  a single  unit. 

When  wrap-around  heat  pipes  have  to  be  selected  for  any  application,  the  savings 
in  energy  as  well  as  the  financial  savings  have  to  be  taken  into  consideration  in  order  to 
make  the  best  selection  for  the  required  application. 

Analysis  of  the  Energy  Requirements 

The  total  energy  required  for  the  cooling  and  dehumidification  processes  with  and 
without  heat  pipes  as  a function  of  the  entering  air  temperature  is  shown  in  Figure  5-5, 
for  two  relative  humidity  values.  For  the  results  shown  in  this  figure,  the  following 
conditions  apply:  ADP=-  11°C,  volumetric  airflow  rate=90  m /min.  The  wrap-around 
heat  pipe  used  was  the  same  model  employed  before.  The  advantage  of  using  heat  pipes 
in  terms  of  the  energy  required  for  cooling  and  dehumidification  is  illustrated  in 
Figure  5-5.  From  this  figure  it  can  be  seen  that  the  energy  required  for  cooling  and 
dehumidification  (with  or  without  heat  pipes)  increases  when  the  relative  humidity 


increases. 
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Figure  5-5.  Energy  required  for  cooling  and  dehumidification  with  and  without  heat 
pipes  for  two  different  relative  humidities 

Feasibility  of  Using  Pampered  Coils  during  the  Defrosting  Process 
As  was  stated  earlier,  this  study  is  concerned  with  the  feasibility  of  using 
dampered  fan-coil  units  and  their  potential  in  energy  conservation  by  reducing  the 
defrosting  frequency.  In  order  to  achieve  this  goal,  the  defrost  efficiency  and  the  amount 
of  melt  collected  at  the  end  of  the  defrosting  process  were  computed  for  a given  set  of 
freezer  conditions.  Typically,  when  an  undampered  coil  is  defrosted,  part  of  the 
defrosting  heat  input  goes  into  the  freezer  space  in  the  form  of  evaporated  water  and 
sublimated  frost  instead  of  being  used  to  melt  the  frost.  However,  if  dampered  coils  are 
used,  the  heat  that  goes  into  the  freezer  space  can  be  minimized  and  more  heat  can  be 
applied  to  the  actual  melting  of  the  frost.  This  will  necessarily  result  in  an  increase  in  the 
defrost  efficiency  and  the  subsequent  reduction  in  the  frequency  of  coil  defrosting. 
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The  results  presented  in  this  section  are  significant  for  freezers  with  more  than 
one  cooling  coil.  This  is  because  when  one  coil  is  in  the  defrosting  mode;  all  others  are 
usually  in  the  refrigeration  mode.  Thus,  a defrosting  process  with  high  efficiency,  as 
manifested  by  a larger  amount  of  melt  and  a smaller  amount  of  evaporation  and 
sublimation,  will  result  in  a smaller  burden  on  all  other  coils,  which  may  be 
simultaneously  operating  in  the  refrigeration  mode  inside  the  freezer.  Any  additional 
moisture  generated  by  the  coil  being  defrosted  will  immediately  deposit  on  the  remaining 
coils  as  undesirable  frost  and  will  necessarily  increase  the  frequency  of  their  defrosting 
cycle.  If  dampered  coils  are  used,  on  the  other  hand,  the  additional  moisture  generated 
from  the  coil  being  defrosted  will  be  trapped  by  the  dampers  and  will  eventually  collect 
into  the  melt  pan  and  be  discharged  outside  the  freezer.  Also,  when  a significant  amount 
of  moisture  is  introduced  into  the  freezer  space  due  to  improper  defrosting  practices,  it  is 
likely  that  freezer  air  may  become  supersaturated.  Formations  resulting  from  exposure  to 
supersaturated  air  have  been  shown  by  Sherif  et  al.  (2001)  to  be  snow-like  in  crystal 
structure  and  much  more  extensive  in  amount  than  traditional  formations  when  the  air  is 
subsaturated. 

Results  and  discussion.  In  order  to  analyze  the  performance  of  dampered  coils, 
either  or  both  the  coil  suction  and  discharge  were  covered  with  plexiglass  and  supported 
by  Velcro  material  as  shown  in  Figures  5-6  and  5-7.  The  way  the  Plexiglass  was 
installed  in  the  coil  suction  and  discharge  is  shown  in  Figures  5-8  and  5-9.  Experiments 
were  performed  for  three  different  scenarios:  undampered  coils,  fully-dampered  coils, 
and  partially-dampered  coils  (at  the  coil  suction  only). 
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Figure  5-6.  Coil  suction  covered  with  plexiglass  and  supported  by  Velcro  material 


Figure  5-7.  Coil  discharge  covered  with  plexiglass  and  supported  by  Velcro  material 
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Figure  5-8.  Coil  suction  being  covered  by  plexiglass 


Figure  5-9.  Coil  discharge  being  covered  by  plexiglass 

The  initial  conditions  of  the  freezer  air  were  held  between  -14.4°C  and  -13.3°C 
dry-bulb  temperature  and  between  78%  and  80%  relative  humidity  for  all  tests  in  this 
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category.  The  coil-face  velocity  was  3.8  m/s.  The  length  of  the  experiments  was  6 hours 
under  refrigeration  and  15  minutes  under  defrosting.  After  each  experiment  was 
completed,  the  melt  from  the  defrosting  process  was  collected  in  order  to  determine  the 
defrost  efficiency. 

For  the  undampered  coil  scenario,  steady  state  conditions  of  freezer  air  at  the 
beginning  were  DBT=-  14°C  and  RH=  78%.  Figure  5-10  shows  an  image  of  the  coil 
before  starting  the  defrosting  process.  The  amount  of  melt  collected  after  completing  the 
defrosting  process  was  1.27  kg.  The  total  heat  input  for  the  defrosting  process  was  730  kJ 
and  the  total  heat  used  for  melting  the  frost  was  470.9  kJ,  thus  yielding  a defrost 
efficiency  of  64.5%. 


Figure  5-10.  Conditions  of  the  coil  at  the  beginning  of  the  defrosting  process 

For  the  fully-dampered  coil  scenario,  steady  state  conditions  of  freezer  air  at  the 
beginning  were  DBT=-  14.4°C  and  RH=  80%.  Again,  a visual  image  of  the  coil  before 
starting  the  defrosting  process  is  shown  in  Figures  5-1 1 and  5-12.  Figures  5-13  and  5-14 
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show  the  appearance  of  the  coil  eight  minutes  into  the  defrosting  process,  while  the 
appearance  of  the  coil  after  the  defrosting  process  has  been  completed  can  be  seen  in 
Figures  5-15  and  5-16.  The  amount  of  melt  collected  after  completing  the  defrosting 
process  was  1.81  kg,  while  the  total  heat  input  was  730  kJ.  The  total  heat  used  for 
melting  the  frost  was  671 .6  kJ,  thus  giving  a defrost  efficiency  of  92%. 

For  the  partially-dampered  coil  scenario,  steady  state  conditions  of  freezer  air  at 
the  beginning  were  DBT=-\4A°C  and  RH=  80%.  A visual  image  of  the  coil  before 
starting  the  defrosting  process  is  shown  in  Figure  5-17.  Figure  5-18  shows  the 
appearance  of  the  coil  eight  minutes  into  the  defrosting  process,  while  Figure  5-19  shows 
the  appearance  after  the  defrosting  process  has  been  completed.  The  amount  of  melt 
collected  after  completing  the  defrosting  process  was  1 .50  kg.  The  total  heat  input  during 
the  defrosting  process  was  729.1  kJ,  while  the  total  heat  used  for  melting  the  frost  was 
552.6  kJ,  thus  yielding  a defrost  efficiency  of  75.8%. 


Figure  5-11.  Conditions  of  the  coil  at  the  beginning  of  the  defrosting  process  (suction 
side) 
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Figure  5-12.  Conditions  of  the  coil  at  the  beginning  of  the  defrosting  process  (discharge 
side) 


Figure  5-13.  Conditions  of  the  coil  after  8 minutes  into  the  defrosting  process  (suction 
side) 
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Figure  5-14.  Conditions  of  the  coil  after  8 minutes  into  the  defrosting  process  (discharge 
side) 


Figure  5-15.  Conditions  of  the  coil  after  15  minutes  into  the  defrosting  process  (suction 
side) 
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Figure  5-16.  Conditions  of  the  coil  after  15  minutes  into  the  defrosting  process 
(discharge  side) 


Figure  5-17.  Conditions  of  the  coil  at  the  beginning  of  the  defrosting  process  (suction 
side) 
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Figure  5-18.  Conditions  of  the  coil  after  8 minutes  into  the  defrosting  process  (suction 
side) 


Figure  5-19.  Conditions  of  the  coil  after  15  minutes  into  the  defrosting  process  (suction 
side) 


Figure  5-20  presents  a graphical  display  of  the  defrost  efficiency  for  all  three 
scenarios.  As  can  be  seen,  the  defrost  efficiency  increases  when  a partially-dampered 
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coil  is  used  relative  to  an  undampered  coil.  Also,  the  highest  defrost  efficiency  is 
obtained  when  a fully-dampered  coil  is  used.  This  means  that  the  use  of  dampers  helps 
increase  the  amount  of  melt  collected  and  thus  the  defrost  efficiency.  Visual  observations 
made  inside  the  freezer  for  the  undampered  coil  indicated  that,  the  sublimated  frost  was 
visibly  diffusing  into  the  freezer  space,  resulting  in  a significant  alteration  of  the 
conditions  of  the  freezer  air.  When  dampered  coils  were  used,  on  the  other  hand,  the 
freezer  remained  clear  because  most  of  the  sublimated  frost  was  trapped  in  the  enclosure 
between  the  dampers  and  the  coil. 


Undampered  coil  Partiallydampered  coil  Fully  dampered  coil 


Figure  5-20.  Defrost  efficiency  for  the  three  scenarios  analyzed 

Visual  observations  made  inside  the  freezer  for  the  undampered  coil  indicated 
that,  the  sublimated  frost  was  visibly  diffusing  into  the  freezer  space,  resulting  in  a 
significant  alteration  of  the  conditions  of  the  freezer  air.  When  dampered  coils  were 
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used,  on  the  other  hand,  the  freezer  remained  clear  because  most  of  the  sublimated  frost 
was  trapped  in  the  enclosure  between  the  dampers  and  the  coil. 

Operating  Freezer  Coils  at  High  Face  Velocities 

In  this  section  the  effect  of  coil-face  velocity  is  investigated  for  a coil  operating 
under  supersaturated  conditions.  In  this  study  the  amount  of  melt  collected  during 
defrosting  and  the  defrost  efficiency  at  several  coil  face  velocities  were  determined.  This 
study  also  examines  the  differences  between  the  types  of  frost  that  forms  at  different  coil- 
face  velocities  during  the  operation  of  freezers  near  saturated  freezer  air  conditions. 
Some  of  the  quantities  analyzed  include  the  defrost  efficiency,  the  average  coil 
performance  load  rate  ( ACPLR ),  and  the  fan  defrost  contribution  ( FDC ).  All  these 
quantities  were  described  in  Chapter  4. 

Results  and  discussion.  In  order  to  analyze  the  effect  of  the  coil-face  velocity  for 
freezers  operating  under  supersaturated  conditions,  two  coil-face  velocities  were  tried; 
3.7  m/s  and  6.5  m/s.  Both  scenarios  were  performed  employing  the  same  freezer  initial 
conditions  for  five  hours  during  refrigeration  and  fifteen  minutes  during  defrosting. 
Initial  freezer  air  conditions  were  DBT=-  14.6°C  and  RH=  75%,  with  an  apparatus  dew- 
point of  -22°C.  During  refrigeration,  0.38  kg/h  of  steam  were  injected  into  the  freezer 
space.  Supersaturated  air  existed  throughout  the  refrigeration  portion  of  the  R/D  cycle. 

For  the  first  scenario  (case  a),  the  coil-face  velocity  was  3.7  m/s.  The  freezer  air 
temperature  stayed  between  -12°C  and  -14°C  for  a period  of  four  hours  during  the 
refrigeration  portion  of  the  cycle.  Shortly  afterward,  air  temperature  in  the  freezer  space 
started  to  rise  due  to  frost  deposition  on  the  coil  and  the  subsequent  reduction  in  cooling 
capacity  (complete  cycle  is  graphically  represented  in  Figure  5-21). 
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Figure  5-21.  Freezer  and  coil  temperatures  during  both  refrigeration  and  defrosting 
(Case  a) 

For  the  second  scenario  (case  b),  the  coil-face  velocity  was  6.5  m/s.  Again, 
freezer  air  was  supersaturated  during  the  five  hours  of  refrigeration.  The  complete  R/D 
cycle  for  this  case  is  shown  in  Figure  5-22.  As  can  be  observed,  freezer  air  temperature 
stayed  constant  for  about  3.8  hours,  shortly  after  which  it  started  to  rise  for  the  same 
reasons  explained  in  the  first  case. 

Successive  images  of  the  coil  during  the  defrosting  portion  of  the  R/D  cycle  for 
both  scenarios  were  taken  and  are  presented  in  Figures  5-19  through  5-22.  The  effect  of 
the  coil-face  velocity  can  be  understood  by  comparing  the  corresponding  images  of  the 
two  cases.  The  initial  frost  thickness  in  Case  a was  5 mm,  while  that  for  Case  b was  4.5 
mm.  As  can  be  seen  in  Figure  5-19,  the  frost  in  Case  a looks  more  compact  and  dense 


than  that  of  Case  b. 
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Figure  5-22.  Freezer  and  coil  temperatures  during  both  refrigeration  and  defrosting 
(Case  b) 

Five  minutes  into  defrosting,  Figure  5-24  indicates  that  in  Case  a the  frost  at  the 
top  of  the  coil  falls  down  instead  of  melting.  This  is  contrary  to  what  is  observed  for 
Case  b,  where  the  frost  at  the  top  of  the  coil  was  observed  to  melt  instead  of  falling  in 
chunks.  This  may  say  something  about  the  degree  of  cohesiveness  of  the  crystal  structure 
of  the  two  cases  and  the  associated  effect  of  a higher  or  lower  coil-face  velocity  on  the 
formation  pattern. 

Ten  minutes  into  defrosting,  Figure  5-25  indicates  that  in  Case  a the  frost  at  the 
top  of  the  coil  continues  to  fall  down  instead  of  melting.  For  Case  b,  the  frost  at  the  top 
of  the  coil  continues  to  melt,  with  only  a small  amount  of  frost  falling  in  chunks.  It  can 
also  be  observed  that  the  frost  associated  with  Case  a is  thicker  than  that  of  Case  b. 

Fifteen  minutes  into  defrosting,  Figure  5-26  shows  that  in  Case  a the  frost  at  the 
top  of  the  coil  continues  to  fall  down  in  chunks  instead  of  melting,  with  a significant 
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amount  of  frost  still  attached  to  the  coil  surface.  At  the  higher  coil-face  velocity  (Case  b), 
the  frost  at  the  top  of  the  coil  continues  to  melt,  with  little  frost  falling  down  in  chunks. 
Although  in  Case  b the  coil  shows  more  frost  than  Case  a,  the  frost  layer  in  the  latter  case 
is  thinner  than  that  at  the  lower  coil-face  velocity  (Case  a).  This  is  because  a significant 
amount  of  the  removed  frost  at  the  higher  coil-face  velocity  has  already  melted  instead  of 
falling  in  chunks.  Melted  frost  is  drained  outside  the  freezer,  whereas  frost  falling  in 
chunks  eventually  evaporates  or  sublimates  into  the  freezer  space  and  is  recycled  back 
into  frost  on  coils  operating  in  the  vicinity  of  the  coil  being  defrosted,  thus  increasing  the 
frequency  of  defrosting  these  coils. 

Results  of  this  experiment  are  rather  interesting  as,  in  both  cases;  the  fifteen 
minutes  of  defrosting  were  not  sufficiently  long  to  completely  defrost  the  coil.  The 
amount  of  melt  collected  at  the  conclusion  of  defrosting  in  Case  a was  0.36  kg.  The  total 
heat  input  in  the  defrosting  process  was  649.5  kJ  and  the  total  heat  used  in  the  melting 
process  was  141.7  kJ,  thus  yielding  a defrost  efficiency  of  21.8%.  For  Case  b,  the 
amount  of  melt  collected  was  1.18  kg.  The  total  heat  input  in  the  defrosting  process  was 
661.6  kJ  and  the  total  heat  used  in  the  melting  process  was  465.3  kJ,  thus  yielding  a 
defrost  efficiency  of  70.4%.  Obviously,  the  defrost  efficiency  for  Case  b is  higher  than 
that  for  Case  a.  This  necessarily  implies  that  frost  depositing  at  higher  coil-face 
velocities  is  more  likely  to  collect  in  the  melt  pan  than  to  sublimate  into  the  freezer  space. 
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B 

Figure  5-23.  Comparison  of  the  coil  condition  after  5 hours  of  refrigeration  (same  as  the 
beginning  of  defrosting)  for  different  coil-face  velocities:  A)  Case  a:  3.7 
m/s  and  B)  Case  b:  6.5  m/s. 
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B 

Figure  5-24.  Comparison  of  coil  condition  5 minutes  into  defrosting  for  different  coil- 
face  velocities:  A)  Case  a:  3.7  m/s  and  B)  Case  b:  6.5  m/s. 
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A 


Figure  5-25.  Comparison  of  coil  condition  10  minutes  into  defrosting  for  different  coil- 
face  velocities:  A)  Case  a:  3.7  m/s  and  b)  Case  b:  6.5  m/s. 
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A 


B 


Figure  5-26.  Comparison  of  coil  condition  15  minutes  into  defrosting  for  different  coil- 
face  velocities:  (a)  Case  a:  3.7  m/s  and  (b)  Case  b:  6.5  m/s 
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In  order  to  be  able  to  create  a meaningful  plot  of  the  defrost  efficiency  as  a 
function  of  coil-face  velocity,  two  more  experiments  were  conducted.  The  coil-face 
velocities  used  in  these  experiments  were  1.6  m/s  and  5.23  m/s.  Both  experiments  were 
performed  under  the  same  initial  freezer  conditions  and  the  same  rate  of  steam  injection 
as  the  previous  experiments.  Results  for  these  two  coil-face  velocities  along  with  the 
previous  ones  are  shown  in  Table  5-1.  The  corresponding  graphical  representation  is 
shown  in  Figure  5-27. 


Table  5-1.  Defrost  efficiency  and  melt  collected  in  the  defrosting  process  for  different 
coil-face  velocities 


Coil-face 

velocity 

(m/s) 

Heat  input  into 
defrosting 

(U) 

Heat  use  in  the 
melting  process 

qj) 

Melt 

collected 

(kg) 

Defrost 

efficiency 

(%) 

1.60 

642.0 

115.6 

0.27 

18.0 

3.69 

649.5 

141.7 

0.36 

21.8 

5.23 

651.0 

331.0 

0.84 

50.8 

6.50 

661.6 

465.3 

1.18 

70.4 

As  can  be  seen  in  Figure  5-27,  for  coil-face  velocities  below  3.5  m/s,  the  defrost 
efficiency  is  almost  constant.  However,  at  higher  coil-face  velocities,  the  defrost 
efficiency  can  be  seen  to  increase  linearly  with  the  velocity.  This  may  be  explained  by 
attributing  the  higher  amount  of  melted  frost  in  the  higher  velocity  case  to  a more 
coherent  crystal  structure  that  prevents  frost  falling  in  chunks  in  the  higher  velocity  case. 
It  may  also  cause  frost  to  be  more  physically  attached  to  the  coil  surface,  thus  making  it 
more  susceptible  to  the  hot  gas  heat  emanating  from  the  inside  of  the  coil. 
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Figure  5-27.  Defrost  efficiency  vs.  coil-face  velocity 

It  is  important  to  note  here,  however,  the  potential  implications  of  operating 
freezer  coils  at  higher  face  velocities  vis-a-vis  the  overall  cost  of  the  complete  R/D  cycle. 
While  we  may  have  a significantly  larger  defrost  efficiency  at  higher  face  velocities,  the 
additional  fan  power  required  at  higher  velocities  may  negate  the  advantages  associated 
with  the  higher  defrost  efficiency.  This  may  be  particularly  true  since  the  duration  of  the 
refrigeration  portion  of  the  R/D  cycle  is  typically  several  folds  longer  than  that  of  the 
defrosting  portion.  For  example,  in  the  experiments  performed  here,  the  fan  horsepower 
for  the  3.7  m/s  case  was  1.12  kW  smaller  that  of  the  6.5  m/s.  Since  the  refrigeration 
portion  of  the  cycle  lasted  five  hours,  the  extra  energy  required  for  the  6.5  m/s  case  is  5.6 


kWh. 
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On  the  other  hand,  operating  at  lower  face  velocities  will  result  in  a higher 
defrosting  frequency  because  of  the  associated  larger  amount  of  moisture  transfer  into  the 
freezer  space.  Thus,  in  a given  24-hour  period,  a coil  operating  at  a lower  face  velocity 
will  be  defrosted  several  times  more  than  one  operating  at  a higher  velocity.  This  is  an 
additional  source  of  energy  expenditure  that  should  be  accounted  for  in  deciding  which 
face  velocity  to  operated  the  coil  at.  Furthermore,  defrosting  coils  employing  the  hot-gas 
refrigerant  method  frequently  requires  artificial  generation  of  hot  gas  in  order  to  have  a 
sufficient  amount  of  it  to  successfully  perform  the  defrosting  process.  This  hot-gas 
generation  process  may,  for  example,  require  the  use  of  electric  energy  in  the  form  of 
band  heaters  wrapped  around  the  liquid-vapor  separator  in  a liquid  overfeed  refrigeration 
system.  A case-in-point  is  the  system  employed  in  this  experimental  investigation, 
whereby  up  to  10  kW  of  electric  band  heater  power  could  be  used  to  boil  the  refrigerant 
liquid  in  the  liquid-vapor  separator  for  the  duration  of  at  least  part  of  the  defrosting 
process.  Obviously,  more  frequent  defrosting  cycles  (as  in  the  lower  face  velocity  case) 
would  make  this  option  less  attractive  than  the  higher  velocity  scenario.  This  naturally 
needs  to  be  factored  into  the  decision  making  process. 

Table  5-2  gives  numerical  values  of  the  fan-defrost  contribution,  while  the 
average  coil  performance  load  rate  values  are  provided  in  Table  5-3.  The  values  given  in 
these  tables  are  for  coil-face  velocities  of  1.6  m/s,  3.7  m/s  and  6.5  m/s. 

Tabulated  values  of  the  fan-defrost  contribution  ( FDC)  are  significantly  larger 
than  the  values  provided  by  ASHRAE  (1998).  ASHRAE  suggests  that  the  theoretical 
contribution  of  the  fan  and  defrost  heat  loads  to  the  net  refrigeration  load  is 
approximately  15%  or  more.  This  suggests  that  a larger  contribution  should  be 
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considered  in  cases  involving  a significant  latent  load.  Also,  the  average  coil 
performance  load  rate  for  all  three  coil-face  velocities  is  remarkably  similar.  This  latter 
result  gives  confidence  in  the  accuracy  of  the  experimental  data  collected. 


Table  5-2.  Fan-defrost  contribution  for  coil-face  velocities  of  1.6  m/s,  3.7  m/s,  and  6.5 
m/s 


Coil-face  velocity 

1.6  m/s 

3.7  m/s 

6.5  m/s 

GLr,  (kJ) 

39703.1 

40388.9 

41654.9 

QtotaU  (kJ) 

5719.4 

5719.4 

5719.4 

Q.otal  (W) 

45422.6 

46108.3 

47374.3 

Qfan  OH) 

9495.0 

10180.8 

11446.8 

Qne,  (W) 

35927.6 

35927.6 

35927.6 

(kJ) 

642.0 

649.5 

661.6 

Qmell  (kJ) 

115.6 

141.7 

465.3 

Qdefjoad  (kJ) 

526.4 

507.8 

196.3 

FDC  (%) 

27.9 

29.8 

32.4 

Table  5-3.  Average  coil  performance  load  rate  for  coil-face  velocities  of  1.6  m/s,  3.7  m/s, 

and  6.5  m/s 

Coil-face  velocity 

1.6  m/s 

3.7  m/s 

6.5  m/s 

V (s) 

18000 

18000 

18000 

V (s) 

900 

900 

900 

* total  00 

18900 

18900 

18900 

Qtota,  (kJ) 

45422.6 

46108.3 

47374.3 

QtotaiA  kW) 

2.40 

2.44 

2.51 

Qdefjoad  (kJ) 

526.4 

507.8 

196.3 

Qdef  .load  ,av  (kW) 

0.03 

0.03 

0.03 

ACPLR  (kW) 

2.43 

2.47 

2.52 

Conclusions 

In  this  chapter  three  alternative  energy  conservation  technologies  in  industrial 

freezers  and  refrigerated  facilities 

have  been  presented. 

The  use  of  heat  pipes  is 

advantageous  for  freezer  operation,  because  the  energy  required  for  the  cooling  and 
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dehumidification  process  decreases.  This  reduction  in  energy  can  be  translated  into 
financial  savings.  The  use  of  partially-dampered  coils  was  shown  to  increase  the  defrost 
efficiency  by  18%,  while  the  use  of  fully-dampered  coils  increases  the  efficiency  by  43%. 
These  dampers  prevent  evaporated  water  and  sublimated  frost  from  going  back  into  the 
freezer  space,  thus  helping  reduce  the  defrosting  frequency  of  other  coils.  Frost 
deposited  at  higher  coil-face  velocities  is  more  cohesive  and  thus  more  likely  to  collect  in 
the  melt  pan  during  defrosting  than  to  sublimate  back  into  the  freezer  space.  At  lower 
coil  face  velocities,  the  relative  lack  of  cohesiveness  of  the  frost  is  likely  to  result  in  a 
significant  amount  of  moisture  transfer  to  the  freezer  space  and/or  a structural  collapse  of 
chunks  of  frost  on  the  freezer  floor.  This  necessarily  results  in  a smaller  amount  of 
collected  melt  and  a correspondingly  smaller  defrost  efficiency.  It  appears  that  a 
significant  departure  to  a cohesive  enough  layer  (that  significantly  affects  the  defrost 
efficiency)  occurs  at  a coil  face  velocity  of  3.5m/s.  Frost  deposited  at  face  velocities  of 
6.5  m/s  results  in  more  than  a four-fold  increase  in  the  defrost  efficiency  relative  to  that 
deposited  at  1.6  m/s.  Tabulated  values  of  the  fan-defrost  contribution  ( FDC)  also  need  to 
take  into  account  the  presence  of  latent  loads  in  the  freezer  space.  This  represents  a 
marked  departure  from  the  previous  ASHRAE  guidelines  regarding  the  so-called  fan- 
defrost  contribution.  Furthermore,  the  average  coil  performance  load  rate  for  all  three 
coil-face  velocities  is  remarkably  similar.  This  latter  result  gives  confidence  in  the 
accuracy  of  the  experimental  data  collected. 


CHAPTER  6 

FORMULATIONS  OF  THERMODYNAMIC  PROPERTIES  OF  SUPERSATURATED 

AIR 

Objectives 

The  most  exact  formulations  of  thermodynamic  properties  of  moist  air  in  the 
temperature  range  of  -100°C  to  200°C  are  based  on  the  study  of  Hyland  and  Wexler 
(1983a,  b).  More  recent  studies  by  Sauer  et  al.  (2001)  and  Nelson  et  al.  (2001)  provide 
psychrometric  data  for  moist  air  in  the  temperature  range  of  200°C  to  320°C  and 
humidity  ratios  from  0 to  1 kg/kga  at  pressures  0.07706,  0.101325,  0.2,  1.0,  and  5 MPa. 
Both  studies  developed  the  psychrometric  data  using  the  most  current  values  of  the  virial 
coefficients,  enthalpy,  and  entropy  of  both  air  and  water  vapor.  Other  psychrometric  data 
were  generated  by  Stewart  et  al.  (1983),  who  created  psychrometric  charts  in  SI  units  at 
low  pressures.  More  psychrometric  charts  are  available  in  the  ASHRAE  brochure  on 
psychrometry  (Olivieri  et  al.  1996). 

While  psychrometric  tables  and  charts  exist  for  subsaturated  moist  air  for  a 
relatively  large  temperature  range,  no  such  data  is  available  under  supersaturated 
conditions.  This  deficiency  in  the  body  of  data  has  historically  impeded  our  ability  to 
analyze  and  study  frost  and  ice  formation  processes  on  coils  exposed  to  supersaturated 
air.  These  coils  may  exist  in  a variety  of  engineering  situations  such  as  industrial  walk- 
in  freezers  and  loading-docks,  to  name  a few.  While  every  attempt  should  be  made  to 
avoid  operating  coils  under  those  conditions,  developing  this  body  of  data  should  enable 
the  refrigeration  engineer  and  system  designer  to  better  predict  the  performance  of  the 
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system  as  well  as  design  a more  optimal  system.  The  objective  of  this  section  is  thus  to 
provide  a systematic  way  of  calculating  key  air  properties  under  supersaturated 
conditions  for  the  purpose  of  improving  our  prediction  capability  under  those  conditions. 


In  this  section,  the  development  of  psychrometric  formulations  as  well  as 
thermodynamic  properties  for  supersaturated  air  from  -40  to  40°C  will  be  presented. 
Supersaturated  air  is  a mixture  of  saturated  air  and  suspended  liquid  water  droplets  or  ice 
crystals  depending  on  whether  the  prevailing  air  temperature  is  larger  or  smaller  than  the 
freezing  point  temperature,  respectively.  Therefore  the  properties  of  supersaturated  air 
may  be  determined  using  the  individual  properties  of  those  of  saturated  air  and  those  of 
the  suspended  liquid  droplets  or  ice  crystals.  These  properties  are  developed  employing 
the  virial  equation  of  state  and  other  relevant  methods.  The  virial  equation  can  be  written 
as 


In  Equation  6-1,  B is  the  second  virial  coefficient  for  gas  “z”,  r is  the  third 
virial  coefficient  for  gas  “z”,  vm  is  the  molar  specific  volume,  and  R is  the  universal  gas 
constant  (8.31441  kPa  m3/(kmol  K)).  The  subscripts  “z”  are  equal  to  “a”  for  dry  air  and 
“w”  for  water  vapor.  The  subscript  “m  ” will  replace  the  z’s  subscripts  to  indicate  the  air- 
water  vapor  mixture  (Olivieri  et  al.  1996) 


Properties  of  Supersaturated  Moist  Air 


(6-1) 


(6-2) 


(6-3) 
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where  xa  and  xw  are  the  mole  fractions  of  air  and  water  vapor  respectively,  B and  B 
are  the  second  virial  coefficients  for  dry  air  and  water  vapor,  respectively,  C and  Cuuw 
are  the  third  virial  coefficients  for  dry  air  and  water  vapor,  respectively,  and  B , Caaw, 
and  Caww  are  the  cross  viral  coefficients  for  moist  air. 

The  second  and  third  virial  coefficients  for  dry  air,  B (cm2/mol)  and  r 

aa  ' ^ aaa 

(cm6/mol2),  respectively,  are  given  by  Olivieri  et  al.  (1996)  according  to  the  following: 

D 6687.72  0.210141xl07  0.924746xl08  (6-4) 

5 =34.9568 = + 


^ 0.1 90905x1 06  0.632467xl08 

=1259.75 + 


(6-5) 


The  second  and  third  virial  coefficients  for  water  vapor,  B (cm2/mol)  and  Cw 


6 2 

(cm  /mol  ),  respectively,  are  expressed  as  (Olivieri  et  al.  1996) 


B =7?T 


0.70x1 0”8  - 0. 1471 84x1 0”8  exp 


1734.29 


7 


(6-6) 


Cwww=(RTy  |o.104x10'14  - 0.335297x1 0~17  exp 


f 3645.09' 

V T j 


+ 


0.70x1  O'8  - 0.1 47 1 84x1  O’8  exp 


f 1734.29  V 

l T [ 

/ 

(6-7) 


The  cross  virial  coefficients  for  moist  air,  B (cm2/mol),  Caaw  (cm6/mol2),  and 


Caww  (cm6/mol2),  are  expressed  as  (Olivieri  et  al.  1996) 


99 


5 =32.366097- 


0.141 138x1 05 


482.737  + 


0.105678x1 06 


0. 319327x10" 


0.1244535xl07  0.2348789x10'° 

rji2  iji3 

0.656394x1 08  0.294442x10" 

? 1 i 

J»2  rpi 


(6-8) 

(6-9) 


exp 


-10.728876x10  + 


2 . 0.347802xl04  0.383383xl06 


(6-10) 


0.33406xl08 


The  above  coefficients  have  been  calculated  using  Equations  6-4  through  6-10 
and  are  summarized  in  Table  6-1.  The  temperature  derivatives  have  also  been  calculated 
and  are  given  in  Table  6-2. 


Table  6-1.  Value  of  the  virial  coefficients 


t 

(°C) 

Baa 

(cm3/mol) 

r 

aaa 

(cm6/mol2) 

Bww 

(cm3/mol) 

c 

WWW 

(cm6/mol2) 

Baw 

(cm3/mol) 

r 

aaw 

(cm6/mol2) 

r 

aww 

(cm6/mol2) 

-40 

-25.0889 

1604.4441 

-4837.4598 

-54249456.48 

-51.8588 

971.079 

-795863.7234 

-30 

-21.6586 

1544.3847 

-3711.5848 

-30621660.01 

-47.4018 

941.8068 

-556734.3014 

-20 

-18.5521 

1492.552 

-2911.788 

-18137967.62 

-43.3786 

913.3703 

-401252.0739 

-10 

-15.7288 

1447.6264 

-2329.5237 

-11211706.74 

-39.7299 

886.3439 

-296660.4159 

0 

-13.1543 

1408.5352 

-1896.2673 

-7198239.87 

-36.4067 

861.0182 

-224216.3255 

10 

-10.7993 

1374.4001 

-1567.5254 

-4780658.789 

-33.3679 

837.5067 

-172750.726 

20 

-8.6387 

1344.497 

-1313.6677 

-3272857.046 

-30.5792 

815.815 

-135368.5722 

30 

-6.6509 

1318.225 

-1114.5057 

-2302551.893 

-28.0114 

795.8843 

-107679.4618 

40 

-4.8173 

1295.0824 

-955.9962 

-1660201.825 

-25.6398 

777.6198 

-86810.3558 
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Table  6-2.  Temperature  derivatives  of  the  virial  coefficients 


t 

(°C) 

dBjdT 

(cm3/mol  K) 

dCaJdT 

(cm6/mol2  K) 

dB^/dT 

(cm3/mol  K) 

dCwm//dT 
(cm6/mol2  K) 

dBjdT 
(cm3/mol  K) 

dCaaJdT 

(cm6/mol2  K) 

dCamv/dT 
(cm6/mol2  K) 

-40 

0.3608 

-6.4687 

134.0212 

4741841.882 

0.4697 

-2.9396 

29764.0056 

-30 

0.3261 

-5.5702 

94.0267 

2485668.989 

0.4229 

-2.8981 

19018.5847 

-20 

0.2959 

-4.8182 

67.6965 

1370887.445 

0.3827 

-2.7799 

12603.6273 

-10 

0.2694 

-4.1847 

49.8733 

790868.1648 

0.3479 

-2.6207 

8617.1634 

0 

0.246 

-3.6481 

37.5051 

474898.199 

0.3175 

-2.4427 

6052.6767 

10 

0.2254 

-3.1909 

28.7286 

295559.8373 

0.2908 

-2.2596 

4352.5921 

20 

0.2071 

-2.7996 

22.3742 

189952.7577 

0.2674 

-2.0798 

3195.4595 

30 

0.1908 

-2.4631 

17.6890 

125666.3880 

0.2466 

-1.908 

2389.3459 

40 

0.1762 

-2.1724 

14.1768 

85341.8427 

0.2281 

-1.7468 

1816.0418 

Thermodynamic  Properties 

Specific  volume:  The  following  equations  represent  the  specific  volume  for  dry,  moist, 
and  supersaturated  air,  and  for  saturated  ice  and  water: 


■ Dry  air:  The  specific  volume  for  dry  air  in  m3/kmol  can  be  determined  from 

Equation  6-1  by  changing  the  subscripts  by  “a”  and  vm  by  va 


v„  =■ 


RTf 


K 

v„ 


+ - 


C„ 


(6-11) 


Moist  subsaturated  air:  The  specific  volume  for  moist  air  in  m3/kmol  can  be 
determined  from  Equation  6-1  by  changing  the  subscripts  “z”  by  “m” 

(6-12) 


v„  = 


RT(Bm  Cm ^ 

m i ttl 


J 


The  specific  volume  can  be  expressed  per  unit  mass  of  dry  air  in  m3/kg  as 


v„  =■ 


28.9645  x„ 


(6-13) 


At  saturation  vm  becomes  vs , and  xa  and  xw  become  xas  and  xws , respectively. 
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■ Saturated  ice:  The  specific  volume  of  saturated  ice  in  m3/kg  (173.1 5K  < T < 273.15) 
is  obtained  from  (Olivieri  et  al.  1996) 

vice  =0.1070003xl0'2-0.249936xl0-7r  + 0.37161U10'9r2  (6-14) 


Saturated  water:  The  specific  volume  of  saturated  water  in  m3/lcp  is 

obtained  from 

(Olivieri  et  al.  1996) 

a6  +aj 

(6-15) 

y liquid 

±v. 

i-0 

where  ao  = -2403.360201 

a4  = 0.373497936*1  O'6 

a/ = -1.40758895 

as  = -2. 1203787xl0'10 

a2  = 0.1068287657 

a6  = -3.424442728 

a2  = -0.291449235  IjcIO"3 

a7  = 0.01619785 

■ Moist  supersaturated  air:  In  order  to  determine  the  specific  volume  of  supersaturated 
air,  it  is  necessary  first  to  introduce  the  definition  of  the  humidity  ratio  of 
supersaturated  air  (Wsuper).  This  can  be  thought  of  as  the  humidity  ratio  of  air  at 

saturation  (ifjplus  the  mass  of  excess  moisture  beyond  saturation  (either  liquid  or 
ice)  per  unit  mass  of  dry  air  ( mdrop /ma).  This  can  be  expressed  by 


W + mdrop 

v sup  er 

m„ 


- w | mdr°P 

’ S 


(6-16) 


The  specific  volume  of  supersaturated  air  can  be  expressed  per  unit  mass  of  mixture 
as  follows: 


maVg  +mdropVdroP 


super 


m 


super 


(6-17) 
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where  msuper  is  the  total  mass  of  the  mixture  and  it  can  be  expressed  as 


m,uPer  = ™a  + + mdrop 


(6-18) 


where  ma  is  the  mass  of  dry  air,  mws  is  the  mass  of  saturated  water  vapor  , and 

is  the  mass  of  excess  moisture  beyond  saturation. 

Introducing  Equation  6-18  into  Equation  6-17  gives 


m 


sup  er 


-v„  +- 


m„ 


m 


ma  + mws  + mdrop  ma  + 


vra+- 


drop 


ma+mws+mdrop 


v drop 


drop 


(6-19) 


Equation  6-19  can  be  written  as: 


sup  er 


i +ws  + tr^E- 

m„ 


-v„  +- 


W. 


v v, 

ws  ^ drop 


i+ws  + rn^  1 +ws  + rn^ 

m„  s m. 


(6-20) 


Using  the  relation  for  the  supersaturated  humidity  ratio  (Equation  6-16)  Equation 
6-20  can  be  expressed  as  follows: 


sup  er 


= (V*  + V“rop  kup,  - K )] 

1 ' YY  sup  er 


(6-21) 


The  specific  volume  of  supersaturated  air  can  be  also  expressed  per  unit  mass  of  dry 
air  as  follows: 


m, 


v = v 

super  super 


sup  er 

m„ 


= v£ 


super 


(i+O 


(6-22) 


Using  Equations  6-30  and  6-31  the  specific  volume  of  supersaturated  air  can  be 
expressed  as 


v = v 

sup  er  j 


+ V. 


drop 


k 


super 


(6-23) 
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where  vdrop  can  either  be  vliquid  or  vice  depending  on  the  air  temperature. 


Specific  enthalpy:  The  following  virial  equations  represent  the  specific  enthalpy  for  dry 
and  moist  air,  and  for  saturated  ice  and  water. 

■ Dry  air:  The  specific  enthalpy  of  dry  air  can  be  expressed  in  kJ/kmol  by 


l.='Zb,T‘+ST 

;=o 


n dB 

\ . 

f dC  ^ 

Baa  ~T — 

v + 

C 0.5T 

v~2 

[{  dT 

i “ 

Trj-f 

V dT  ) 

a 

(6-24) 


where  bo  = -7907.8691 
b/  = 28.709015 
&2  = 0.002643 1805 


65  = -0.1 0405863x1  O'4 
b4  = 0.1866041xl0'7 
b5  = -9.7843331  xl0‘12 


■ Moist  subsaturated  air:  The  specific  enthalpy  of  moist  subsaturated  air  can  be 

expressed  in  kJ/kmol  by 


( 5 


\i=0 


+ RT 


X eiT  + *.,0  + X fiT‘  + C.0 

V-=0 


X 

) 


B-T 


dB  V 

m 1 

dT 


v‘  + 


( dC  ^ 

C -0.51-—^- 

v dT  j 


vj 


(6-25) 


where 


e0  = -6.3290874 
ej  = 28.709015 
e2  = 0.0026431805 
e5  = -0.000010405863 
e4  = 0.00000001866041 
e5  = -9.7843331  xlO'12 
ia.o=-  7.9141982  kJ/kmol 


fo  = -0.005008 
f,=  32.491829 
/2  = 0.0065576345 
f3  = -0.000026442147 
f4  = 0.000000051751789 
f5  = - 3.1541624  jcIO11 
iw,o  = 35.99417  kJ/kmol 


The  enthalpy  of  moist  air  per  unit  mass  kJ/kg  is  given  by 
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i = 


28.9645  x 


(6-26) 


At  saturation  im  becomes  is,  andxa  ,xw,  and  vm  become  xas,xws,  and  vf, 
respectively. 

■ Saturated  ice:  The  enthalpy  of  saturated  ice  water  in  kJ/kg  is  given  by 


(6-27) 


where  g0  = -647.595  g3  = 0.1 083437x1  O'5 

g,  = 0.274292  g4  = 0.107  xlO'5 

g2  = 0.002910583 

■ Saturated  water:  The  enthalpy  of  saturated  liquid  water  in  (kJ/kg)  is  given  by 


i liquid  a 0.01214  + TV liquid 


dp^ 

dT 


jcIO 


-3 


(6-28) 


where  a is  given  by  Equation  6-29  for  the  range  of  273.15  K <T<  373.15  K (Olivieri 
et  al.,  1996) 

a = -0.563473xl0-°036(r-27315)  (6'29) 

i=0 

where  j0  = -0.1 141 138xl04  j3  = 0.1451 133xl0'6 

j,  = 4.1930463  j4  = -0.100523xl0‘9 

j2  = -0.8134865X10"4 

In  order  to  determine  the  enthalpy  of  saturated  water,  the  derivative  of  the  saturation 
pressure  with  temperature  is  needed.  The  saturation  pressure  over  ice  for  the 
temperature  range  of  173. 15K  < T<  273. 15K  is  given  by  (ASHRAE  2001) 
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log  ,(P,)  = J + h +l1T+IsT2+l,r+lsT,+lsLog,(T)  (6'30) 

where:  l0  = -0.567453  59  x 104  l4=  0.20747825  x 10'8 

/y  - -5.152305  jc  101  l5  = -0.9484024  x 1012 

12  = -0.9677843  x 1 O'2  l6=  0.41635019x10* 

1 3 = -0.62215701  x 10'6 

The  saturation  pressure  over  liquid  water  for  the  temperature  range  of  273. 15K  < T < 
473. 15K  can  be  written  as  (ASHRAE  2001) 


logf  (ps)  = Y + n \+n2T  + n3T2  + "J*  + nsLoSe  (T) 


(6-31) 


where  n0  = -0.58002206  x 104  n3  = 0.41764768  x 10'4 

n,  = -0.5516256x  10*  n4  = -0.14452093  x 10‘7 

n2  = -0.48640239  x 10*  n5  = 0.65459673  x 10* 

In  both  Equations  6-30  and  6-31,  Ps  is  in  Pa. 

■ Moist  supersaturated  air:  The  enthalpy  of  supersaturated  air  can  be  expressed  per  unit 
mass  of  mixture  as  follows: 


...  maia+m¥aiw  + mdropidr0p 


sup  er 


m. 


sup  er 


(6-32) 


where  msuper  is  the  total  mass  of  the  mixture  (Equation  6-18) 
With  the  use  of  Equation  6-18  Equation  6-32  can  be  written  as 


m„ 


+- 


m.. 


m 


drop 


(6-33) 


"super  *a  ' , , *w  ' , 1 drop 

ma+mw+mdroP  ma+mws+mdrop  ma+mws+  mdrop 


or 
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^sup  er 


Ir.  +' 


W. 


' H"  - 


YYl  YYl 

1 + W.+-4&-  1 + W+ 


mdroP/ma 
TTL 

l + WM+-±^ 
m„ 


(6-34) 


tdrop 


Using  the  relation  for  the  supersaturated  humidity  ratio  (Equation  6-16) 
Equation  6-34  can  be  expressed  as  follows: 


super 


1 + JV 


[is+idroP(W^er-Ws\ 


sup  er 


(6-35) 


The  enthalpy  of  supersaturated  air  can  be  also  expressed  per  unit  mass  of  dry  air  as 
follows: 


m 


hup  er  hup  er 


sup  er 

m„ 


= T (l  + W ) 

Lsup  er  \A  r sup  er  ) 


(6-36) 


Using  Equations  6-35  and  6-36  the  enthalpy  of  supersaturated  air  per  unit  mass  of  dry 
air  can  be  expressed  as: 

i =i+i,  (w  -w)  (6-37) 

sup  er  s drop  \ sup  er  s / K 7 


where  idrop  can  either  be  iliquid  or  iice  depending  on  the  air  temperature. 

Specific  entropy:  The  following  equations  represent  the  specific  entropy  for  dry  and 
moist  air: 

■ Dry  air:  The  specific  entropy  of  dry  air  in  kJ/kmol  K can  be  determined  from  Olivieri 
et  al.  (1996) 


sa  =YjoiT‘+o5  log eT  + R 


i=0 


V 

f p ] 

f pva  Yl 

-log. 

+ log. 

„v 

V 1 atm  J 

l K1  JJ 

(6-38) 


B..  +- 


dB. 


1 1 

■ + — 


dT  )va  2 


C +T- 


dCaaa 

dT 


v2 
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where  o0  = -0.16175159xl03  03=0.24880547x1  O'7 

o ,=0. 52863609x1  O'2  o4=- 0. 1 22304 1 6x  1 O'10 

O2=-0. 1 5608795x1 0‘4  o5= 0.28709 1 0 1 5x1 02 

■ Moist  subsaturated  air:  The  specific  entropy  of  moist  air  in  kJ/kmol  K can  be 
determined  from 


sm  = x„ 


( 4 

YjU‘T‘  +U5  1o§.  T + Sa,0 

V<= 0 
|Y 


( 5 


+ x„ 


YjWiT‘  +VV5  l0ge  T + S 

V i=o 


w,0 


+ 


(6-39) 


R< 


-log, 


r p \ 


W 

r 


B_  +■ 


p 

V 1 atm  J 

dB_ 


+ xa  log, 


f r\ — \ 

tX, 

\XaRT  J 


1 ( 


dT 


■ + 


C_  +T- 


dCm 

dT 


log, 


2v : 


Pv. 


\\ 


\x*rtjj 


where  u0  = -0.34373874x1 02 
uj= 0.52863609x1  O'2 
w?=-0. 1 5608795x1  O'4 
W3=0. 24880547x1  O'7 
w^=-0. 1223041 6xlO'10 
«5=0.287090 1 5x1  02 
sa  o = -0.196125465  kJ/kmol  K 
sw,0  = -0.06331449  kJ/kmol  K 


w0  = 0.21966034x10' 
W;=0. 1 974381 9x10"' 
W2=-0.70128225x10‘4 
w5=0. 14866252x1  O'6 
w4=-0.14524437x10'9 
W5=0. 55663583x10'° 
W6=0. 32284652x1 02 


The  specific  entropy  per  unit  mass  of  dry  can  be  expressed  air  in  kJ/kg  K as 


s_  = ■ 


(6-40) 


28.9645  x„ 


At  saturation  sm  becomes^,  andxfl  and  xw  become  xas and  xWJ , respectively. 
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■ Saturated  ice:  The  specific  volume  of  saturated  ice  in  kJ/kg  K in  the  range  (173.1 5K  < 
T < 273.15)  is  obtained  from 

sice  =-4.470727 + 0.5021 09x1 0-2r  + 0.56241 55x1 0~5r2  - 0.274292  loge  T (6-41) 
+ (-  0.249936x1  O’7  + 0.743288x1 0"9r)(l  01 325  - Pm  )xl0-3 

■ Saturated  water:  The  specific  volume  of  saturated  water  in  kJ/kg  K in  the  range 
(273. 15K  <T<  373.15)  is  obtained  from 

snquid  = E ziT  + z4  !oge  T + z5  exp[z6 (T  - 273 . 1 5)]+  vliquid  ^ xl  0'3 

/=0  Ui 


where  z0  = -234.707325  z4  = 2.4188095 

z,  = -0.1 177858x1 0'3  z5  = -0.1976361xl0‘2 

z2  = 0.1501808xl0'6  zg=-0. 0865209 

zj  = - 0.8946646x1  O'10 

• Moist  supersaturated  air:  The  specific  entropy  per  unit  mass  of  mixture  of 

supersaturated  air  can  be  expressed  as  follows: 


WA+V,  + V^ 


(6-43) 


super 


m, 


super 


where  msuper  is  the  total  mass  of  the  mixture.  Using  a procedure  similar  to  the  one 

employed  to  determine  the  specific  volume  and  specific  enthalpy;  it  can  be  shown 
that  the  specific  entropy  of  supersaturated  air  per  unit  mass  of  dry  air  can  be 
expressed  as  follows: 

,''  sup  er  — Ss  ^ drop  l/f  up  er  ~ ^ s ) (6-44) 


where  sdrop  can  either  be  sliquid  or  sice  depending  on  the  air  temperature. 
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Specific  heat: 

■ Moist  subsaturated  air: 

+ WmCpw  (6-45) 

where  Cpa  and  Cpw  are  the  specific  heat  of  dry  air  and  water  vapor,  respectively,  and 
Wm  is  the  humidity  ratio  of  moist  air.  At  saturation  Cpm  becomes  Cps  and  Wm 
becomes  Ws . 

■ Moist  supersaturated  air: 

C„.~C„+C^r(Wm,-W,)  (6-46) 

where  Cps  is  the  specific  heat  of  air  at  saturated  conditions,  and  Cpdrop  is  the  specific 
heat  of  saturated  water,  which  can  either  be  CpHquid  or  Cpice  depending  on  the  air 
temperature. 

Psvchrometric  Formulations 

Mole  fraction  of  air  and  water  vapor.  Moist  air  is  considered  to  be  saturated  when  its 
relative  humidity  has  a value  of  100%.  In  this  state  the  mole  fraction  of  air  is  given  by 
xas , while  that  of  water  vapor  is  given  by  xm  . These  can  be  expressed  as 

r (6-47) 

« p 

_ =^_  (6-48) 

WS  p 

where  £ is  the  enhancement  factor  expressed  by  Equation  6-49,  Pws  is  the  partial  pressure 

of  water  vapor  at  saturation  and  P is  the  total  pressure.  The  enhancement  factor  (f)  is  a 
correction  parameter  that  takes  into  account  the  effect  of  dissolved  gases  and  pressure  on 
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the  properties  of  the  condensed  phase,  as  well  as  the  effect  of  intermolecular  forces  on  the 
properties  of  moisture  itself  (ASHRAE  2001).  Olivieri  et  al.  (1996)  proposed  the 
following  equation  to  determine  this  factor: 


loge(s)  = 


RT 


vc  + ln(l  - (ixasP)- 


1 n 

X P 


RT 


B. 


2 xJP 

R 

’(p-p,-xjp) 

R 

'lxJ{l-xJP2-Pj] 

RT 

^ aw 

RT 

2(rtJ 

+ 


x 3P2 

as  r 

r2t2 


c„ 


C - 


c„ 


(rtJ 


2 (RTJ 

~2xJ(2-3xJP 2 


BB - 


{rtJ 


BB 


_j_ 

6Xas2(l~Xas)2P2 

R R 

3xJP2 

D2  _ 

3-^flJ2(l-JCaJ)P2l 

1 

1^ 
1 

U / f 

ww  aw 

2 R2T2 

" aa 

1 

'sr 

"-rS 

i 

C 


2xJ(i-x,Jl-3xJP2' 

R2 

' Pj(Mx^~xJ  P2' 

1 

"aw 

2 (RTJ 

b: 


(6-49) 


where  vcis  the  specific  molar  volume  of  liquid  water  in  m3/kmol,  xas  and  xws  are  the 
mole  fractions  of  air  and  water  vapor  at  saturation,  respectively. 

Equation  6.49  must  be  solved  using  a numerical  method  because  xas  and  xws  are  both 

functions  of  e.  The  term  /?  is  the  Henry’s  Law  constant,  and  it  accounts  for  air  dissolved 
in  liquid  water.  This  term  can  be  expressed  in  (1/Pa)  as 


/?  = lxlO-4 — — — 
101325 


(6-50) 


where  [5a  is  the  Henry’s  law  constant  for  air.  This  term  can  be  determined  from 
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1 _ x°2  | xn2  (6-51) 

Pa  Poi  Pn2 


where  x(h  - 0.22  and  xN  = 0.78  are  the  mole  fractions  of  oxygen  and  nitrogen  in  the 

air.  Himmelblau  (1960)  gives  a set  of  equations  to  compute  fJ0?  and  fiN 

pQi  = exp(- 19.786773  + 23.2930487  -9.8984983T2  +2.2363 17273  (6-52) 

-0.2961843474  + 0.0170849375) 

and 


-288743.088  + 1 1 1 167.5947  + 
529382.752 


+ 4444.44537  + 


89505.6879 


+ 16.1931 1277 


491323.518 

7 

328398.259 

73 


- 27988. 296472 
-404.94098074 


(6-53) 


The  term  yis  the  isothermal  compressibility  of  liquid  water  and  it  can  be  determined  from 
(Nelson  et  al.  2001)  in  (1/Pa) 

7. 1201893  jc  1 0‘7  + 5.2000728  x 10_97  (6-54) 

(l  + 0.0 1 093297  - 0.3305  x 1 O'4  72  J 

where  7 <450  K and  where  the  pressures  are  less  than  IMPa. 

Equation  6-49  has  been  numerically  solved  for  a temperature  ranging  from  -60°C  to  60°C 
and  the  results  are  presented  in  Table  6-3.  In  order  to  propose  a new  simplified  equation 
for  the  enhancement  factor  the  data  presented  in  Table  6-3  are  fitted  into  a polynomial 
represented  by  Equation  6-55. 

s = 1.00391723  - 7.8220553  x 10~6t  + 6.94681594  x lO'V  (6-55) 

+ 3.040587  xlO'Y  -2.7852  xlO’V  -5.656 xl0"13t5 
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where  t is  the  temperature  in  °C. 


Temperature  (°C) 

Enhancement  factor  (s) 

-60 

1.0063 

-50 

1.0057 

-40 

1.0051 

-30 

1.0047 

-20 

1.0043 

-10 

1.0041 

0 

1.0039 

10 

1.0039 

20 

1.0041 

30 

1.0044 

40 

1.0048 

50 

1.0053 

60 

1.0058 

Air  humidity  ratio: 

■ Moist  subsaturated  air:  This  can  be  expressed  by 

sP 


W.  = 0.62198- 


P-sP, 


(6-56) 


where  Pw  is  the  partial  pressure  of  water  vapor  and  P is  the  total  pressure.  At 
saturation  Wm  becomes  Ws , and  Pw  becomes  Ps . 

■ Moist  supersaturated  air:  This  can  be  determined  from  Equation  6-16. 

Degree  of  saturation: 

■ Moist  subsaturated  air:  The  degree  of  saturation  (jud  )for  subsaturated  air  is  given  by 


Pd 


wm 
~ w. 


(6-57) 
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where  Wm  is  the  humidity  ratio  of  the  given  moist  air  sample  and  Ws  is  the  humidity 

ratio  of  saturated  air  at  the  same  temperature  and  pressure.  The  degree  of  saturation 
may  have  a value  between  0 and  1 for  air  in  the  subsaturated  zone. 

■ Moist  supersaturated  air:  The  degree  of  saturation  {judsuper) for  supersaturated  can  be 
expressed  as 

u (6-58) 

r*d  ,sup  er  yy 

r s 

The  degree  of  saturation  for  supersaturated  air  has  a value  larger  than  1 . 

Results 

A sample  collection  of  property  data  under  supersaturated  conditions  in  the 
temperature  range  -40°C  to  40°C  and  humidity  ratios  between  0 and  0.1  kg/kga  is  shown 
in  Table  6-4.  The  data  presented  include  the  specific  volume,  specific  enthalpy,  specific 
entropy,  and  humidity  ratio  for  degrees  of  saturation  from  100%  to  200%.  These 
properties  are  determined  at  atmospheric  pressure  and  are  presented  in  two  psychrometric 
charts:  one  for  temperatures  from  -40°C  to  0°C,  and  the  other  from  0°C  to  40°C  (Figures 
6-1  and  6-2).  Properties  of  saturated  air  determined  in  this  chapter  have  a maximum 
deviation  of  ±0.1%  compared  to  the  data  developed  by  Hyland  and  Wexler  (1983a,  b). 

Figure  6-1  shows  the  psychrometric  chart  for  temperatures  from  -40°C  to  0°C  and 
for  humidity  ratios  between  0 and  0.0076  kgw/kga,  while  Figure  6-2  shows  the 
psychrometric  chart  for  temperatures  from  0°C  to  40°C  and  for  humidity  ratios  between  0 
to  0.1  kgw/kga.  In  these  figures,  lines  of  constant  enthalpy,  degree  of  saturation, 
humidity  ratio,  dry-bulb  temperature,  and  specific  volume  are  represented  in  the 
supersaturated  zone  of  the  psychrometric  chart. 
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Figures  6-3  and  6-4  show  examples  of  state  points  placed  on  the  psychrometric 
chart.  For  example,  if  the  dry-bulb  temperature  is  16°C  and  the  degree  of  saturation  is 
150%  for  Point  1,  the  specific  enthalpy,  humidity  ratio,  and  the  specific  volume  can  be 
determined  using  Figure  6-1.  This  yields  a specific  enthalpy  of  -13.963  kJ/kg,  a humidity 
ratio  of  0.001396  kgw/kga,  and  specific  volume  of  0.729  m3/kg.  In  a similar  manner,  for  a 
dry-bulb  temperature  of  -4°C  and  a humidity  ratio  of  0.002979  kgw/kga  for  Point  2,  the 
other  properties  can  be  determined  from  the  same  chart.  This  yields  a degree  of 
saturation  of  110%,  a specific  enthalpy  of -2.635  kJ/kg,  and  a specific  volume  of  0.7653 
m /kg.  These  two  points  are  represented  in  Figure  6-3. 

Similarly,  using  Figure  6-2,  if  the  dry-bulb  temperature  is  8°C  and  the  degree  of 
saturation  is  120%  for  Point  3,  the  specific  enthalpy,  humidity  ratio,  and  the  specific 
volume  can  be  determined.  This  gives  a specific  enthalpy  of  -24.897  kJ/kg,  a humidity 
ratio  of  0.00802  kgw/kga,  and  a specific  volume  of  0.8046  m3/kg.  An  example  at  higher 
temperatures  is  represented  by  Point  4,  which  has  a dry-bulb  temperature  of  22°C  and  a 
humidity  ratio  of  0.02679  kgw/kga.  After  locating  this  point  on  the  chart,  the  properties 
can  be  read  as  degree  of  saturation  of  160%,  specific  enthalpy  of  -65.587  kJ/kg,  and 
specific  volume  of  0.853 1 m3/kg. 

Conclusions 

Psychrometric  data  for  supersaturated  air  have  been  developed  and  presented  in 
the  form  of  psychrometric  charts  for  degrees  of  saturation  from  100%  to  200%.  This 
method  extends  the  current  psychrometric  data,  which  are  presently  only  available  in  the 
subsaturated  zone  of  the  psychrometric  chart.  The  properties  developed  in  this 
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investigation  should  enable  the  refrigeration  engineer  and  system  designer  to  better 


predict  the  performance  of  the  system  as  well  as  design  a more  optimal  system. 

Table  6-4.  Thermodynamic  properties  of  moist  supersaturated  air  (standard  atmospheric 


pressure,  101325  Pa) 


Md  = 

100% 

Hd,  super 

= 120% 

t 

(°C) 

ws 

Vs 

is 

Kup  er 

V 

super 

1 super 

^sup  er 

(kgw/kga) 

(m3/kg) 

(kJ/kg) 

(kJ/kgK) 

(kg\v/kga) 

(m3/kg) 

(kJ/kg) 

(kJ/kg  K) 

-40 

0.0000793 

0.659700 

-40.229 

-0.1584 

0.0000952 

0.659700 

-40.236 

-0.1584 

-38 

0.0000992 

0.665400 

-37.976 

-0.1496 

0.0001190 

0.665400 

-37.984 

-0.1496 

-36 

0.0001237 

0.671200 

-35.905 

-0.1408 

0.0001484 

0.671200 

-35.915 

-0.1408 

-34 

0.0001536 

0.676900 

-33.820 

-0.1320 

0.0001843 

0.676900 

-33.832 

-0.1320 

-32 

0.0001902 

0.682600 

-31.718 

-0.1233 

0.0002282 

0.682600 

-31.733 

-0.1234 

-30 

0.0002346 

0.688400 

-29.597 

-0.1145 

0.0002815 

0.688400 

-29.615 

-0.1146 

-28 

0.0002883 

0.694100 

-27.454 

-0.1057 

0.0003460 

0.694100 

-27.476 

-0.1058 

-26 

0.0003533 

0.699900 

-25.282 

-0.0969 

0.0004240 

0.699900 

-25.309 

-0.0970 

-24 

0.0004314 

0.705700 

-23.078 

-0.0880 

0.0005177 

0.705700 

-23.111 

-0.0881 

-22 

0.0005251 

0.711500 

-20.834 

-0.0790 

0.0006301 

0.711500 

-20.874 

-0.0791 

-20 

0.0006373 

0.717300 

-18.545 

-0.0699 

0.0007648 

0.717300 

-18.593 

-0.0701 

-18 

0.0007711 

0.723100 

-16.201 

-0.0607 

0.0009253 

0.723100 

-16.258 

-0.0609 

-16 

0.0009303 

0.729000 

-13.793 

-0.0513 

0.0011164 

0.729000 

-13.861 

-0.0516 

-14 

0.0011191 

0.734900 

-11.311 

-0.0416 

0.0013429 

0.734900 

-11.392 

-0.0419 

-12 

0.0013425 

0.740900 

-8.742 

-0.0318 

0.0016110 

0.740900 

-8.838 

-0.0322 

-10 

0.0016062 

0.746900 

-6.072 

-0.0215 

0.0019274 

0.746900 

-6.186 

-0.0219 

-8 

0.0019166 

0.753000 

-3.283 

-0.0100 

0.0022999 

0.753000 

-3.417 

-0.0105 

-6 

0.0022811 

0.759100 

-0.357 

0.0000 

0.0027373 

0.759100 

-0.515 

-0.0006 

-4 

0.0027081 

0.765300 

2.728 

0.0115 

0.0032497 

0.765301 

2.543 

0.0108 

-2 

0.0032074 

0.771700 

5.995 

0.0236 

0.0038489 

0.771701 

5.778 

0.0228 

0 

0.0037895 

0.778100 

9.473 

0.0364 

0.0045474 

0.778101 

9.220 

0.0355 

0 

0.0037895 

0.778100 

9.473 

0.0364 

0.0045474 

0.778101 

9.479 

0.0364 

2 

0.0043810 

0.784500 

12.982 

0.0492 

0.0052572 

0.784501 

12.989 

0.0492 

4 

0.0050540 

0.791100 

16.696 

0.0627 

0.0060648 

0.791101 

16.713 

0.0628 

6 

0.0058180 

0.797800 

20.644 

0.0769 

0.0069816 

0.797801 

20.673 

0.0770 

8 

0.0066830 

0.804600 

24.852 

0.0919 

0.0080196 

0.804601 

24.897 

0.0921 

10 

0.0076610 

0.811600 

29.352 

0.1078 

0.0091932 

0.811602 

29.416 

0.1080 

12 

0.0087660 

0.818800 

34.179 

0.1248 

0.0105192 

0.818802 

34.267 

0.1251 

14 

0.0100120 

0.826200 

39.370 

0.1430 

0.0120144 

0.826202 

39.488 

0.1434 

16 

0.0114130 

0.833800 

44.963 

0.1620 

0.0136956 

0.833802 

45.116 

0.1625 

18 

0.0129890 

0.841700 

51.008 

0.1830 

0.0155868 

0.841703 

51.204 

0.1837 

20 

0.0147580 

0.849800 

57.555 

0.2060 

0.0177096 

0.849803 

57.803 

0.2069 

22 

0.0167410 

0.858300 

64.660 

0.2300 

0.0200892 

0.858303 

64.969 

0.2311 

24 

0.0189630 

0.867100 

72.385 

0.2560 

0.0227556 

0.867104 

72.767 

0.2573 

26 

0.0214480 

0.876400 

80.798 

0.2840 

0.0257376 

0.876404 

81.266 

0.2856 

28 

0.0242260 

0.886000 

89.976 

0.3150 

0.0290712 

0.886005 

90.545 

0.3170 

30 

0.0273290 

0.896200 

100.006 

0.3480 

0.0327948 

0.896205 

100.694 

0.3504 

32 

0.0307930 

0.909600 

110.979 

0.3840 

0.0369516 

0.909606 

111.805 

0.3869 

34 

0.0346600 

0.918300 

123.011 

0.4240 

0.0415920 

0.918307 

123.999 

0.4274 

36 

0.0389710 

0.930300 

136.209 

0.4670 

0.0467652 

0.930308 

137.384 

0.4710 

38 

0.0437780 

0.943100 

150.713 

0.5140 

0.0525336 

0.943109 

152.107 

0.5188 

40 

0.0491410 

0.956800 

166.683 

0.5650 

0.0589692 

0.956810 

168.329 

0.5706 
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Table  6-4.  Continued. 


super 

= 140% 

^dy  super 

= 160% 

t 

(°C) 

^sup  er 

V 

super 

*super 

“^sup  er 

^sup  er 

v 

K super 

*sup  er 

C 

°sup  er 

(kgw/kga) 

(m3/kg) 

(kJ/kg) 

(kJ/kg  K) 

(kgw/kga) 

(m3/kg) 

(kJ/kg) 

(kJ/kg  K) 

-40 

0.0001110 

0.659700 

-40.242 

-0.1584 

0.0001269 

0.659700 

-40.249 

-0.1585 

-38 

0.0001389 

0.665400 

-37.992 

-0.1497 

0.0001587 

0.665400 

-38.000 

-0.1497 

-36 

0.0001732 

0.671200 

-35.925 

-0.1409 

0.0001979 

0.671200 

-35.935 

-0.1409 

-34 

0.0002150 

0.676900 

-33.845 

-0.1321 

0.0002458 

0.676900 

-33.857 

-0.1321 

-32 

0.0002663 

0.682600 

-31.748 

-0.1234 

0.0003043 

0.682600 

-31.763 

-0.1235 

-30 

0.0003284 

0.688400 

-29.634 

-0.1146 

0.0003754 

0.688400 

-29.652 

-0.1147 

-28 

0.0004036 

0.694100 

-27.499 

-0.1059 

0.0004613 

0.694100 

-27.521 

-0.1059 

-26 

0.0004946 

0.699900 

-25.336 

-0.0971 

0.0005653 

0.699900 

-25.364 

-0.0972 

-24 

0.0006040 

0.705700 

-23.144 

-0.0882 

0.0006902 

0.705700 

-23.177 

-0.0884 

-22 

0.0007351 

0.711500 

-20.913 

-0.0793 

0.0008402 

0.711500 

-20.953 

-0.0794 

-20 

0.0008922 

0.717300 

-18.640 

-0.0703 

0.0010197 

0.717300 

-18.688 

-0.0704 

-18 

0.0010795 

0.723100 

-16.315 

-0.0611 

0.0012338 

0.723101 

-16.372 

-0.0613 

-16 

0.0013024 

0.729000 

-13.929 

-0.0518 

0.0014885 

0.729001 

-13.997 

-0.0521 

-14 

0.0015667 

0.734900 

-11.473 

-0.0422 

0.0017906 

0.734901 

-11.554 

-0.0425 

-12 

0.0018795 

0.740901 

-8.934 

-0.0325 

0.0021480 

0.740901 

-9.030 

-0.0329 

-10 

0.0022487 

0.746901 

-6.300 

-0.0223 

0.0025699 

0.746901 

-6.413 

-0.0228 

-8 

0.0026832 

0.753001 

-3.551 

-0.0110 

0.0030666 

0.753001 

-3.685 

-0.0115 

-6 

0.0031935 

0.759101 

-0.673 

-0.0012 

0.0036498 

0.759101 

-0.830 

-0.0017 

-4 

0.0037913 

0.765301 

2.358 

0.0101 

0.0043330 

0.765302 

2.173 

0.0095 

-2 

0.0044904 

0.771701 

5.562 

0.0220 

0.0051318 

0.771702 

5.345 

0.0212 

0 

0.0053053 

0.778102 

8.968 

0.0345 

0.0060632 

0.778102 

8.715 

0.0336 

0 

0.0053053 

0.778102 

9.486 

0.0364 

0.0060632 

0.778102 

9.492 

0.0364 

2 

0.0061334 

0.784502 

12.997 

0.0493 

0.0070096 

0.784503 

13.004 

0.0493 

4 

0.0070756 

0.791102 

16.730 

0.0628 

0.0080864 

0.791103 

16.747 

0.0629 

6 

0.0081452 

0.797802 

20.703 

0.0771 

0.0093088 

0.797803 

20.732 

0.0772 

8 

0.0093562 

0.804603 

24.942 

0.0922 

0.0106928 

0.804604 

24.987 

0.0924 

10 

0.0107254 

0.811603 

29.481 

0.1083 

0.0122576 

0.811605 

29.545 

0.1085 

12 

0.0122724 

0.818804 

34.356 

0.1254 

0.0140256 

0.818805 

34.444 

0.1257 

14 

0.0140168 

0.826204 

39.605 

0.1438 

0.0160192 

0.826206 

39.723 

0.1443 

16 

0.0159782 

0.833805 

45.270 

0.1631 

0.0182608 

0.833807 

45.423 

0.1636 

18 

0.0181846 

0.841705 

51.401 

0.1844 

0.0207824 

0.841708 

51.597 

0.1851 

20 

0.0206612 

0.849806 

58.051 

0.2077 

0.0236128 

0.849809 

58.298 

0.2086 

22 

0.0234374 

0.858307 

65.278 

0.2322 

0.0267856 

0.858310 

65.587 

0.2333 

24 

0.0265482 

0.867108 

73.149 

0.2587 

0.0303408 

0.867111 

73.531 

0.2600 

26 

0.0300272 

0.876409 

81.733 

0.2873 

0.0343168 

0.876413 

82.201 

0.2889 

28 

0.0339164 

0.886010 

91.114 

0.3190 

0.0387616 

0.886015 

91.682 

0.3209 

30 

0.0382606 

0.896211 

101.381 

0.3528 

0.0437264 

0.896216 

102.069 

0.3552 

32 

0.0431102 

0.909612 

112.631 

0.3897 

0.0492688 

0.909619 

113.457 

0.3926 

34 

0.0485240 

0.918314 

124.987 

0.4308 

0.0554560 

0.918321 

125.974 

0.4342 

36 

0.0545594 

0.930316 

138.560 

0.4751 

0.0623536 

0.930324 

139.735 

0.4791 

38 

0.0612892 

0.943118 

153.501 

0.5236 

0.0700448 

0.943126 

154.895 

0.5283 

40 

0.0687974 

0.956820 

169.975 

0.5763 

0.0786256 

0.956830 

171.622 

0.5819 
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Table  6-4.  Continued. 


^d, super 

= 180% 

r*dy  super 

= 200% 

t 

(°C) 

^sup  er 

V 

sup  er 

huper 

C 

° super 

Wsup  er 

V 

super 

^super 

C 

° super 

(kgw/kga) 

(nrVkg) 

(kJ/kg) 

(kJ/kg  K) 

(kgw/kga) 

(m3/kg) 

(kJ/kg) 

(kJ/kg  K) 

-40 

0.0001427 

0.659700 

-40.255 

-0.1585 

0.0001586 

0.659700 

-40.262 

-0.1585 

-38 

0.0001786 

0.665400 

-38.008 

-0.1497 

0.0001984 

0.665400 

-38.016 

-0.1498 

-36 

0.0002227 

0.671200 

-35.945 

-0.1409 

0.0002474 

0.671200 

-35.955 

-0.1410 

-34 

0.0002765 

0.676900 

-33.869 

-0.1322 

0.0003072 

0.676900 

-33.882 

-0.1322 

-32 

0.0003424 

0.682600 

-31.778 

-0.1235 

0.0003804 

0.682600 

-31.794 

-0.1236 

-30 

0.0004223 

0.688400 

-29.671 

-0.1148 

0.0004692 

0.688400 

-29.689 

-0.1148 

-28 

0.0005189 

0.694100 

-27.544 

-0.1060 

0.0005766 

0.694100 

-27.566 

-0.1061 

-26 

0.0006359 

0.699900 

-25.391 

-0.0973 

0.0007066 

0.699900 

-25.418 

-0.0974 

-24 

0.0007765 

0.705700 

-23.210 

-0.0885 

0.0008628 

0.705700 

-23.243 

-0.0886 

-22 

0.0009452 

0.711500 

-20.993 

-0.0796 

0.0010502 

0.711501 

-21.032 

-0.0797 

-20 

0.0011471 

0.717301 

-18.736 

-0.0706 

0.0012746 

0.717301 

-18.783 

-0.0708 

-18 

0.0013880 

0.723101 

-16.429 

-0.0615 

0.0015422 

0.723101 

-16.486 

-0.0617 

-16 

0.0016745 

0.729001 

-14.065 

-0.0523 

0.0018606 

0.729001 

-14.134 

-0.0526 

-14 

0.0020144 

0.734901 

-11.635 

-0.0428 

0.0022382 

0.734901 

-11.716 

-0.0431 

-12 

0.0024165 

0.740901 

-9.127 

-0.0332 

0.0026850 

0.740901 

-9.223 

-0.0336 

-10 

0.0028912 

0.746901 

-6.527 

-0.0232 

0.0032124 

0.746902 

-6.641 

-0.0236 

-8 

0.0034499 

0.753002 

-3.820 

-0.0120 

0.0038332 

0.753002 

-3.954 

-0.0125 

-6 

0.0041060 

0.759102 

-0.988 

-0.0023 

0.0045622 

0.759102 

-1.146 

-0.0029 

-4 

0.0048746 

0.765302 

1.987 

0.0088 

0.0054162 

0.765303 

1.802 

0.0081 

-2 

0.0057733 

0.771703 

5.129 

0.0204 

0.0064148 

0.771703 

4.912 

0.0196 

0 

0.0068211 

0.778103 

8.462 

0.0327 

0.0075790 

0.778104 

8.210 

0.0318 

0 

0.0068211 

0.778103 

9.498 

0.0364 

0.0075790 

0.778104 

9.505 

0.0364 

2 

0.0078858 

0.784504 

13.011 

0.0493 

0.0087620 

0.784504 

13.019 

0.0493 

4 

0.0090972 

0.791104 

16.764 

0.0629 

0.0101080 

0.791105 

16.781 

0.0630 

6 

0.0104724 

0.797805 

20.761 

0.0773 

0.0116360 

0.797806 

20.790 

0.0774 

8 

0.0120294 

0.804605 

25.032 

0.0925 

0.0133660 

0.804607 

25.076 

0.0927 

10 

0.0137898 

0.811606 

29.609 

0.1087 

0.0153220 

0.811608 

29.674 

0.1090 

12 

0.0157788 

0.818807 

34.532 

0.1261 

0.0175320 

0.818809 

34.621 

0.1264 

14 

0.0180216 

0.826208 

39.841 

0.1447 

0.0200240 

0.826210 

39.959 

0.1451 

16 

0.0205434 

0.833809 

45.576 

0.1642 

0.0228260 

0.833811 

45.730 

0.1647 

18 

0.0233802 

0.841710 

51.793 

0.1858 

0.0259780 

0.841713 

51.989 

0.1865 

20 

0.0265644 

0.849812 

58.546 

0.2095 

0.0295160 

0.849815 

58.794 

0.2104 

22 

0.0301338 

0.858313 

65.896 

0.2344 

0.0334820 

0.858317 

66.205 

0.2354 

24 

0.0341334 

0.867115 

73.913 

0.2614 

0.0379260 

0.867119 

74.295 

0.2627 

26 

0.0386064 

0.876417 

82.668 

0.2905 

0.0428960 

0.876422 

83.136 

0.2922 

28 

0.0436068 

0.886019 

92.251 

0.3229 

0.0484520 

0.886024 

92.820 

0.3249 

30 

0.0491922 

0.896222 

102.756 

0.3575 

0.0546580 

0.896227 

103.444 

0.3599 

32 

0.0554274 

0.909625 

114.282 

0.3954 

0.0615860 

0.909631 

115.108 

0.3983 

34 

0.0623880 

0.918328 

126.962 

0.4376 

0.0693200 

0.918335 

127.950 

0.4410 

36 

0.0701478 

0.930331 

140.910 

0.4832 

0.0779420 

0.930339 

142.086 

0.4872 

38 

0.0788004 

0.943135 

156.289 

0.5331 

0.0875560 

0.943144 

157.682 

0.5379 

40 

0.0884538 

0.956840 

173.268 

0.5875 

0.0982820 

0.956850 

174.914 

0.5931 
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Figure  6-1.  Psychrometric  chart  in  the  supersaturated  zone  for  a temperature  range  from 
-40  to  0 °C 
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Figure  6-2.  Psychrometric  chart  in  the  supersaturated  zone  for  a temperature  range  from 
0 to  40  °C 
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Figure  6-3.  Psychrometric  chart  in  the  supersaturated  zone  for  a temperature  range  from 
-40  to  0 °C,  showing  two  sample  points 
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Figure  6-4.  Psychrometric  chart  in  the  supersaturated  zone  for  a temperature  range  from 
0 to  40  °C,  showing  two  sample  points 
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CHAPTER  7 

HEAT  AMD  MASS  TRANSFER  ON  A CYLINDER  SURFACE  IN  CROSS  FLOW 
UNDER  SUPERSATURATED  FROSTING  CONDITIONS 

Objectives 

This  chapter  presents  a calculation  procedure  for  a cylinder  surface  exposed  to 
supersaturated  air.  In  this  procedure  the  psychrometric  formulations  for  supersaturated 
moist  air  will  be  used  to  predict  frost  formation  and  heat  transfer  rates  on  the  cylinder 
surface  in  question.  The  cylinder  surface  was  selected  because  it  could  represent  the 
tube  in  a freezer  coil.  While  this  is  not  identical  to  a finned-tube  freezer  coil  per  se,  it 
serves  as  an  illustration  of  how  the  techniques  developed  here  can  be  utilized  to  predict 
the  performance  of  a real  freezer  coil  in  the  supersaturated  zone  of  the  psychrometric 
chart.  This  geometry  has  also  been  selected  so  that  results  of  the  supersaturated  model 
could  be  compared  to  other  published  works  in  the  subsaturated  zone  of  the 
psychrometric  chart,  such  as  Andrichak  (1962),  Chung  and  Algren  (1958a,  b),  and  Parish 
and  Sepsy  (1972),  among  others.  This  comparison  is  crucial,  as  the  characteristics  of  the 
frost-formation  process  under  supersaturated  conditions  should  be  referenced  to  those  of 
a frosted  surface  with  a predictable  performance  (e.g.,  under  subsaturated  conditions). 

The  calculation  procedure  developed  here  simulates  the  frost-formation  process  as 
a quasi-steady-state  one.  In  a quasi-steady-state  approach,  time-dependent  phenomena 
are  treated  as  steady-state  ones  for  sufficiently  small  time  intervals.  In  other  words,  the 
overall  solution  is  treated  as  a collection  of  a series  of  steady-state  solutions  at  different 
time  increments.  So,  while  only  steady-state  equations  are  employed,  the  overall  solution 
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is  time-dependent.  The  accuracy  of  this  approach  can  be  almost  as  high  as  the  real-time 
solution  of  the  time-dependent  problem  provided  that  the  time  increment  chosen  is 
sufficiently  small.  The  size  of  the  time  step  is  experimented  with  and  sensitivity  of  the 
solution  to  the  step  size  will  be  established  and  checked. 

The  problem  of  frost  formation  on  cylinders  in  cross  flow  is  a complicated  one 
because  of  the  following  factors: 

• The  boundary  layer  solution  is  transient.  The  thermal  resistance  of  the  frost  layer 
changes  with  both  time  and  position  when  the  frost  layer  is  growing  in  thickness. 
Therefore  the  frost  surface  temperature  also  changes  with  time  and  position.  The 
thermal  and  diffusion  boundary  layers  also  change,  resulting  in  a change  in  the  heat 
and  mass  transfer  coefficients  with  time. 

• The  frost  density  and  thermal  conductivity  (parameters  needed  in  the  modeling)  vary 
continuously  during  the  development  of  the  frost  layer,  although  eventually  a quasi- 
steady state  is  achieved. 

• The  frost-formation  process  involves  simultaneous  heat  and  mass  transfer.  The 
model  presented  in  this  chapter  employs  empirical  correlations  for  heat  and  mass 
transfer  for  flow  past  a circular  cylinder.  The  frost  layer  is  not  uniform  around  the 
cylinder  during  the  formation  process,  which  complicates  the  situation  due  to  the  fact 
that  the  geometry  is  no  longer  a cylindrical  one. 

• There  is  no  unique  behavior  for  the  heat  and  mass  transfer  coefficients  past  the 
separation  point  (Giedt  1949).  Thus  the  complexity  of  the  modeling  effort  is 
expected  to  be  significantly  higher  in  that  region. 
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For  the  reasons  mentioned  above,  modeling  of  this  problem  will  be  carried  out 
using  two  methods.  First,  the  frost  will  be  assumed  to  grow  uniformly  around  the 
cylinder,  thus  neglecting  the  angular  dependency  of  the  variables.  This  allows  the 
employment  of  empirical  correlations  for  the  average  heat  and  mass  transfer  coefficients. 
The  second  method  accounts  for  the  angular  dependency  using  local  properties  as  well  as 
the  local  heat  and  mass  transfer  coefficients.  This  part  of  the  analysis  will  be  performed 
only  from  the  stagnation  point  to  the  separation  point.  This  would  allow  the  prediction  of 
the  spatial  and  temporal  variations  of  the  frost  thickness  and  temperature. 

Analysis 

As  was  mentioned  earlier,  the  analysis  presented  here  is  a quasi-steady  state  one, 
whereby  steady  state  conditions  are  allowed  to  exist  for  sufficiently  small  intervals  of 
time.  The  analysis  presented  will  illustrate  a mechanism  for  computing  the  frost 
formation  and  heat  transfer  rates  as  functions  of  both  time  and  angular  position  on  the 
cylinder  surface.  Use  is  made  of  empirical  correlations  for  the  heat  transfer  coefficient 
and  for  frost  properties.  A section  of  the  cylinder  used  is  shown  in  Figure  7-1. 

The  energy  balance  at  any  angular  position  on  a small  area  segment  on  the 
cylinder  surface  (Figure  7-1)  in  the  vicinity  of  the  angular  position  (0),  can  be  written  as 
follows: 

SQconv.d  + 3Qdiff,e  ~ 8Qcond,e  = 0 (7’  1 ) 

where  8Qcomfi , 8Qdiff,e » and  ^QCOnd,e  are  the  heat  transfer  rates  by  convection,  due  to  mass 
diffusion,  and  by  conduction,  respectively,  to  the  finite  segment  of  the  cylinder  surface  in 
the  vicinity  of  the  angular  position  6. 
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Figure  7-1.  Small  area  segment  on  the  cylinder  surface  in  the  vicinity  of  the  angular 
position  ( 0 ) 

The  rate  of  heat  convected  can  be  determined  from 

=K,.,^.->p)se  (7-2) 

where  hc  o 9 , tfg , and  rfe  are  the  heat  transfer  coefficient,  the  frost  surface  temperature, 
and  frosted  cylinder  radius  at  angular  position  9,  respectively,  and  where  A and  tm  are 

the  cylinder  length  and  the  temperature  of  the  free  stream  moist  air,  respectively.  The 
quantity  SO  represents  a small  but  finite  angular  deviation  around  6. 

The  local  Nusselt  number  (Nu0 ) and  heat  transfer  coefficient  [hcfi  0 ) over  the 

forward  portion  of  the  cylinder  (0<6*<80o)  can  be  determined  employing  the  equation 
proposed  by  Martinelli  et  al.  (1943) 


Nu , 


_ K,o,e{^rfe)  _ | j 


14  Re°*m„Pr. 


0.4 

Df, super*  'super 


sup  er 


1- 


f e ^ 


V90  , 


(7-3) 


where  ksuper  is  the  thermal  conductivity  of  supersaturated  air,  while  ReDf  swer  and 


Df, super  ' 


TV  are  the  Reynolds  and  Prandtl  numbers  evaluated  under  supersaturated  conditions, 


respectively. 
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The  heat  transfer  rate  due  to  mass  diffusion  is  given  by 

time  = - V,Aw  k - )»  P-4) 

where  hmass6  and  Wfe  are  the  local  mass  transfer  coefficient  and  the  humidity  ratio  of  air 
evaluated  at  tfe  (assuming  saturation  conditions  in  the  vicinity  of  the  air-frost  interface), 
respectively,  and  where  Lse  and  Wm  are  the  latent  heat  of  sublimation  of  the  frost  and 
the  humidity  ratio  of  the  free  stream  moist  air. 

The  humidity  ratio  of  the  free  stream  moist  air  (Wm) can  be  determined  from 


(ASHRAE  2001) 


W = = 0.62198-  e(P'®  ('7'5^ 

ma  P-e(Pj) 

where  mw  and  ma  are  the  mass  of  water  vapor  in  air  and  the  mass  of  dry  air,  respectively, 
Ps  is  the  saturation  pressure  of  the  water  vapor,  P is  the  total  pressure,  </>  is  the  relative 

humidity  of  the  air,  and  e is  the  enhancement  factor,  which  can  be  determined  from 
Equation  6-52. 

Equation  7-4  can  be  modified  to  determine  the  diffusion-driven  heat  transfer  rate 
under  supersaturated  conditions  replacing  Wm  with  the  humidity  ratio  of  supersaturated 
air  i}Vsuper),  which  can  be  expressed  as  follows: 


W 

super 


_ mWs+mdrop 

m„ 


= W.+- 


m 


drop 


(7-6) 


m„ 


where  (fVs)  is  the  humidity  ratio  of  air  at  saturation  and  md  is  the  mass  of  excess 


moisture  beyond  saturation  (either  liquid  or  ice,  depending  on  whether  the  dry  bulb 
temperature  of  the  air  is  above  or  below  the  water  freezing  point,  respectively). 
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The  local  latent  heat  of  sublimation  for  frost  can  be  determined  from  (modified 
from  Parish  1970) 

Ls0  = [-0.1 083(l  .8  tfg  + 32)+  2833]  (7-7) 

where  Ls0  in  Equation  7-7  is  expressed  in  J/kg. 

The  local  mass  transfer  coefficient  at  any  angular  position  0 can  be  determined  as 
follows: 


h = 

mass ,0  j-  2/3 

^ p,super^e 


(7-8) 


The  Lewis  number  is  a relevant  parameter  to  situations  involving  simultaneous 
heat  and  mass  transfer.  It  is  defined  as 


Le=a**er 


(7-9) 


D 


diff 


The  mass  diffusion  coefficient  can  be  determined  in  terms  of  the  molecular 


weights  and  molar  volumes  of  the  dry  air  and  water  vapor  (Holman  1 997) 


Ddiff  = 0.04357  (rj/2 

Uv^+v-)2 


(7-10) 


where  Ddiff  is  in  m2/s,  Patm  is  in  Pa,  Ma  and  Mv  are  the  molecular  weights  of  dry  air 


(28.9)  and  water  vapor  (18),  respectively,  and,  Vaand  Vvare  the  molar  volumes  of  dry 
air  (29.9)  and  water  vapor  (18.8),  respectively. 

The  rate  of  heat  transfer  by  conduction  is  determined  using  the  Fourier’s  heat 
conduction  equation 


& 


cond  ,i 9 


kfgMffe  tp ) 
L°Se{rfe/rp) 


-SO 


(7-11) 


126 


where  kfe  is  the  local  thermal  conductivity  of  frost. 

There  are  several  equations  to  determine  the  frost  density  and  thermal 
conductivity,  but  for  the  purpose  of  this  analysis  only  two  of  them  are  used.  For  the  frost 
density  [pfe) an  empirical  expression  proposed  by  Hayashi  et  al.  (1977)  is  applied.  This 

expression  is  a sole  function  of  the  frost  surface  temperature  as  follows: 

pfg  =65Oexp(o.227t/0)  (7-12) 

where  pfe  is  in  kg/m3.  This  equation  should  be  applied  under  the  following  conditions: 

frost  surface  temperature  between  -25°C  and  0°C;  airstream  velocities  between  2m/s  and 
6 m/s,  and  airstream  humidity  ratio  equal  to  0.0075  kgw/kga.  While  these  conditions  are  a 
bit  restrictive  for  supersaturated  air,  it  was  felt  that  the  error  resulting  from  using  the 
correlation  at  slightly  higher  humidity  ratios  than  the  ones  suggested  by  Hayashi  et  al. 
(1977)  may  be  partially  compensated  for  by  the  fact  that  the  frost  surface  temperature  in  a 
supersaturated  air  scenario  would  necessarily  be  higher  than  the  corresponding  frost 
surface  temperature  in  a subsaturated  air  scenario,  assuming  everything  else  being 
unchanged.  Higher  frost  surface  temperatures  produce  higher  frost  densities,  a fact 
consistent  with  experimental  observations  made  by  Sherif  et  al.  (2002)  for  snow-like  frost 
resulting  from  a supersaturated  air  scenario.  It  is  thus  argued  that  the  manner  with  which 
the  density  is  made  dependent  on  the  frost  surface  temperature  in  Equation  7-12  should 
make  the  use  of  the  equation  justified. 

For  the  frost  thermal  conductivity  ( kfg ),  an  empirical  correlation  proposed  by 
Yonko  and  Sepsy  (1967)  is  applied.  This  correlation  is  given  by  the  following  equation: 
kfe  =0.024248  + (0.00072311/^)+  (o.OOOOO  11 83p/ff2)  (7-13) 
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where  pfe  is  in  kg/m3  and  kfg  is  expressed  in  W/(m-K).  The  above  equation  is  only 

valid  for  frost  densities  less  than  573  kg/m3.  This  limiting  value  of  the  frost  density  is 
likely  to  be  a good  upper  bound  for  some  of  the  most  severe  scenarios  observed  in  a 
parallel  experimental  investigation  performed  by  the  authors  when  supersaturated  frost 
deposition  was  taking  place.  Thus,  it  was  felt  that  the  use  of  Equation  7-13  is  justified. 
The  local  frost  deposition  rate  per  unit  area  can  now  be  computed  using 

™f,e=hmassMm-Wfe)  CM 4) 

While  the  local  increase  in  frost  thickness  can  be  computed  from 

AV  _"V,*Ar  (7-15) 
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where  Aris  the  time  increment  during  which  steady  state  conditions  were  assumed  to 
exist  to  make  the  quasi-steady  state  assumption  valid. 

The  new  frost  thickness  at  the  new  time  is  related  to  that  at  the  previous  time 
according  to  the  equation 


^ ffi  ~ X fjg  + ffi 


(7-16) 


The  frost  surface  temperature  can  now  be  obtained  by  substituting  Equations 
7-2,  7-4,  and  7-1 1 into  Equation  7-1  to  yield 

^ _ Kofif/e L°ge irfe !ri‘ kf0tp  + k er^eh  s{fVm  — Wj-S ^Loge [r^ //y ) (7-17) 
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The  calculation  procedure  can  proceed  as  follows:  the  time  and  the  frost 

thickness  (Xf0)  are  first  initialized  to  zero,  and  the  frost  surface  temperature  (tfS) is  set 


equal  to  that  of  the  cylinder  surface  (tp ) . Air  properties  are  evaluated  at  the  air-frost 
boundary  layer  film  temperature.  Using  Equations  7-3  and  7-8  the  local  heat  and  mass 
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transfer  coefficients  are  computed.  Equations  7-12,  7-13,  and  7-14  are  then  used  to 
determine  the  frost  density,  frost  thermal  conductivity,  and  the  mass  deposition  rate  per 
unit  area,  respectively.  The  increase  in  frost  thickness  and  the  new  frost  thickness  can 
then  be  computed  using  Equations  7-15  and  7-16,  respectively.  Using  the  parameters 
determined  above,  a first  approximation  for  the  frost  surface  temperature  is  determined 
employing  Equation  7-17.  This  approximate  value  of  tfg,  is  then  used  to  determine  the 

new  boundary  layer  film  temperature  in  order  to  re-compute  the  local  heat  and  mass 
transfer  coefficients,  the  local  frost  density  and  thermal  conductivity,  as  well  as  a new 
approximate  value  for  the  frost  surface  temperature.  This  calculation  procedure  is 
repeated  until  the  value  of  the  new  frost  surface  temperature  computed  is  within  a 
difference  of  0.000 1°C  of  the  previous  value.  A numerical  procedure  describing  the  flow 
of  calculations  is  illustrated  in  Figure  7-2. 

A more  simplified  set  of  equations,  whereby  the  angular  dependency  of  the 
variables  is  neglected  (implying  that  the  frost  is  uniformly  building  up  on  the  surface), 
can  be  employed  if  the  average  heat  and  mass  transfer  coefficients  are  used  in  lieu  of  the 
local  coefficients.  For  example,  Churchill  and  Bernstein  (1977)  proposed  the  following 
formula  for  the  average  Nusselt  number: 


hcoDf 

c.o  J 


(7-18) 


This  equation  is  recommended  for  (^eD/sup<,,.EVsuper  )>0.2,  where  all  the  properties 


are  evaluated  at  the  film  temperature  between  the  air  and  the  frost  surface.  All  other 
equations  remain  the  same  provided  that  variables  with  angle  dependency  are  replaced 
with  their  average  values  integrated  over  the  whole  cylinder. 
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Figure  7-2.  Flow  chart  of  the  numerical  procedure  used 
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Results  and  Discussion 

The  procedure  described  in  the  previous  section  was  applied  to  cylinders  under 
forced  convection  for  both  subsaturated  and  supersaturated  air,  and  comparisons  with 
experimental  data  were  made  whenever  possible.  However,  before  generating  results 
from  the  model,  the  stability  of  the  solution  was  checked  in  terms  of  its  dependence  on 
the  time  step.  This  check  is  important  as  we  do  not  want  the  results  to  depend  on  the  size 
of  the  time  step.  That  dependence  is  shown  in  Figure  7-3,  which  shows  the  variation  of 
the  frost  thickness  with  time  for  three  values  of  the  time  step.  As  can  be  seen,  there  are 
no  discemable  differences  among  the  three  cases. 


Figure  7-3.  Effect  of  the  time  step  on  the  stability  of  the  solution:  — Ax=  Is, zlz=3s, 

and  — Az=5s 

For  all  further  results  shown  in  this  chapter,  the  time  step  used  is  one  second. 
This  means  that  steady  state  conditions  were  assumed  to  exist  for  periods  of  one  second 
at  a time.  Once  the  solution  converged  at  a particular  point  in  time,  the  generated  results 
were  used  as  initial  conditions  for  the  subsequent  point  in  time. 
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Figures  7-4  and  7-5  show  the  variation  of  the  frost  thickness  with  time  for  both 
subsaturated  and  supersaturated  air  at  the  stagnation  point  and  at  the  separation  point 
(corresponding  to  angular  position  of  80°),  respectively.  Results  of  the  present  model  are 
compared  with  those  of  Andrichak  (1962)  and  Parish  and  Sepsy  (1972)  for  the 
subsaturated  cases.  As  can  be  seen,  results  from  the  present  model  agree  well  with  those 
of  other  researchers.  The  difference  in  thickness  for  supersaturated  and  subsaturated 
conditions  is  notable,  considering  the  small  span  of  time  during  which  frost  thickness  is 
computed.  Comparing  Figures  7-4  and  7-5,  it  is  easy  to  observe  that  the  frost  thickness 
at  the  stagnation  point  is  higher  than  that  at  the  separation  point,  with  all  other  variables 
remaining  unchanged. 
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Figure  7-4.  Variation  of  the  frost  thickness  with  time  at  the  stagnation  point:  — present 
work  (subsaturated  conditions),  ■ Adrichack  (1962),  o Parish  and  Sepsy 
(1972),  — present  work  (supersaturated  conditions) 
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Figure  7-5.  Variation  of  the  frost  thickness  with  time  at  the  separation  point:  — present 
work  (subsaturated  conditions),  ■ Adrichack  (1962),  o Parish  and  Sepsy 
(1972),  — present  work  (supersaturated  conditions) 

Figures  7-6  and  7-7  show  the  variation  of  the  frost  surface  temperature  and  frost 
thickness  with  angular  position  from  the  stagnation  point  to  the  separation  point, 
respectively,  for  two  different  times.  Results  of  the  present  model  are  compared  with 
those  of  Raju  and  Sherif  (1993)  for  times  of  3 and  6 minutes.  As  can  be  seen,  results 
from  the  present  model  agree  well  with  those  of  Raju  and  Sherif  (1993).  Other 
researchers  have  shown  that  the  heat  and  mass  transfer  rates  decrease  as  the  angular 
position  changes  from  the  stagnation  point  to  the  separation  point  (Schmidt  and  Wenner 
1941  and  Giedt  1949). 

This  study  also  includes  results  where  the  angular  dependency  of  the  variables  is 
neglected  (implying  that  the  frost  is  uniformly  building  up  on  the  surface).  Figure  7-8 
shows  the  effect  of  the  Reynolds  number  on  the  frost  thickness  for  subsaturated  and 
supersaturated  air  scenarios. 
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Figure  7-6.  Variation  of  the  frost  surface  temperature  with  the  angle  ( 9 ) for  different 

times:  a)  3 min:  — present  work  and  ■ Raju  and  Sherif  (1993),  b)  6 min: 
— present  work  and  o Raju  and  Sherif  (1993) 


Figure  7-7.  Variation  of  the  frost  thickness  with  the  angle  (0)  for  different  times:  a)  3 
min:  — present  work  and  b)  6 min:  — present  work 
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Figure  7-8.  Effect  of  the  Reynolds  number  on  the  frost  thickness:  a)  Re=2592: 
— subsaturated  conditions,  o supersaturated  conditions,  and  b)  Re=9965: 
▲ subsaturated  conditions  and  — supersaturated  conditions 

A close  examination  of  the  figure  reveals  the  significant  differences  between  the 
frost  thickness  at  a certain  point  in  time  for  subsaturated  air  versus  supersaturated  air.  It 
is  natural  to  expect  that  there  would  be  more  frost  when  the  air  is  supersaturated.  It  is 
also  natural  to  expect  a larger  thickness  for  higher  Reynolds  numbers.  Observations 
made  earlier  in  this  investigation  confirm  these  results.  Sherif  et  al.  (1993)  made  some 
interesting  observations  regarding  the  Reynolds  number  dependency.  They  argued  that 
there  are  two  opposing  mechanisms  at  play  vis-a-vis  the  Reynolds  number.  While  a 
higher  Reynolds  number  is  associated  with  a higher  free  stream  velocity  and  a 
correspondingly  higher  film  heat  transfer  coefficient,  higher  Reynolds  numbers  also 
result  in  a higher  frost-air  interface  temperature  (see  Figure  7-9).  For  a constant  free- 
stream  air  temperature,  the  temperature  driving  force  between  the  free  stream  air  and  the 
cylinder  surface  is  reduced  as  the  Reynolds  number  is  increased,  resulting  in  a lower  heat 
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transfer  rate.  This  seems  to  reduce  the  overall  positive  effect  of  a higher  Reynolds 
number  from  a heat  transfer  rate  vantage  point.  Sherif  et  al.  (1993)  also  argued  that  the 
reduced  temperature  difference  driving  force  is  associated  with  a reduced  humidity  ratio 
difference  driving  force,  which  results  in  a reduction  in  the  frost  deposition  rate.  This 
contributes  to  the  tapering  off  of  the  frost  thickness-time  curves  shown  in  Figure  7-8  as 
time  increases.  An  examination  of  this  phenomenon  for  supersaturated  air  cases  indicates 
that  the  effect  in  question  is  less  pronounced  than  in  subsaturated  air  cases.  This  may  be 
an  intuitive  result  given  the  fact  that  frost  formation  in  supersaturated  air  is  driven  by 
both  diffusion  and  convection  and  is  thus  less  susceptible  to  being  affected  by  a single 
variable  such  as  the  Reynolds  number. 

In  order  to  verify  the  validity  of  the  model  described  here,  results  form  Padki  et 
al.  (1989)  and  Chung  and  Algren  (1958a,  b)  for  cylinders  in  humid  air  cross  flow  under 
subsaturated  conditions  are  utilized.  The  cylinder  employed  in  both  studies  was  25.4  mm 
in  outside  diameter.  Variation  of  the  frost  surface  temperature  with  time  for  two 
Reynolds  numbers  is  given  in  Figure  7-9  for  both  subsaturated  and  supersaturated  air. 
Results  of  the  present  model  are  compared  with  those  of  Padki  et  al.  (1989)  for  the 
subsaturated  cases.  From  the  figure,  one  can  easily  see  that  results  of  both  models  agree 
very  well.  For  the  same  Reynolds  number,  the  frost  surface  temperature  is  significantly 
higher  when  supersaturated  air  is  used.  This  can  be  explained  using  the  results  of  Figure 
7-8.  A larger  frost  thickness  necessarily  contributes  to  increasing  the  frost  surface 


temperature. 
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Figure  7-9.  Variation  of  the  frost  surface  temperature  with  time  for  different  Reynolds 
numbers:  a)  Re=2592:  — present  work  (subsaturated  conditions),  • Padki  et 
al.  (1989),  and  + present  work  (supersaturated  conditions),  b)  Re=9965:  — 
present  work  (subsaturated  conditions),  o Padki  et  al.  (1989),  and  x 
(supersaturated  conditions) 

A similar  comparison  with  the  data  of  Padki  et  al.  (1989)  and  Chung  and  Algren 
(1958a,  b)  can  be  performed  employing  the  total  frost  mass  accumulated  on  the  cylinder 
surface.  Figure  7-10  shows  how  the  total  mass  changes  for  two  Reynolds  numbers  under 
both  subsaturated  and  supersaturated  air.  In  this  figure,  results  from  Padki  et  al.  (1989) 
and  Chung  and  Algren  (1958a,  b)  were  included  for  the  subsaturated  case.  Form  the 
figure,  one  can  see  that  the  frost  mass  computed  based  on  the  present  model  is  slightly 
larger  than  those  of  Chung  and  Algren  (1958a,  b),  but  smaller  than  those  of  Padki  et  al. 
(1989)  for  the  subsaturated  case.  For  the  same  Reynolds  number,  the  amount  of  frost 
deposited  on  the  cylinder  increases  when  supersaturated  air  exists.  The  gap  between  the 
subsaturated  and  supersaturated  cases  widens  in  so  far  as  the  frost  mass  at  higher 
Reynolds  numbers  is  concerned. 
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Figure  7-10.  Variation  of  the  mass  of  frost  with  time  for  different  Reynolds  numbers:  a) 
Re=2592:  — present  work  (subsaturated  conditions),  ▲ Chung  and  Algren 
(1958),  • Padki  et  al.  (1989),  and  — present  work  (supersaturated 

conditions),  b)  Re=9965: present  work  (subsaturated  conditions), 

□ Padki  et  al.  (1989),  o Chung  and  Algren  (1958),  and  x (supersaturated 
conditions) 

Conclusions 

A semi-empirical,  quasi-steady  state  calculation  procedure  capable  of  computing 
the  frost  thickness  and  surface  temperature  along  with  the  mass  of  frost  deposited  on 
surfaces  in  the  vicinity  of  supersaturated  air  was  presented.  Model  results  were  compared 
with  existing  experimental  data  for  the  subsaturated  cases  and  were  found  to  agree  very 
well.  This  model  intended  to  illustrate  a relatively  simple  but  accurate  methodology  to 
predict  frost  formation  and  heat  transfer  rates  on  cold  surfaces.  It  should  open  doors  for 
studying  coil  frosting  problems  in  the  supersaturated  zone  of  the  psychrometric  chart  for 
more  complex  geometries. 


CHAPTER  8 

THEORY  OF  FROST-FORMATION  MECHANISMS 
Formation  Mechanism  in  the  Subsaturated  Zone 
While  the  mechanism  of  frost  formation  remains  essentially  the  same  in  both  the 
subsaturated  and  supersaturated  zones,  the  formation  mechanism  under  the  latter 
condition  is  much  more  complex.  The  additional  complication  arises  mainly  due  the 
presence  of  supersaturated  air  with  airborne  ice  crystals  that  tend  to  directly  deposit  on 
neighboring  surfaces,  including  the  freezer  walls  and  the  coils.  The  traditional  frost- 
formation  problem  is  complicated  in  its  own  right  beyond  the  traditional  moving 
boundary  problem  with  the  simultaneous  transfer  of  heat  and  mass.  This  is  primarily  due 
to  the  fact  that  the  formation  process  cannot  solely  be  described  by  determining  the 
location  of  the  boundary.  Other  complicating  factors  include  the  fact  that  frost  properties 
as  well  as  the  frost-air  interface  temperature  are  both  spatially  and  temporally  dependent. 
The  change  in  the  frost-air  interface  temperature  is  primarily  attributed  to  the  change  in 
the  frost  thermal  resistance.  As  the  interface  temperature  changes,  the  partial  pressure  of 
water  vapor  in  the  vicinity  of  the  frost  surface  also  changes,  thus  causing  a change  in  the 
structure  of  the  thermal  and  diffusion  boundary  layers.  Such  changes  eventually  affect 
both  the  heat  transfer  and  frost  deposition  rates.  In  cases  where  the  surface  temperature 
of  the  frost  is  in  the  vicinity  of  the  freezing  temperature,  additional  complications  may 
arise  due  to  the  occurrence  of  repeated  cycles  of  melting  and  re-freezing  which  tend  to 
change  the  structure  of  the  frost  layer  and  consequently  change  the  frost  density  and 
thermal  conductivity  with  both  time  and  position  without  appreciably  increasing  the  frost 
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thickness  (Sherif  1991).  This  is  more  likely  to  occur  in  high-humidity  and/or  high- 
temperature  environments  and  is  likely  to  result  in  an  uneven  distribution  of  the  frost 
density,  at  least  in  the  initial  phases  of  the  melting  and  re-freezing  processes.  The 
unevenness  in  the  frost  density  profile  is  primarily  attributed  to  the  fact  that  the  frost 
deposited  will  partially  melt  and  seep  through  the  pores  of  the  porous  frost  layer, 
eventually  reaching  the  metal  surface  where  it  freezes  into  ice,  and  partially  stays  as  frost, 
thus  contributing  to  an  increase  in  the  frost  thickness  (Sherif  et  al.  1993).  Repeated 
occurrence  of  this  phenomenon  will  cause  the  frost  density  to  be  larger  near  the  metal 
surface  relative  to  that  near  the  upper  frost  layers  and  can  eventually  transform  the  whole 
layer  into  an  ice  block  after  very  long  deposition  times.  The  transformation  into  a block 
of  ice  can  usually  be  attributed  to  the  fact  that  the  frost  layer  loses  most  of  its  porosity  as 
the  melt  continues  to  freeze  while  seeping  downward  through  the  porous  layer.  With  the 
pores  mostly  closed,  the  melt  will  be  forced  to  collect  at  layers  that  are  successively 
closer  to  the  frost-air  interface  and  away  from  the  metal  surface,  thus  causing  the  would- 
be  ice  layer  to  form  near  the  frost-air  interface.  This  trend,  if  continued,  would  contribute 
to  creating  a more  uniform  density  distribution  in  the  evolving  ice  layer  after  a significant 
amount  of  time  has  elapsed  since  the  start  of  the  melting  and  re-freezing  phenomenon. 

While  the  above  issues  seem  to  be  somewhat  intuitive,  other  issues  remain 
partially  unresolved  and  perhaps  somewhat  ambiguous.  One  of  these  issues  has  to  do 
with  the  mechanism  that  determines  how  much  of  the  water  vapor  transforming  into  frost 
actually  deposits  on  the  upper  frost  layers  and  how  much  of  it  diffuses  inward  into  the 
pores  of  the  porous  frost  layer.  White  and  Cremers  (1974)  report  that  approximately  half 
the  amount  of  water  vapor  transferred  into  the  frost  layer  diffuses  inward,  condenses,  and 
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increases  the  density  of  the  layer,  while  the  other  half  deposits  at  the  surface  and 
contributes  directly  to  increasing  the  thickness  without  appreciably  increasing  the 
density.  This  diffusion  would  normally  take  place  long  before  any  melting  and  re- 
freezing had  started.  While  this  theory  seems  to  be  a viable  one,  the  exact  fraction  of 
water  vapor  depositing  onto  versus  that  diffusing  into  the  frost  layer  is  perhaps  a bit  more 
contentious. 

One  of  the  issues  that  may  need  further  investigation  has  to  do  with  the  nature  of 
the  structure  of  the  frost  layer.  In  the  case  of  snowflakes,  for  example,  recent  studies 
such  as  that  of  Besant  (1999)  report  that  there  are  80  different  crystal  configurations  each 
of  which  start  and  grow  at  dust  particle  nucleation  sites  in  air  under  slightly  different 
temperature  and  humidity  conditions.  Besant  (1999)  explains  that  frost  growth  is  similar 
to  that  of  snow  and  is  thus  initiated  at  nucleation  sites  as  well.  When  water  surface 
tension  forces  are  relatively  small  such  as  the  case  on  aluminum  or  copper  surfaces, 
Besant  (1999)  argues  that  microscopic  droplets  of  water  would  form  first  at  nucleation 
points  on  the  surfaces.  Those  droplets  would  grow  in  size,  which  differ  depending  on 
temperature,  before  the  water  would  freeze.  Water  freezing  would  cause  frost  crystals  to 
form  in  the  vicinity  of  the  location  where  freezing  took  place  last  on  the  water  droplet 
surface  farthest  from  the  solid  surface.  This  issue  seems  not  to  have  been  researched  to  a 
significant  enough  degree  and  some  research  should  prove  useful  in  providing  additional 
insight  into  the  mechanisms  at  play  as  far  as  frost  crystal  formation  patterns  are 
concerned.  These  formation  patterns  are  necessarily  more  mysterious  when  the 
conditions  prevailing  in  the  freezer  correspond  to  the  supersaturated  regime. 
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Other  unresolved  issues  have  to  do  with  the  mechanism  that  causes  the  heat 
transfer  performance  degradation  in  frosted  coils.  It  was  believed  all  along  that  this 
degradation  is  caused  by  the  ever-increasing  thermal  resistance  of  the  frost  layer 
(Stoecker  1960).  While  this  may  be  the  main  reason  in  the  case  of  plain  surfaces,  it  may 
not  be  the  ultimate  cause  when  it  comes  to  finned  surfaces  according  to  Barrow  (1985) 
who  argued  that  the  decrease  in  the  air-side  coefficient  of  heat  transfer  is  a result  of 
geometrical  blockage  by  the  frost  layer  as  opposed  to  an  increase  in  the  frost  resistance. 

While  the  above  observations  represent  a reasonable  description  of  the  current 
state  of  knowledge  in  so  far  as  the  frost- formation  mechanism  in  the  subsaturated  zone  is 
concerned,  a significant  amount  of  research  is  needed  in  order  to  provide  a comparable 
amount  of  knowledge  concerning  the  mechanism  of  formation  in  the  supersaturated  zone. 
Prior  to  the  conduct  of  this  research  program,  information  pertaining  to  the  frost- 
formation  mechanism  in  the  supersaturated  zone  was  sketchy.  Here  we  attempt  to  devise 
a viable  theory  based  on  experimental  evidence  in  the  laboratory  as  well  as  field 
observations  (Smith  1992).  This  will  be  the  subject  of  the  following  section. 

Formation  Mechanism  in  the  Supersaturated  Zone 
Frost  formation  in  the  supersaturated  zone  was  probably  first  reported  by  Smith 
(1992)  who  presented  some  psychrometric  calculations  for  freezer  coils  to  avoid  coil 
operation  in  the  supersaturated  zone  (Smith  1992)  and  discussed  issues  pertaining  to 
latent  heat,  equipment-related  loads  and  applied  psychrometrics  in  typical  industrial 
freezers  (Smith  1992).  Supersaturated  conditions  in  freezers  usually  exist  under  two 
most  common  scenarios.  The  first  scenario  pertains  to  mixing  of  warm  humid  air  from 
outside  the  freezer  with  cold  dry  air  from  within  the  freezer.  This  scenario  is  likely  to 
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take  place  when  the  freezer  door  is  opened  and  will  result  in  the  mixing  process  line 
between  the  two  conditions  crossing  the  saturation  curve  on  the  psychrometric  chart  into 
the  supersaturated  zone.  This  scenario  was  investigated  by  Sherif  et  al.  (1997)  and  Al- 
Mutawa  (1997).  This  scenario  is  very  common  in  walk-in  freezers  where  air  mixing  is 
likely  near  the  freezer  door. 

In  a typical  interaction  between  the  freezer  and  the  anteroom,  an  air  exchange 
takes  place  such  that  the  cold  air  from  the  freezer  leaves  at  the  bottom  and  mixes  with  the 
anteroom  air  creating  supersaturated  conditions  outside  the  freezer  (Sherif  et  al.  1997). 
This  causes  suspended  liquid  water  droplets  to  exist  in  the  air,  which  result  in  fog 
formation  at  the  lower  anteroom  air  layers.  Concurrently,  the  anteroom  air  at  the  upper 
layers  rushes  inside  the  freezer  to  replace  the  exfiltrated  air  and  mixes  with  the  freezer  air 
in  the  upper  layers  of  the  freezer.  This  creates  supersaturated  conditions  inside  the 
freezer  with  airborne  ice  crystals  and  an  accelerated  formation  of  snow-like  frost  on  the 
freezer  walls  and  coils.  Snow-like  frost  is  different  than  the  traditional  frost  that  forms 
under  subsaturated  conditions  in  that  it  has  less  favorable  heat  transfer  characteristics  and 
is  denser.  The  latter  characteristic  causes  snow-like  frost  to  require  a larger  amount  of 
energy  to  defrost  since  the  process  takes  longer  to  complete  as  reported  by  Sherif  et  al. 
(1997),  Al-Mutawa  (1997),  Al-Mutawa  et  al.  (1998a,  1998b,  1998c,  1998d)  , and  Al- 
Mutawa  and  Sherif  (1998).  This  makes  the  formation  of  snow-like  frost  even  more 
undesirable  because  of  the  double  energy  penalty  associated  with  its  presence  and 
removal. 

The  second  scenario  under  which  supersaturated  air  exists  is  associated  with  air- 
cooling  processes  inside  the  freezer.  If  the  difference  between  the  air  entering 
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temperature  (or  temperature  upstream  of  the  coil)  and  the  coil  refrigerant  temperature 
inside  the  freezer  is  large  enough  and  if  a significant  latent  load  exists  in  the  freezer,  the 
cooling  process  line  is  likely  to  cross  into  the  supersaturated  zone  on  the  psychrometric 
chart.  The  temperature  difference  mentioned  above  is  referred  to  as  the  coil  TD  in 
industrial  refrigeration  circles  and  is  frequently  used  as  a measure  of  how  well  the  coil  is 
operating.  As  mentioned  above,  the  crossing  of  the  cooling  process  line  into  the 
supersaturated  zone  is  frequently  accompanied  by  the  presence  of  a significant  latent  load 
in  the  freezer. 

As  reported  earlier,  the  frost-formation  process  in  the  supersaturated  zone  is 
probably  driven  by  two  distinct  mechanisms.  The  first  mechanism  is  the  traditional 
diffusion  mechanism  as  is  the  case  in  the  subsaturated  zone.  In  this  case,  the  water  vapor 
concentration  (or  the  humidity  ratio)  driving  force  between  the  free  stream  and  the  air  in 
the  vicinity  of  the  cold  surface  acts  as  the  main  driving  mechanism  for  the  transfer  of 
water  vapor  from  the  free  stream  to  the  cold  surface  in  the  form  of  solid  frost  (bypassing 
the  liquid  phase).  The  water  vapor  difference  is  predominantly  the  result  of  the 
temperature  difference  that  exists  between  the  free  stream  and  the  cold  surface.  Since  the 
cold  surface  is  often  cooler  than  the  free  stream,  the  mass  transfer  driving  force  is  often 
towards  the  cold  surface.  However,  the  mass  transfer  driving  force  is  not  solely  a 
function  of  the  temperature  difference.  Frost  will  form  if  the  free  stream  air  in  the  freezer 
has  a larger  humidity  ratio  than  the  air  present  in  the  vicinity  of  the  cold  surface  provided 
that  the  temperature  of  the  latter  is  less  than  the  dew-point  temperature  of  water  vapor  in 
air  and  also  below  the  freezing  point.  In  supersaturated  conditions,  the  above 
mechanisms  are  still  applicable  but  with  a greater  intensity  since  the  moisture  content  in 
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the  free  stream  is  significantly  larger  than  in  cases  involving  subsaturated  air.  Since  the 
excess  moisture  (beyond  saturation)  in  supersaturated  air  is  incapable  of  existing  in  the 
vapor  phase,  it  would  have  to  exist  in  one  of  two  possible  formations;  either  as  suspended 
liquid  droplets  if  the  air  temperature  is  above  freezing  or  as  airborne  ice  crystals  if  the  air 
temperature  is  below  freezing.  The  presence  of  airborne  ice  crystals  contributes  to  an 
additional  amount  of  frost  accumulation  by  merely  depositing  on  any  cold  surface  they 
encounter,  including  the  coil  surface  and  the  freezer  walls  and  ceiling.  Thus,  while  one 
of  the  frost-formation  mechanisms  in  the  supersaturated  zone  is  identical  to  that  in  the 
subsaturated  zone,  a second  mechanism  driven  by  direct  precipitation  of  airborne  ice 
crystals  on  the  cold  surface  kicks  in.  This  mechanism  is  driven  by  convection  currents 
inside  the  freezer.  Airborne  ice  crystals  originate  as  frozen  suspended  liquid  water 
droplets  and,  as  a consequence,  are  fundamentally  different  in  their  crystalline  structure 
than  the  more  traditional  frost  crystals.  Frozen  liquid  water  droplets  have  a larger 
specific  heat  and  are  denser  and  thus  require  a larger  amount  of  time  to  defrost  relative  to 
the  more  traditional  frost  crystals.  Also,  because  those  airborne  ice  crystals  are  bigger  in 
size  than  the  traditional  frost  crystals,  they  tend  to  clog  the  passages  of  the  freezer  coil 
much  faster  and  with  a greater  intensity. 


CHAPTER  9 

SOFTWARE  FOR  THE  DESIGN  AND  OPERATIONAL  PHASES  OF 
REFRIGERATED  FACILITIES 

Objective 

The  objective  of  this  chapter  is  to  develop  two  pieces  of  software  that  can  be 
utilized  in  both  the  design  and  operational  phases  of  refrigerated  facilities. 

Software  for  Predicting  the  Onset  of  Frost  Formation 

This  software  will  be  based  on  the  analytical  model  described  in  Chapter  7.  It 
will  be  used  to  determine  the  frost  formation  and  heat  transfer  on  a cylinder  surface  under 
forced  convection.  Once  more  it  is  important  to  point  out  that  although  the  results 
obtained  using  the  software  are  for  a cylinder,  they  are  a good  approximation  for  tubes  in 
a freezer  coil.  This  software  performs  analysis  under  both  subsaturated  and 
supersaturated  freezer  conditions. 

User  manual.  In  this  section  the  functionality  of  the  simulation  software  will  be 
described.  This  will  be  done  by  presenting  the  different  possible  screens  that  will  appear 
during  the  operation  of  the  software.  When  the  software  is  started  it  will  show  the  Main 
Menu  screen  (Figure  9-1).  This  is  the  screen  that  will  allow  the  users  to  operate  the 
software.  Here  the  user  can  perform  several  tasks  such  as:  input  simulation  data,  see 
simulation  results  (numerical  and  graphical),  and  save  and  print  data  and  results  of  the 
simulation. 
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Figure  9-1.  Main  Menu  of  the  software 

0 If  the  “Data”  command  is  chosen,  the  software  will  show  a screen  where  the  user 
can  input  all  the  data  required  for  the  simulation  (Figure  9-2).  Here,  the  user  can  select 
the  type  of  simulation  to  be  performed:  a)  considering  average  properties  around  the 
cylinder,  and  b)  at  any  specific  angle  along  the  circumference  of  the  cylinder.  After  the 
user  inputs  the  data  of  the  simulation,  the  results  can  be  obtained  in  two  ways: 
numerically  or  graphically  by  selecting  the  command  “Results”  or  “Graphical  Results,” 
respectively. 

@ If  the  command  “Results”  is  selected  a new  screen  showing  the  results  will  appear 
(Figure  9-3).  The  results  that  are  displayed  on  this  screen  include:  frost  surface 
temperature,  mass  of  frost,  frost  thickness,  Reynolds  number,  and  the  convective  heat  and 
mass  transfer  coefficients. 

0 If  the  command  “Graphical  Results”  is  selected,  a new  screen  showing  a sub-menu 
(Figure  9-4)  will  appear.  In  this  screen  several  simulation  scenarios  can  be  selected.  This 
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includes  frost  thickness  vs.  time,  frost  surface  temperature  vs.  time,  mass  of  frost 
deposited  vs.  time,  and  the  convective  heat  and  mass  transfer  coefficients  vs.  time. 
Figure  9-5  shows  some  of  these  scenarios  on  the  screen. 


Type  of  Simulation 


Considering  Average  Conditions  Around  the 
Cylinder’s  Surface 

C At  Any  Specific  Point  on  the  Cylinder  Surface 


Data  for  the  Simulation 

Entering  Air  T emperature  (C) 
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Figure  9-2.  Screen  to  input  the  simulation  data 

0 If  the  command  “Air  Properties”  is  selected,  a new  screen  (Figure  9-6)  will  appear. 
This  screen  displays  air  properties  such  as  thermal  conductivity,  viscosity,  thermal 
diffusivity,  Prandlt  number,  and  the  specific  heat  for  any  temperature  that  the  user  inputs. 
@ If  the  command  “Data  and  Results  File”  is  selected,  a file  with  all  the  simulation 
information  such  as  data  and  results  will  be  saved  under  the  name  “Results.” 

@ This  command  can  be  used  to  print  the  data  and  results  of  the  simulation.  This  is  the 
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print  version  of  the  file  created  with  the  command  “Data  and  Results  File.” 

This  command  is  used  to  exit  the  software. 

If  this  command  is  selected  a new  screen  showing  the  information  about  the 
software  as  well  as  the  information  about  the  authors  will  be  displayed. 


Figure  9-3.  Screen  showing  the  numerical  results 


Figure  9-4.  Graph  menu  to  see  the  graphical  results 
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Figure  9-5.  Screen  showing  graphical  results:  A)  frost  surface  temperature  vs.  time  B) 
frost  thickness  vs.  time 
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In  Figure  9-5  (a)  and  (b)  if  the  command  “Go  Back”  is  selected  the  Graph  Menu 
screen  will  be  displayed  again.  Also,  the  graphs  show  in  Figure  9-5  can  be  printed  if  the 
“Print”  command  is  selected. 


Figure  9-6.  Screen  showing  air  properties 

Software  for  Simulating  Freezer  Operation 
The  second  software  will  be  used  to  guide  freezer  operators  in  avoiding 
supersaturated  operating  conditions.  The  idea  is  for  the  operators  to  input  the  system  data 
and  let  the  software  provide  them  with  information  as  to  where  the  system’s  operating 
point  is  on  the  psychrometric  chart. 

User  manual.  In  this  section  the  functionality  of  the  simulation  software  will  be 
described.  This  will  be  done  by  presenting  the  different  screens  that  appear  during  the 
execution  of  the  software.  When  the  software  is  started  it  will  show  the  Main  Menu 


screen  (Figure  9-7). 
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Figure  9-7.  Main  Menu  of  the  software 

0 If  the  command  “Data”  is  chosen,  the  software  will  display  a screen  where  the  user 
can  input  the  required  data  for  the  simulation  (Figure  9-8). 

@ If  the  command  “Results”  is  selected,  the  software  will  display  a screen  with  the 
results  (Figure  9-9).  As  can  be  seen  from  the  figure,  the  user  will  obtain  results 
numerically  and  graphically.  In  this  screen  the  solid  straight-line  represents  the  cooling 
and  dehumidification  process,  while  the  curved  line  represents  the  saturation  curve,  the 
dotted  straight-line  being  a reference  line.  The  interception  between  the  solid  and  dotted 
straight-lines  represents  the  air  leaving  conditions.  From  Figure  9-9(a)  it  can  be 
observed  that  the  system  is  operating  under  subsaturated  conditions.  This  is  shown  in  the 
box  named  “Note”  where  a message  is  presented.  In  Figure  9-9(b)  we  can  see  how  the 
software  displays  a warning  message  in  the  box  named  “Note,”  indicating  that  the  system 
is  operating  under  supersaturated  conditions.  The  critical  conditions  where  the  cooling 
and  dehumidification  straight-line  process  crosses  the  saturation  curve  are  shown  in  the 


box  named  “Critical  Conditions.” 
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Figure  9-8.  Simulation  data  screen 

(3)  If  the  command  “Properties  of  Supersaturated  Air”  is  selected  a screen  where  the 
properties  of  supersaturated  air  (specific  enthalpy,  specific  volume,  degree  of 
supersaturation,  and  supersaturated  humidity  ratio)  can  be  computed  is  displayed. 

(4)  This  command  is  used  to  exit  the  software. 
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Figure  9-9.  Some  results  screen  of  the  software:  A)  showing  operation  under  subsaturated 
conditions,  and  B)  showing  operation  under  supersaturated  conditions 


CHAPTER  10 
CONCLUSIONS 

Here  is  a summary  of  some  the  conclusions  of  this  investigation: 

• The  aforementioned  discussions  and  observations  underline  the  importance  of 
applying  psychrometric  theory  in  the  design  and  operation  of  freezers.  Without 
proper  understanding  of  the  theory,  one  can  easily  see  how  freezer  operation  can 
be  dominated  by  the  presence  of  unwanted  ice  crystal  formation,  snow-like  frost 
accumulation,  and  severe  degradation  in  the  coil  heat  transfer  performance. 

• In  this  investigation  the  existence  of  a demarcation  line  between  the  type  of  frost 
that  forms  under  subsaturated  conditions  and  that  frost  that  forms  under 
supersaturated  conditions  was  established.  With  the  use  of  a mathematical 
model,  the  location  where  the  transition  from  one  condition  to  the  other  occurs 
was  also  established. 

• Operating  freezers  under  supersaturated  conditions  is  most  costly  because  it 
increases  the  frequency  and  duration  of  defrosting  and  the  amount  of  defrosting 
heat  input  required. 

• The  use  of  heat  pipes  is  advantageous  for  freezer  operation  as  the  burden  of 
cooling  and  dehumidifying  the  air  is  redistributed  such  that  the  heat-pipe  coil 
performs  the  sensible  cooling,  while  the  dehumidifying  process  is  carried  out  by 
the  freezer  coil. 
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• The  use  of  partially-dampered  coils  was  shown  to  increase  the  defrost  efficiency 
by  18%,  while  the  use  of  fully-dampered  coils  increases  the  efficiency  by  43%. 

• The  use  of  dampered  coils  prevents  evaporated  water  and  sublimated  frost  from 
going  back  into  the  freezer  space,  thus  helping  reduce  the  defrosting  frequency  of 
other  coils  that  may  be  operating  in  the  refrigeration  mode  in  the  vicinity  of  the 
coil  being  defrosted. 

• Increasing  the  coil-face  velocity  during  refrigeration  increases  the  defrost 
efficiency.  This  may  be  partially  explained  by  the  more  cohesive  nature  of  the 
frost  forming  at  higher  face  velocities  and  the  subsequent  improvements  in  melt 
collection  and  the  corresponding  reduction  in  evaporated  water  and/or  sublimated 
frost  into  the  freezer  space. 

• A semi-empirical,  quasi-steady  state  calculation  procedure  capable  of  computing 
the  frost  thickness  and  surface  temperature  along  with  the  mass  of  frost  deposited 
on  surfaces  in  the  vicinity  of  supersaturated  air  was  presented.  Model  results 
were  compared  with  existing  experimental  data  for  the  subsaturated  cases  and 
were  found  to  agree  very  well.  This  model  should  open  doors  for  studying  coil 
frosting  problems  in  the  supersaturated  zone  of  the  psychrometric  chart  for  more 
complex  geometries. 

• Psychrometric  data  for  supersaturated  air  have  been  developed  and  presented  in 
the  form  of  psychrometric  charts  for  degrees  of  saturation  from  100%  to  200%. 
This  method  extends  the  current  psychrometric  data,  which  are  presently  only 
available  in  the  subsaturated  zone  of  the  psychrometric  chart.  The  properties 
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developed  in  this  investigation  should  enable  the  refrigeration  engineer  and 
system  designer  to  better  predict  the  performance  of  the  system  as  well  as  design 
a more  optimal  system. 

• Two  pieces  of  software  whose  primary  purpose  is  to  guide  refrigeration  system 
design  engineers  and  freezer  operators  in  their  mission  were  developed. 


APPENDIX 

UNCERTAINTY  ANALYSIS 


The  method  used  for  the  uncertainty  calculations  is  the  one  described  by  Kline 
and  McClintock  (1953).  Suppose  that  a result  “R”  is  a function  of  a number  of 
independent  variables  XI,  X2  , X3  , Xn. 

R = R(X{,X2,Xr...,X J (A-l) 

And  each  of  those  variables  has  an  uncertainty^  ,^x  ,—,%x  )•  The  total 

uncertainty  in  the  result,  ”,  is  found  through  the  following  equation: 
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In  this  investigation  properties  of  R-22  were  determined  using  software  developed 
by  Gallagher  et  al.  (1993).  The  properties  of  the  water  were  determined  using  software 
developed  by  Intellipro  (1998).  The  uncertainties  of  these  properties  were  calculated 
using  the  same  method  describe  above.  For  example,  the  enthalpy  is  a function  of 
temperature  and  pressure;  i = i(t,P),  therefore  the  uncertainties  of  the  enthalpy  were 


determined  using  Equation  A-2. 
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In  this  equation  the  terms  di/dt  and  di/dP  were  determined  numerically  using 


the  software,  while  the  uncertainty  of  the  temperature  (£,)  and  pressure  (<^p)  were  taken 


from  Table  A-l.  The  uncertainties  of  these  properties  are  shown  in  Table  A-2. 
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Table  A- 1 . Uncertainty  of  experimental  measurements 


Quantity  Instrumentation  Method 

Uncertainty 

Temperature  Copper-constantan  (Type  T)  thermocouples 

Frost  surface  FLIR  Systems  thermacam  PM  695  infrared 

temperature  camera 

Relative  humidity  Mamac  HU-224  humidity  transducer. 

Pressure  Mamac  PR-274  and  PR-264  pressure 

transducers 

Coil-face  velocity  System  controls  and  instrumentation 

M4000E  digital  HVAC  analyzer 

Mass  of  the  melt  Arlyne  scale  Co.  610L  precision  digital 

scale 

£ = ±0.2°C 
^ = ±2°C 

%RH  = -2% 

£p=  ±1%  of  the 
measured  pressure 
<^u  = ±2%  of  the 
measured  velocity 
£m  = ±0.005  kg 

5 mmeU  b 

Table  A. -2  Uncertainty  of  the  refrigerant  and  water  properties 

Quantity 

Uncertainty 

Enthalpy  of  refrigerant,  iR 

t.  = ±2% 

~‘r 

Density  of  refrigerant,  pR 

^,-±2% 

Dynamic  viscosity  of  refrigerant,  fJ.R 

Enthalpy  of  liquid  water,  iliquid 

^“±2% 

Enthalpy  of  ice,  iice 

4»=±2% 

Table  A-3.  Uncertainty  of  the  results 

Quantity 

Uncertainty 

Refrigerant  volumetric  flow  rate,  VR 

^ = ±2% 

Refrigerant  mean  velocity,  um  R 

Reynolds  number,  Re 

^=±3-5% 

Discharge  coefficient,  Cd 

#c,=±1% 

Flow  coefficient  of  the  orifice  meter,  Korijice 

tK  = ±1% 

Refrigerant  vapor  mass  flow  rate,  mR  vap 

t.  = ±2.5% 

Z’mR,vap 

Defrosting  heat  input  rate,  Qdej  in 

Zn  =±4% 

xZdef  ,in 

Defrosting  heat  input,  Qde f in 

to*,-** 

Rate  of  heat  remove  with  the  FCU-T  melt,  Qmelt 

Heat  remove  with  the  FCU-T  melt,  QmeU 

N© 

ox 

+1 

II 

vj 

Defrost  efficiency,  T]dcf 

L = ±4% 

' 7Tldef 

Average  coil  performance  load  rate,  ACPLR 

^ ACPLR  ~ ±4% 

Fan  defrost  contribution,  FDC 

% FDC  = -4% 

159 


Table  A-3.  Continued 


Quantity 

Uncertainty 

Air  volumetric  flow  rate,  Va 

^a  = ±2% 

Air  velocity,  ua 

Air  density,  pa 

Saturation  pressure,  Ps 

£ps  = ±2% 

Partial  pressure  of  water  vapor,  Pw 

#Pw  = ±3% 

Air  humidity  ratio,  Wa 

4w.=± 3% 

Air  enthalpy,  ia 

4 -±3% 

Table  A-4.  Function  and  partial  derivative  for  the  calculated  quantities 
Function  Partial  derivative 


Vr 

(m3/s) 


VR  = 4. 1 87657x1 0‘5  +(6.3786  lxl  0‘4Ap) 
- (2.6509  lxl  O'4  AP2) 

+ (7.27778x1 0‘5  AP3) 

-(l. 0033x1 0‘5AP4) 

+ (5.31 96x1 0-7  AP5 ) 


= [6.3786 lxl 0“4  -(5.30 182x1  O'4  AP 
SAP  L v 

+ (2 1.8333x1 0‘5  AP2) 

+ (4.0 132x1  O'5  AP3) 

+ (26.598x1  O'7  AP4)] 


um,R 

(m/s) 


“ orifice 


1 orifice 


= 0.0001393  m2 


1 


A 


orifice 


Re 


Re  = 


PRum,R^\ 

PR 


, D\  =13.386  mm 


dRe  urnD\  . dRe  D\Pr  . dRe  D lpRum 

dpR  nR  ’ dum  Pr  ’ dpR  p R 


cd  Cd  = [3.5959  + 0.03 12(p2-‘  -0.184v/8 

+ 0.039  - 0.0 1 58y/3  + 9 1 1 

1-/  Re0'75 

y/  = 0.5 


R orifice 


K 


orifice 


a/i- 


•;  i//=0.5 


dCd  68.7825y/2'5 

SPe  Pe175 


dK  Orifice  _ 1 

5C<t  V1-^4 


(m3/s) 


= w 


a 


™a=PaK 

(kg/s) 


y . 
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Table  A-4.  Continued 

Function 

^ R,vap  ril R.vap  ~ ^ orifice  ^ orifice  ‘P  'J  ~PPr P water 

(kg/s) 

g = 9.81  m/s2 

P water  = 1000  kg/™3 


Partial  derivative 


dm 


R.vap  _ q ■yQyj  ^ orifice  ^ orifice  W gP  water 


dpR 

dm 


V SPRP  water ^ 


R.vap  _ q J()l  ] K°r‘fice  ^ orifice  P gP  water  P R 


dAP 
dm 


-JgpRPwater^P 


dK. 


R'mP-  = A orifice  ll/2yj2gpRpH 

orifice 


,AP 


Qdef  in  Qdef , in  ™ R.vap  (fin  ‘ out 

(kW) 


Qdef  .in  Qdef. in  rh  R.vap  vin  ‘ out  r def 

(kJ) 


Qmelt 

(kW) 


Qmelt 


m 


melt  yliq 


(fliq  ‘ice ) 


1 def 


Qmelt  Qmelt  ^ melt  yliq  ‘ice 

(kJ) 


ride/ 


Vdef  = 


Qmelt 
Qdef , in 


ACPLR 

(kW) 


ACPLR=- 


(.Qtotal  ^ Q def  .in  Qmelt  ) 


1 total 


FDC 

(%) 


FDCSQfm  +Qdef-in  Qme“ \|Q0) 


dQdef.ii 

dm. 


(fin  ‘out  ) > 


1 R.vap 

dQdef.in  dQJef h 


di:. 


diin 


= m 


R.vap 


dQdef.i 


(fin  ‘ out  V” def  ’ 


dmRyap 

dQdefjn 

df def 

dQdef.in  _ dQdej  h 


™ R.vap  O', in  ‘ out  ) 


dilm 


di,. 


m R ,vap^  def 


dQmelt 

yliq  hce 

).  dQmel, 

» met, 1 

darnel, 

Tdef 

dr  def 

dQmelt 

dQmelt 

«. melt 

dinq 

^ice 

Tdef 

dQmelt 

(i  i 

1.  dQmell 

dQmelt 

dmme„ 

yliq  lice 

)>  m. 

dhiq 

^ice 

dridef 

1 

dridef 

Qmelt 

dQmelt 

Qdef , in 

dQdefjn 

Qdef , in 

dACPLR 

dACPLR  1 

dQdefin  dQmelt  ^ total 


dFDC  dFDC  100 

dQdefjn  dQmel,  Qnet 
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Table  A-4.  Continued 

Function 

Pa  P P 

(m3 /kg)  Pa  RaTa(\  + \.6mWa) 


Partial  derivative 

dpa  _ 1.6078/3 

dra  Rja2{\+\.emwtt) 
dpa  _ 1.6078 P 

dWa  RaTa{\  + \Mima)2 


Ps 

(kPa) 


Ps  = exp 


I«i(rr+«6 1°ge(7’) 


for  [263.15K<r^273.15K] 

P,  =exp  +b5  l°g e{T) 

o 

for  [273. 15K  <T<  473. 15K] 


dT 


= P, 


■ + fl2  — + 3^47*^  +4 h — — 


dP±_ 

dT 


- + b2  +2bJ  + 3bJ2  + 


bA 


Pw  PSRH 

(kPa)  w 100 


dPw  RH  DPW  _ Ps 

dPs  100  ’ dRH  100 


(kgw/kga) 


Wn  = 0.622- 


P-P,, 


dWa 

dPw 


= 0.622 


_J 

P-P 

1 1 W 


+ 


Pw 

{P-Pj 


ia  ia  = 1.006fa  + (2501  + 1.805/fl  )Wa 

(kJ/kg) 


* 

-^  = 1.006 + 1.8O50'a 

^ a 
8‘a 

dwa 


= 2501  + 1.805/ 
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